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@ _ Konfluente hypergeometrische Funktionen 


Zusammenfassender Bericht 


Von Francesco G. Tricomt, Turin!) 
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§ 1. Allgemeine Ubersicht 


Man kann wohl sagen, dass die meisten fiir die Anwendungen wichtigen 
_ speziellen Funktionen der Analysis nur zwei Klassen angehéren: elliptische 
_ Funktionen und hypergeometrische Funktionen mit Einschluss eines sehr 
_ wichtigen Grenzfalles dieser letzteren: die konfluenten hypergeometrischen Funk- 
“ tionen, wovon hier die Rede sein wird. Insbesondere stellen die Zylinderfunk- 
-tionen, die Laguerreschen (und Hermiteschen) Polynome, die unvollstandige 
_ Gammafunktion (mit ihren zahlreichen Sonderfallen: Fehler-Integral, Sinus- 
_ integral und verwandte Funktionen) usw. nichts anderes als Sonderfalle der 
_ konfluenten hypergeometrischen Funktionen dar, die so unter anderem die 
- Betrachtung von allen diesen Funktionen von einem einheitlichen Standpunkt 


= Ubrigens erlauben diese «konfluenten Funktionen», die lineare Differential- 


4 
j 
4 d®y dy nth i 
; 7 (dy +.by x) + (a1 + O14) Gy t (a2 t bax) 9 = 0 (1) 
$ 1) Istituto di Analisi Matematica dell’ Universita. 
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in allen Fallen geschlossen zu integrieren, was sehr interessant ist, weil diese 
Gleichung und ihre Transformierten in den verschiedenartigsten Anwendungen | 
vorkommen: zum Beispiel im Keplerschen Problem der Wellenmechanik. _ 

Trotz all diesem und der Tatsache, dass die grundlegende konfluente Funk-_ 
tion: die ganze Funktion 


co 


+1)...(a+n—1 is 
Ba, 63 2) = Fila; 0; 2) = ye Ey i 2) 


schon 1836 bei KuMMER [1]+) vorkommt, ist bis jetzt die Kenntnis der Theorie 

der allgemeinen konfluenten Funktionen nur verhaltnismdssig wenig verbreitet, 

besonders unter denen, welche sie am meisten brauchen: den angewandten 

Mathematikern. 
Mége dieser Bericht zur Beseitigung dieses Ubelstandes helfen! 


§ 2. Kummersche und Whittakersche Funktionen 


Die alteste der konfluenten Funktionen ist zweifelsohne die Kummersche 
Funktion (2), die einen Grenzfall der GauBschen hypergeometrischen Funktion 
F(a, b; c; 2) darstellt, weil 

HEA 1G SAS ES 
D(a, c; x) = lim F(a, b; ¢; >) (3) 
ist, wobei die zwei singularen Stellen x = b und x = co der Funktion rechts ins’ 
Unendliche zusammenjliessen. Daher der Name: konfluente hypergeometrische 
Funktion. _ 
Diese Funktion stellt eine erste Lésung der «konfluenten Gleichung»: 
d*y 


ad 
% agx + (0-4) Gay =0 (4) 


dar, einer Gleichung, welche als eine kanonische Form der Gleichung (1) be- 
trachtet werden kann). Noch mehr: solange der Parameter c keine ganze 
Zahl ist, besitzt (4) auch die von @ linear unabhangige zweite Lésung 


x1-° O(a —c+1,2—c¢; x); 
so dass dann das allgemeine Integral der konfluenten Gleichung in der Form 
y = Cy O(a, c; x) + Cy x1-° O(a —c + 1,2—c; x) (5) 
(C, und C, willkiirliche Konstanten) geschrieben werden kann. 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 273. 
2) Fiir die Reduktion von (1) (und anderen Differentialgleichungen) zur kanonischen Form (4), 
siehe [3] oder [4] oder [5]. Ahnliche Bemerkungen werden im folgenden meistens unterdriickt. 


Pear ee? 


a ON Pee eee 


_ 


he ee i eT er eee 


4 


4 


. 


Via Gere - tee. se A a teak —— . _- 
—¢ < ae x bis fxs re . ‘ ° as : 
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___ Trotzdem — abgesehen von der Schwierigkeit, dass c auch eine ganze Zahl 
sein kann — ist.es viel besser, die ganze Theorie auf einem Paar von Lésungen 


der Grundgleichung (anstatt auf der einzigen @-Funktion) zu bauen, wie es 
zum Beispiel E.T. WHITTAKER gemacht hat [2]. Nur dass dieser hervorragende 
- Forscher — anstatt der Gleichung (4) — die durch die Substitution 
erhaltene neue Gleichung (ohne Glied mit 2’) 
dz 1 x 1— 4p? 
Sr (4+tetght) <0 ° 


zugrunde legt; was mit sich zieht, dass — anstatt © — die nicht eindeutige 
Funktion A 
ME ipa eee O(u—x+5,2u+ is x) (8) 
in den Vordergrund gestellt wird und dass neben ihr eine zweite Losung W,, , 
auftritt, welche — solange 2 w keine ganze Zahl ist — durch die Formel 


a I'(—2 w) I'(2 p) 
Wan = P(1/2— w— x) M4) P(1/2 + w— x) M,, 4A) ©) 


berechnet werden kann. . 

Vielleicht hat gerade diese unndtige Mehrdeutigkeit der Grundfunktion 
M,,,, oder (noch mehr) die unbequeme Bezeichnung mit den Indizes x und yu 
(die manchmal recht komplizierte Ausdriicke werden kénnen) dazu beigetragen, 
dass die konfluenten Funktionen bis jetzt so wenig «volkstiimlich» geworden 
sind, trotz ihrer Beriicksichtigung in einem so verbreiteten und geschatzten 
Werke wie «Modern Analysis» von WHITTAKER-WATSON und ihrer grossen 
Wichtigkeit fiir die Anwendungen! 

Aus diesen und anderen Griinden habe ich seit einigen Jahren versucht, die 
elegante Einfachheit der Kummerschen @-Funktion mit den Vorteilen der 
«zweispurigen» Whittakerschen Theorie so zu verbinden, dass ich neben ® 
eine zweite Lésung (a, c; x) der konfluenten Gleichung (4) eingefiihrt habe, 
namlich eine solche Funktion, dass 

Waser P(u— “+ 52m tl; x) (10) 
ist; das heisst eine Funktion, welche in einem dhnlichen Verhialtnis zur Whitt- 
akerschen Funktion W,,, steht wie ® zu M,,,,. Auf diesen Grundlagen, das 
heisst auf dem Funktionenpaar ®, W, wird hier ebenso wie in den «Higher 
Transcendental functions» des « Bateman-Projects» [3] und in meinen Biichern [4] 
und [5] die ganze Theorie der konfluenten Funktionen aufgebaut. Eine neulich 
erschienene, wertvolle Monographie von H. BUCHHOLZ [6] beniitzt dagegen die 


¢ 


G 
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alten Whittakerschen Bezeichnungen. Wie dem auch sei, die Formeln (8) und 
(10) gestatten jederzeit einen leichten Ubergang von einem System von Be- 
zeichnungen zum andern. 


§ 3. Grundeigenschaft der @-Funktion 


Neben ihrer Potenzreihenentwicklung (2) und der Differentialgleichung (4) 
stellt die elegante Integraldarstellung 
1 


ih 
O(a, ¢; x) = Fare festte-3 (1 — pe-e-4 ds (0<Rea< Rec), (11) 
0 


die durch das klassische Laplacesche Verfahren leicht gewonnen wird, eine der 
Haupteigenschaften der -Funktion dar. Sind die angegebenen Giiltigkeitsbe- 
dingungen nicht erfiillt, dann ist nur das gewohnliche Integral durch ein 
Schleifenintegral auf dem iiblichen Weg zu ersetzen. 

Es gibt aber auch viele andere, teilweise wichtige Integraldarstellungen der 
@-Funktion, die auf den verschiedensten Wegen gewonnen werden kénnen. 
Wir wollen von ihnen hier nur die zwei folgenden erwahnen: 

&+400 ) 
O(a,c;—x) = rar f ae ees dz (Rex>0, 0<E< Re a), 


§—100 


eB 12 
P(a e: =) = ri fe oh Sant td Samed 0 iS 
0 


Li a 
=F Slt 1E._,()] (Rea>0, Res>0). | 


Die zweite Gleichung stellt die ®-Funktion als die Laplace-Transformierte des 
Produktes aus ¢¢—1 und der Zylinderfunktion 


—py/2 ie = bre 
E,() =#-"? J, (2Vt) = 2! Tv : n 2 1) n! (13) 


(fiir y = c — 1) dar, einer Funktion, welche manchmal der tiblichen Funktion 

J, vorzuziehen ist, weil sie auch im Falle einer nichtganzzahligen Ordnung » 

eine eindeutige Funktion bleibt. Diese Funktion besitzt lbrigens die einfachere 

Differentiationsformel Ej = —E,,, und ist eine Lésung der auch etwas ein- 

facheren (im Vergleich mit der von BEssEL) Differentialgleichung vom Typus (1): 
xy" + (1+) +y=0. 

Sehr wichtig sind ferner die Kummersche Formel 


D(a, c; x) = e* D(c — a, c; —x) (14) 
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"A 


und die sechs Differentiationsformeln 


d"D _ (4)n 

a D(atn,c+n;x), 

qd” 

da ee Mey (a), ee D(a 1, er x) , 

qd” 

dx” (2 Di (1a = c), 492"! Dia,c — 0; x) 

(15) 

a” —2x“ n (¢— @)n —2 . 

pane ©) = {=1) aa @ (a,c +n; x), 

qd” 

an (eae oo 42 8-1 D) (6 a),e-* xo 8! O(a +n, C3 %)., 

33 (e-* 4¢—-1 D) = (—1)" (1 — c), e~* 2° *-1 Dla — n,c— Nn; 4), 
wo n = 1, 2, 3, ... und (im allgemeinen) 

V(a+n) _ 

x (@= Ta) a(a+1)...(a+n-—1) 


Setzt man in diesen Differentiationsformeln ” = 1, so erhalt man aus je zwei 
Gleichungen durch Elimination von d®/dx als Rekursionsformeln fiir die ®- 


_ Funktion folgende 6 linearen Beziehungen zwischen ® = ®(4, ¢, x) und je zwei 


“ 
7 


ihrer Nachbarfunktionen: 
P(a+)=D(a+1,c;x), Pa—)=O(a—1,c¢; x); 
(16) 
P(c+) =DP(a,c+1;x), O(c—)=P(a,c—1; 4). 
Diese Formeln sind mit einem Zeichen versehen (was sich in der Praxis als 
sehr niitzlich erweist!), das schematisch zeigt, welche Funktionen (@ = o) 
jedesmal zusammengehéren. Seine Bedeutung ergibt sich von selbst, wenn 
man an eine Schematisierung der (reellen) Ebene denkt, wo a und c die ortho- 
gonalen Koordinaten sind. 


(c —a) D(a—) + (2a—c+x)PB—aP(a+) =0 e—o—e 

c(a + x) B— (c— a) x O(c+) —acP(a+) =0 sa 

c® —cP(a—) — x O(c+) =0 a 

(a—1+ 2) B+ (c— a) G(a—) — (¢— 1) B(e—) = 0 Fi af) 
(a—c+1)®—a@(at) + (c— 1) P(c—) = 0 jaa 

¢(c — 1) Bc -) —c (c—1+2%) G+ (c—a) xP(c+) = 9 i! 
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Zum Schluss wollen wir folgende wichtige Laplace-Transformierte 
1 
Q, [t?-1 B(a, c; )] = (0) s-” F(a, bs c: =) (Reb>0, Res>1) (18) 


erwahnen, zusammen mit ihrem bedeutenden Sonderfall b = c, in dem man 
Q,[t¢-1 O(a, c;t)] = Ic) s-* (kh —'s—4)=*,  (Ree> 0, Res 1) e aiee 
hat. 


§ 4. Definition und Grundeigenschaften der ¥-Funktion 


Die zweite Lésung Y(a, c; x) der konfluenten Gleichung (4) lasst sich besser 
als eine Laplace-Transformierte einer elementaren Funktion als durch eine 
Formel der Art von (9) definieren. Man kann namlich 


PG, 0, x) = ae Q, [¢¢-2 (1 + fe-4-2] 


1 (Rea>0, Rex>0) (19) 


ae fa-—1 (1 aE Ue ae ay : 


0 


setzen, wo — falls der Realteil von a nicht positiv ist - das gewéhnliche Integral 
durch ein Schleifenintegral: 


—ani or 
W a,c; x)= Pega fenmtee3 Lt apret a (a +1, 2,3,...) (20) 


zu ersetzen ist. 
Die somit in der Halbebene Re x > 0 erklarte Funktion WY wird leicht — 
durch Drehung des Integrationsweges um einen Winkel y (—2/2< p< a/2) - 
in der ganzen, entlang der negativen reellen Halbachse aufgeschnittenen kom- 
plexen x-Ebene verlingert. Der Schnitt stellt aber (im allgemeinen) eine Ver- 
zweigungslinie unserer Y-Funktion dar. 
Solange c keine ganze Zahl ist, wird die Gleichung (9) durch folgende, 
fundamentale Beziehung zwischen © und W ersetzt: 
eka) 


(a, ¢; x) = D(a, c; x) + - 


I'(c — 1) 
['(a—c¢+1) 


Ta) - *1~°@(a—c+1,2—c;x). 
(21) 


Diese Formel kann auch fiir die numerische Berechnung von ¥ benutzt wer- 
den’). Ist dagegen c= +1 (n=0, 1, 2, ...) eine positive?) ganze Zahl, dann 


1) Fir die numerische Berechnung von @ leistet die Definitionsreihe (2) gute Dienste. 


*) Der Fall einer nichtpositiven ganzen Zahl kann durch die folgende Formel (25) auf den 
Fall ¢ = » + 1 zuriickgefiihrt werden. 
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¥ 


_ sind die iiblichen Schwierigkeiten zu tiberwinden, die bei einer Differential- 
_gleichung vom Fuchsschen Typus wie (4) vorkommen, wenn der Unterschied 
zwischen zwei Wurzeln der bestimmenden Gleichung sich auf eine ganze Zahl 


; reduziert. Man sieht auf diese Weise ein, dass (21) durch folgende, weniger ein- 
4 fache Formel: 
V(a,n 4-154) = at {[o(a, n +1; x) logx 
+ 3 on Lola + m) — oll +m) — vl tm + ml ser} (22) 
ee ee at 


ersetzt werden muss, wobei p(z) = I’(z)/I'(z) die logarithmische Ableitung der 
Gammafunktion ist und die letzte Summe im Falle n = 0 wegzustreichen ist. 

Im grossen und ganzen sind die Grundeigenschaften der Y-Funktion denen 
~ yon @ ahnlich. Insbesondere gelten sechs ahnliche Differentiationsformeln wie 
_ die Gleichungen (15), wovon die folgenden drei eine Idee geben: 


ga ee re ee ee ee 


7 dx” Y ( 1) (4)n Y (a nm, e+ 1; %) , 
; % 
dx” (x ¥ ) (4)n (a Cc Ds x Y (a +n, C; x) ) (23) 


oil (e-* W) = (—1)"e-* Pla, e+; oh. 


ax” 


eee en eee ee 
' 


Daraus folgen — mit ahnlichen Bezeichnungen wie in (16) — die Rekursionsfor- 


meln fiir V: 
Wia—-)— (2a—c+2)P+ala—c+1) Pa+)=0 e—0-0 | 


(a+ x) W+a(c—a—1) Plat) —* Pe +) =0 oe 
(c— a) VW — xP(c+) + V(a—) =0 nell 
(a—1+2) B— Wa—-) + (@—c+1) Me-) = 9 4 (24) 
W — a P(at+) —Y(c—) =0 ie: 
(¢c—a—1) Pc—) — (¢-1 +4) P +2 Pct) =9 


und die Kummersche Formel findet ihr Gegenstiick in der wichtigen Beziehung 


Wa, c; x) = x1-° W(a—e + 1,2—0;4). (25) 
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Ferner entsprechen den Integraldarstellungen (12) und (18) folgende Glei- 
chungen fiir die Y-Funktion: 
E+1400 f ] 
rf a T(z) (a — 2) [(a—c+1-—2) Ede 
2700 


Ea, ¢; x) = T(a) T(a—c +1) 


(0<&< min [Re a, Re (a—c+ 1)]), 
1 Zeus? a—(c+ /2 
(a, ¢; 7) “fare See [Pree me dev 


(Rea>0, Rea> Rec +1, Res>0), ” 


_ (ob) P(b—¢+1) 
I(a+b—c+1) 


(26) 


eee WG cra)! 


x (Fo, b—e+1;a+b—-c+1;1~4), 
(Re b> 0, Rec< Re b+ 1, Re s> 0). 


Die in der zweiten Formel vorkommende K-Funktion ist die modtfizierte Zylin- 
derfunktion dritter Art nach Watsons Bezeichnungen. 


§ 5. Die verschiedenen Formen des allgemeinen Integrals 
der konfluenten Gleichung 


Durch die Funktion © und ¥ oder, viel besser, durch die Funktion Y und 
die Funktion 
a)» an 


. - = = 3 ( 
O*(aoarS TO D(a, ¢; x) = De I(c+n) x! Gn 


n=0 


(welche auch fiir nichtpositive ganze c ihre Bedeutung nicht verliert) kann 
man sofort folgende vier Partikularlésungen der konfluenten Gleichung (4) 
bestimmen: 


1 = OMG, C) x), 9g xe OtG CY, BS: 


Vg = F(a, t7 x), ty =e” Wea ae —%)} 
die sechs Paare von im allgemeinen linear unabhangigen Integralen unserer 
Gleichung erzeugen. 


Die entsprechenden Wronski-Determinanten haben alle die Form: 


We, = Kye? 2-80 s 123 Ay (28) 


a 


" 


Tae eee ee Ye AD mY 


a SAS 


NES t 
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wo die sechs Konstanten K,, folgende Werte besitzen: 
aes al 1 
Ky.=—sinac, Ky=—-3~, Ky = -=— 
uo A, ; “= L@) ae (1 — a)’ 
| (29) 
cede Y if 


Kys =— Kos == GE _ ern(c—a)i 


I'(c— a)’ I'(a—c+ 1)’ 


und e=sgn(3m x) ist. Die Ausnahmefalle der entsprechenden sechs ver- 


ITO © a 


_ schiedenen Formen des allgemeinen Integrals von (4) ergeben sich von selbst: 


Zum Beispiel die Form y= C,y,+ Cy 2, das heisst (im wesentlichen) die 
Form (5) des Integrals, darf nicht im Falle c= ganze Zahl, die Form 


y=C,4,+ Cz 43 nicht im Falle a= ganze Zahl <0 verwendet werden usw., 
_wiahrend die letzte Form: y = C, y3 + Cy ¥, immer Verwendung finden dart, 


- weil immer K,, + 0 ist. 


Ist eine von diesen sechs Integralformen, zum Beispiel die Form 
Cy Vr — Cy Vs 


festgesetzt, dann kénnen (von den obigen Ausnahmefallen abgesehen) die zwei 


, iibrigen Integrale als lineare Kombinationen von y, und y, ausgedrtickt werden. 


So kommt man auf eine Gruppe von zwolf wichtigen Formeln, wovon (9) ein 


_ Beispiel liefert, wahrend ein anderes Beispiel wie folgt iautet: 


é 


a 1 emct 
¥4~ Sina c ee Pur, P'(c — a) yal, (30) 


- wobei ¢, wie oben, mit sgn (jm x) zu identifizieren ist’). 


§ 6. Entwicklungen nach Zylinderfunktionen und 
Laguerreschen Polynomen 


Unter den vielen a priori méglichen Entwicklungen der konfluenten Funk- 
tionen nach vorbestimmten Funktionen (zum Beispiel trigonometrische Funk- 
tionen, orthogonale Polynome usw.) nimmt eine Entwicklung der ®-Funktion 
nach Zylinderfunktionen eine besondere Stellung ein. Es handelt sich namlich 
um eine Entwicklung, die — unter gewissen Umstanden — so rasch konvergiert, 
dass ihr erstes Glied schon fast den ganzen Betrag der Funktion liefert (cf. § 8). 

Diese vom Verfasser gefundene und in der ganzen komplexen x-Ebene kon- 
vergierende Entwicklung hat folgende Form: 


O* (a, Pie hn 0 Ln gg (4X) (x = — a). (31) 


n=0 


1) In [5], S. 76, findet man die vollstandige Tabelle der angedeuteten zwlf linearen Gleichungen 


vom Typus (30). 
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Sie kann — mit gleichzeitiger Bestimmung der Koeffizienten a,, — durch beid- 
seitige Anwendung der Laplace-Transformation (nach einer Multiplikation mit 
x°—1) und Verwendung der Formeln (12) und (18’) bewiesen werden. 


Man findet so, dass y 
an = 2-" Anl x, +) (32) 


ist, wo die A,(x, A) gewisse, vom Verfasser eingehend studierte [10] Polynome 
in x und 4 sind, welche durch die Rekursionsformel 
(w+ 1) Anygy = (1+ 24-1) A,_1-22%A,_2, (m =2, 3,...) 
und die «Anfangsbedingungen » 
Ag=1, A,=0, A,=A 
rasch bestimmt werden kénnen. Insbesondere findet man 


As=—3%, Ag=FAAt1), Ap=—— x(51 +43). 


Ersetzt man in (31) die eindeutige Zylinderfunktion E, durch die iibliche 
J,, dann ergibt sich 


D¥ (a, 6; x) = e712 (xe x) 0-92 sy A,(x, +) (=) tte (2 Vx). (33) 
n=0 


Mittels der Grundbeziehung (21) zwischen © und ¥ ist es leicht, eine ahn- 
liche Entwicklung auch fiir Y% zu bekommen, aber wir werden sie hier nicht 
angeben, weil bei den entsprechenden Rechnungen keine bemerkenswerte Ver- 
einfachung erzielt wird. Wir erwahnen hier lieber, ihrer grossen Einfachheit 
halber, folgende Entwicklung der Y-Funktion nach Laguerreschen Poly- 
nomen L®); 


eet, 
I(a) Pa, c; x)= S . ar ee (34) 
n=0 ; 


die zwar sehr schlecht konvergiert (wenn iiberhaupt?)), aber immer Abel- 
summierbar ist. Auch die Entwicklung (34) riihrt vom Verfasser her. 


§ 7. Asymptotisches Verhalten fiir x + oo oder x >0 


Das Hinzufiigen der zweiten Lésung ¥Y ist hauptsachlich darum so zweck- 
missig, weil Y ein ganz einfaches asymptotisches Verhalten fiir x -+ co zeigt, 
wahrend das von ® ziemlich verzwickt ist. 


1) Ein Fall, in dem die Reihe (34) sicher konvergiert, ist der Fall a < 1/2,0<el%*<w<oo. 


a 


Vol VI, 1955 Konfluente hypergeometrische Funktionen 267 
q Die Y-Funktion besitzt namlich die ganz einfache asymptotische Entwick- 
_ lung (im Sinne von Poincare) fiir z > oo: 

- 

q es 

; WwW a . aS 4 em la (4)n 

J (a, ¢; x)~ x 2a) . ae (35) 
: woraus sich insbesondere ergibt 

W(a,c; x) = x7 [1+ = 8— 2 + 02-9] (35" 
| 1; é ; ) 


- Um diese Entwicklung zu erhalten, braucht man nur das sogenannte Wat- 
~ sonsche Lemma bei den Definitionsformeln (19)—(20) anzuwenden. 

. Wahrend sich ®(a, c; x) fiir x > 0 ganz einfach verhilt, ist das Verhalten 
_ von Y komplizierter und muss aus den Grundbeziehungen (21)-(22) heraus- 
gefunden werden. Zum Beispiel findet man in dem etwas schwierigen Falle 
-~¢=1s0 die Formel 


—T'(a) Pa, 1; x) =logx + p(a) +2C+O(xlog|x|), (|¥| > 9) (36) 


- wo y dieselbe Bedeutung hat wie in (22) und C die Euler-Mascheronische Kon- 
_ stante bedeutet. 


§ 8. Asymptotisches Verhalten fiir grosse Parameterwerte 


Die obigen asymptotischen Darstellungen gentigen fiir viele Anwendungen 
nicht. Zum Beispiel in dem wichtigen Sonderfalle der Laguerreschen Polynome: 


L(x) = fe i % D(—n, a +1; %)5 (37) 


was am meisten interessiert, ist nicht das Verhalten fiir x oo, sondern das- 
jenige fiir 7 > oo, das heisst fiir a > — oo. 

In einigen Fallen (zum Beispiel im Falle c -> co, in welchem schon die Defi- 
nitionsreihe von ® einen «asymptotischen Charakter» besitzt) ist die Unter- 
suchung des Verhaltens der konfluenten Funktion fiir grosse Werte der Parameter 
leicht, aber im allgemeinen, und insbesondere wenn gleichzeitig |x| > 0o ist 
die Aufgabe nicht trivial. Ubrigens ist die Sache — der grossen Anzahl der 
a priori méglichen Falle wegen — noch nicht ganz abgeschlossen, und die vielen 
asymptotischen Darstellungen, die man schon zur Verfiigung hat, sind keines- 
wegs immer sehr einfach und einleuchtend?). 

Wie dem auch sei, so gibt es wenigstens einen Fall, in welchem sich das 
gesamte Bild des asymptotischen Verhaltens der konfluenten Funktion ® schon 


1) Fiir viele dieser Darstellungen siehe insbesondere [7]. 
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mit klaren Ziigen abzeichnet!), namlich der Fall, in dem a,c und x > 0 veell 
sind und (bei einem beschrankt bleibenden c) a gegen — oo strebt, wahrend x 
entweder beschrankt bleibt oder zusammen mit x = c/2 — a gegen Unendlich 
strebt. Derselbe Fall also, der sich bei den Laguerreschen Polynomen als der 
wichtigste erwiesen hat und welcher mit Verallgemeinerungen der bei den 
Laguerreschen Polynomen erprobten Methoden [8] behandelt werden kann [11]. 

Wir werden uns hier daher auf diesen Fall beschranken, welcher sich in 
folgende vier Unterfalle aufteilt : 

1. Die Variable x bleibt entweder beschrankt, oder sie ist O(x*) mit @ < 1/3. 

2. Man hat Ax Sx <4Bx, wo A und B zwei beliebige feste positive 
Zahlen kleiner als eins bedeuten. 

3. Es ist x = 4 % + O(n"), 

4. Man hat 4 B* x < x, wo B* irgendeine feste positive Zahl grésser als 
eins bedeutet. 

Der erste Unterfall ist rasch erledigt, weil die stets konvergente Entwicklung 
nach Zylinderfunktionen (33) dann auch einen «asymptotischen Charakter» in 
dem Sinne besitzt, dass der Rest der Reihe — falls sie mit einem Glied der 
Ordnungszahl n= 3m+1 (m=0,1,2,.. .) abgebrochen wird — immer von 
einer kleineren Gréssenordnung (fiir x > + oo) als das letzte betrachtete Glied 
ist. Insbesondere hat man fiir m = 0 


®*(a, c; x) = e7?(x x)a-9/2 [ie-1(2Vx x) +0 («-)|, o = min (1 = Oy 2) , 


(38) 


was eine Verbesserung einer klassischen Formel von PERRON [13] darstellt. 
Ferner darf man (38) (ohne Restglied) — falls x ziemlich gross im Vergleich mit 
c und x ist — auch als eine angendherte Formel fiir die Berechnung von @* 
betrachten, die manchmal nur ganz geringe Fehler verursacht. Mit anderen 
Worten: unter den angedeuteten Umstanden wird das erste Glied der Entwick- 
lung (33) fast den ganzen Betrag der Funktion ®* liefern. 

Im zweiten Unterfalle, das heisst in der «Oszillationsstrecken 0< x<4% 
(in welche in der Regel sdmtliche reelle Nullstellen der @-Funktion fallen), 
wendet man das klassische Sattelpunktverfahren an einer passenden, von (11) 
abgeleiteten Integraldarstellung der ®-Funktion an. Man erhalt so die asym- 
ptotische Entwicklung (im Sinne von POINCARE) 


ay pe Spy hla) s Mecos@aey sin (O* + 3 p 2/2) 
ONE Re rine Vecusnae Per ) sina op > (39) 


wobei 


(0, 77/2) % 


6 = arccos Vee » OF = x(2 0 —sin2 6) + (5 ae a) It 


1) Der Fall der Funktion Y wird durch (21) auf den der Funktion © zuriickgefiihrt. 
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ist, und die beiden ersten Koeffizienten A,(#) durch folgende Ausdriicke ge- 
_ geben sind: 


Ag(@)=1, A,(6) = +5 [7 cosec?@ + (3 c*— 6c +2) sin?6 — i]. 


_ Bricht man die Reihe nach dem ersten Glied ab, dann erhalt man unter Bertick- 
-sichtigung der Formel 


ii 2) ed — 6/2) 
F(¢—a) F(x+ e/2) 


ite TO (neo) 


eee ee ee paris ai 5 head) 


die ganz einfache, oftmals schon ausreichende asymptotische Darstellung: 


: ae 2 x cos 8)1-¢ 
; eo" OF (a6 4) = ( ) 
4 a sin 2 0 


| Im dritten Unterjalle — das heisst in der Nahe des «Wendepunktes» x = 4 x - 
_ wendet man das allgemeine asymptotische Verfahren fiir die Lésungen gewisser 
- linearer Differentialgleichungen an, das ich Fubini-Verfahren genannt habe[12]. 
Mit diesem Verfahren — welches eine Verbesserung des sogenannten Liouville- 
- Steckloff-Langerschen Verfahrens darstellt - kommt man auf eine neue asym- 
_ptotische Entwicklung, deren erstes Glied folgende einfache asymptotische 


- Darstellung 


[sinO* + O(x-})]. (40) 


31/3 


en? D*(a, c, x) = aed x)7/5-°[ A,(t) cosa + A,(t) sina + O(x-*”)] (41) 


liefert, wobei man 


a= 43 — (sx), t=O (1) 
- gesetzt hat und A,(¢) und A,(t) die beiden «Airyschen» Funktionen, das heisst 
gewisse Kombinationen von Zylinderfunktionen der Ordnungen =: 1/3, be- 
deuten. 
Im vierten Unterfalle endlich — das heisst in der «Monotoniestrecke» x > 4 x 
— arbeitet man mit ahnlichen Methoden wie im zweiten Unterfalle, aber nun ist 
der Unterschied zum Falle der Laguerreschen Polynome betrachtlicher. Die 
Endformeln enthalten hyperbolische, anstatt Kreisfunktionen, und das Gegen- 
stiick von (40) ist 


e-*? @*(a, c; x) = sinax nme — exp[x Gin2 &— 2 8][1 + O(u-)] (42) 
mw x Slt 


Coj # = ie 


mit 
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§ 9. Nullstellen der konfluenten Funktionen 


Im Reellen, das heisst fiir reelle a, c, kann man die wichtigen reellen Null- 
stellen von ® und Y wie auch die komplexen Nullstellen von V leicht abzahlen. 
Die einfachsten Ergebnisse auf diesem Gebiet betreffen die Anzahl P(a, c) der 
positiven1) Nullstellen von @, falls c >0 ist, und die Anzahl P*(a, c) der 
positiven Nullstellen von Y ohne Einschrankung fiir die Parameter. Man hat 
namlich, fiir c > 0, 


P(a,c)=0, (@21); Pla,c)=—[a], (a<0), 


wobei [a] den ganzen Teil der (negativen) Zahl a (das heisst die grésste ganze 
Zahl S a) bedeutet. Und fiir die positiven Nullstellen von Y hat man 


P*(a,c) = P*(a—¢+1,2—c)= Pla,c), (€=1). 


Was die Lage all dieser Nullstellen anbetrifft, so sind einige interessante 
Angaben in den folgenden Satzen enthalten: 

1. Die positiven Nullstellen von © (méglicherweise mit Ausnahme der 
grdssten) besitzen die obere Schranke 


2u+V4n%+0(2—0), 


welche, solange 4a < c?+ 1, kleiner als 4 x + 1 ist. 

2. Die r-te positive Nullstelle &, von O(a, c; x) ist (fiir c = const) eine ab- 
nehmende Funktion von x =c/2—a. Sie erscheint im Unendlichen fiir 
a=—(r—1). 

3. Es gilt ;2 
é etek 
o> 4%’ 

WO j-4,r die -te positive Nullstelle von J,_,(x) bedeutet und auch 


<— Tea, r = 
2% + V2 WA se eee 


’ 


solange die Quadratwurzel reell ausfallt. 
4, Das Produkt x &, strebt fiir x +00 gegen seinen Grenzwert 1/4 7?_, ,, 
und zwar abnehmend. . 
5. Samtliche Nullstellen von Y(a, c; x) besitzen die obere Schranke 


2u+V402—(c—1)?+14<Vibwt+1, 
solange die Grésse unter dem ersten Quadratwurzelzeichen positiv ausfallt. 


1) Dank der Kummerschen Beziehung kann man die negativen Nullstellen von © auf die 
positiven zuriickfiihren. 


—_ 


ayFr 


ie i 


AV Oe eee ee 


Ne ee ee ey Ee 


‘Vol. VI, 1955 Konfluente hypergeometrische Funktionen 27M 


6. Die Anzahl der reellen wnd komplexen Nullstellen von YW (a, c; x) auf der 
anzen (entlang der negativen reellen Halbachse aufgeschnittenen) x-Ebene ist 
wenigstens fiir reelles a und c) immer endlich. 

Ferner kénnen die positiven Nullstellen von ® mittels der Formeln des 
-vorigen Paragraphen fiir x > oo asymptotisch abgeschatzt werden. Insbeson- 
dere hat man 


wy by FP 6 oF 2 iL 20 % 3 
fr oe ee 


—-> 7Q 


Nu 
| 


wo 
& = 4 % cos? 0° 


und 6° die im Intervall (0, ) liegende Wurzel der transzendenten Gleichung 


26—sin26=(r+a— 7) = (44) 


bedeutet. Dabei ist angenommen, dass fiir x oo das 6° weder gegen Null 


_ noch gegen z strebt. Nach unserer Numerierung der Nullstellen ist diese Be- 


dingung fiir 7 = const nicht erfiillt. Sie wird aber erfiillt, wenn man die Null- 
stellen «von der Mitte ausgehend» numeriert. 


§ 10. Spezialfalle der konfluenten Funktionen 


Wie im Anfang angedeutet, ist einer der Hauptvorziige der Theorie der 
allgemeinen konfluenten Funktionen der, dass sie uns die Betrachtung der 
meisten speziellen Funktionen von einem einheitlichen Standpunkt aus ge- 
stattet. Die folgende kleine Tabelle gibt Auskunft iiber die wichtigen Spezial- 
falle der konfluenten Funktionen: 


reduzieren sich die hypergeometrischen 
konfluenten Funktionen auf 


Wenn die Parameter 
so sind, dass 


Zylinderfunktionen 


Laguerresche Polynome 
Unvollstandige Gammafunktionen 
Funktionen des parabolischen Zylinders 

@=¢,P =I (1 —4, +) 


Dabei ist immer ” = 0,1, 2,.... Ubrigens ist zu bedenken, dass — dank (14) 
und (25) — die Parameterpaare (a, c) und (c — a, c) oder (4, ¢) und (a—c+1, 
2 —c) fiir D bzw. V als nicht wesentlich verschieden zu betrachten sind. Und 
ebenso (mit einigen Vorbehalten) zwei solche Paare (a, c), (a’, c’), wenn a’ —a@ 
und c’ — c ganze Zahlen sind. Zum Beispiel fiihrt der Fall c = 3/2 gleich dem 
Falle c = 1/2 zu Funktionen des parabolischen Zylinders. 
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Die wichtigen unvollstandigen Gammafunktionen: 
x co 


y(a, x) Ee Pdi, Ia, x) =Tla)—ylaex) = fete dt | 


0 x 


oder (besser!) die modifizierte unvollstandige Gammafunktion: 


Y*(, *) = Fegy Vlo, #) = e-* O*(1, 1 +a; x), (45), 


welche eine ganze Funktion von x und « ist, lassen ihrerseits folgende bemer- 
kenswerte Unterfalle zu: 


a=0 Integralsinus und verwandte Funktionen 

«a= 1/2 Fehlerintegral, Fresnelsches Integral und verwandte Funktionen 
a=n+1 = Wien" 2 (x) e,(2) ==! x Ion 

“=—n Dn es 7 ee 


Ausser (37) und (45) sind die wichtigsten Verbindungsformeln: 


Js) = ae - (ZY) ei OCF +1429;272), | 


, 2 v 2 ,t(x-va 1 . 
HP(x) =— i /_ 2 2) +1/2 pi ( (> +914 2;-2i2), 


nr 
O(n +1,%+1; x) = —- . LS pena” "le, x), 


1 0” 


Yin+1,a+1; x) = aig secs ee [ena ™ F\ae,-2) | 


4 ik cy v/2 2 
Y(- Foi) = 27 e* Dia), 


(46) 


Ein «= | (1 — e~‘) st = I(0,*)+log x +C, 


x x 
Ein (iz) = Cin +i Six = /(1— cost) 2 +i [sine &, 
0 0 


1 ae 
Erfix = — Erf (i x) = fe” dt— xO, zi). 


0 
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Es versteht sich von selbst, dass die Grundeigenschaften der vorigen Funk- 
_tionen nichts anderes als Sonderfalle der Grundeigenschaften der allgemeinen 
konfluenten Funktionen sind. Es kommt aber eine Menge von neuen speziellen 
 Eigenschaften hinzu, wie zum Beispiel der wohlbekannte Fall der Zylinder- 
_ funktionen zeigt. Es handelt sich dabei um neue, teilweise sehr interessante 
4 Integraldarstellungen, Reihenentwicklungen usw. 

; So besteht zum Beispiel im Falle der unvollsténdigen Gammafunktion 
_ neben einer Entwicklung des Typus (33) auch eine wesentlich andere Entwick- 
7 lung nach Zylinderfunktionen, namlich die stets konvergente (vom Verfasser 
_ herriihrende) Entwicklung: 
: 
- 
j 


yon 2) =e" Salt) a Basalt) (a0) = Sor) A 
_ wok eine willkiirliche Konstante bedeutet. Ist dazu « = 1/2 (Fall des Fehler- 
integrals), dann kommt noch ein dritter Typus von Entwicklungen nach Zylin- 
derfunktionen in Frage, zum Beispiel eine verbliiffend einfache und rasch kon- 
vergente Entwicklung nach modifizierten Zylinderfunktionen erster Art J,, 
namlich: 


Y/2 Bet (Vz) = Lyla) — Lyle) — Taal) + Zyl) + Iggle) — = (48) 


Hier sind alle Koeffizienten gleich + 1, und zwar Andert sich das Vorzeichen 


_ nach je zwei Gliedern. Diese Entwicklung und die ahnliche fir Erfi (Vx) 
sowie auch (47) sind fiir numerische Rechnungen sehr gut geeignet. 
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Summary 


Short exposition of the essential features of the theory of the confluent 
hypergeometric functions on the basis of the function ,F, = ® of KuMMER, and of 
a second solution ¥ of the basic differential equation. 

Among other things these functions permit of a unitary consideration of most 
special functions of importance for the applications. 

This exposition, while it does not adduce proof, may enable the reader, to a 
certain extent, to use the confluent functions even if no more detailed work is 
available. 


(Eingegangen: 10. Mai 1954.) 


Instationare Randbedingung fiir die durch eine Querschnitts- 
anderung gebildete Ubergangszone in einer instationdren, 
eindimensionalen Gasstromung 


Von YIAN-NIAN CHEN, Winterthur?) 


In einer eindimensionalen Strémung kompressibler Gase wird sich jede 


Stérung als Druck- oder Verdiinnungswelle in beiden Richtungen fortpflanzen. 


Uber die eigentlichen Bewegungen der Gasteilchen iiberlagert sich somit 
eine Stérungserscheinung. Solche instationaren Strémungen kénnen graphisch 
nach dem sogenannten Charakteristiken-Differenzen-Verfahren berechnet wer- 
den. Im Randgebiet mit sprunghafter Querschnittsanderung ist aber eine insta- 
tiondre Behandlung der Strémung sehr schwierig, weshalb diese dort im allge- 
meinen unter Vernachlassigung der mit der Zeit veranderlichen Gréssen als 


1) Gebriider Sulzer AG. 
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quasistationar betrachtet wird. Diese Vereinfachung ist nur bei Stromungen 
mit Unterschallgeschwindigkeit gerechtfertigt, wenn das Zeitintervall des gra- 
phischen Verfahrens so gross gewahlt ist, dass ein gentigender Ausgleich der im 
Randgebiet hin- und herwandernden Wellen sichergestellt wird. Sollte hingegen 
ein einzelner Verlauf der Str6mung im Randgebiet berechnet werden, so muss 
sie unbedingt instationar behandelt werden. 

Bei Strémungen mit Uberschallgeschwindigkeit ist jedoch keine Welle im- 

stande, gegen die Stromrichtung zuriickzuwandern; somit ist auch ein solcher 
Wellenausgleich nicht méglich. Hier kommt nur eine instationare Betrachtung 
der Str6mung in Frage. 
: In der vorliegenden Arbeit wird eine instationare Uberschallstrémung in der 
- Diisenzone eines Rohres zunachst nach dem konventionellen Verfahren stationar 
behandelt, um dann mit dem Ergebnis der instationaren Betrachtung einen 
Vergleich zu erméglichen. Im letzten Abschnitt kommt das Problem der Uber- 
schallstro6mung in einem konischen Rohriibergang zur Besprechung. 


EE PO ee Tee Cee a 


1. Stro6mun¢gs¢gleichungen 


> 


| Die Strémung eines kompressiblen Gases in einem Rohr wird hier als ein- 
dimensional betrachtet; das heisst, seine Geschwindigkeit weise keine Kompo- 


- nenten und keine Veranderungen quer zur Rohrachse auf. Fir das in Figur 1 
wanadruc 
Norma/spannung: Norma/spannung : 
hk 
eV hed Pt Px OX FD (Wyt Wey IX) 
——— ——— > 
™, Wt Wy IX 
+ Ox OX 
Ff 
ax F +l, ox 
Figur 1 


Volumenelement der Rohrstr6mung. 


dargestellte Volumenelement gelten folgende Gleichungen: 
a) Kontinuitatsgleichung: 


fwodt —[fwot(fwondx]d= fax qd; (1) 


Massenzufuhr — Massenabfuhr = Massenanderung 
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b) Impulsgleichung : 


fdxq 5? = —(p),dx-+0(f% w,) dx +p fades y+ wifdxe (2) 


|w| 
Masse x = Resultierende Kraft an der Kontur in der Achsenrichtung aus 
Beschleunigung statischem innerer Reibung Wanddruck Wandreibung 
Druck in den 
Endquerschnitten 


c) Energiegleichung: 


fdxo-<(gut+ )=—@fwedx+n(fwefw) dx+A(fTyede + gah dx 


Anderung der Energie pro = Leistung der Normalspannungen -+ Warmestrom -+ Warme- 
Zeiteinheit (potentielle an den Endquerschnitten durch Leitung zufuhr 
Energie vernachlassigt) von aussen 


‘Hierin bedeuten: / Rohrquerschnitt; g Erdbeschleunigung; p statischer Druck; 
gq von aussen pro Volumen- und Zeiteinheit zugefiihrte Warme; T absolute 
Temperatur; ¢ Zeit; u innere Energie pro Gewichteinheit; w Eigengeschwin- 
digkeit des Gasteilchens; x Weg; 7 Zahigkeitskoeffizient; 2 Warmeleitzahl; 
€ Reibungskoeffizient und @ Dichte. 

Wenn die Glieder fiir Warmeleitung und innere Reibung nicht beriicksich- 
tigt werden, lassen sich die Gleichungen mit folgendem Ergebnis lésen}) : 


a =w fir Lebenslinie, (4) 
a =w+t (=4)" fiir Stérlinien. (5) 


Sie entsprechen drei Scharen von Charakteristiken. 

a) Langs der Lebenslinie dx/dt = w in der (x, t)-Ebene bewegt sich ein Gas- 
teilchen mit der Eigengeschwindigkeit w. Das ist die eigentliche Bewegung des 
Gases. ; 

b) Langs der Stérlinie pflanzt sich die Stérung mit einer Ubergeschwindigkeit 
+ (x p/e)'” iiber das sich bewegende Gasteilchen hinweg fort. Die Uberge- 
schwindigkeit wird Schallgeschwindigkeit genannt und mit a bezeichnet. Die 
absolute Fortpflanzungsgeschwindigkeit der Stérung ist c= w + a. 

Ein Gas wird durch Druck- oder Verdiinnungswellen adiabatisch kompri- 
miert bzw. expandiert, wenn die Warmeleitung innerhalb der Gasmasse und 
die innere Reibung vernachlassigt werden. Die Schallgeschwindigkeit hangt in 
diesem Fall nur von der absoluten Temperatur des Gases ab. Aussere Reibung, 
dussere Warmezufuhr und Querschnittsdénderung haben keinen Einfluss auf diese 
Tatsache, soweit von deren Wirkung auf die innere Reibung abgesehen wird. 


1) R.Saver, Nichtstationdre Probleme der Gasdynamik. Ein Auszug daraus findet sich in 
E. Jenny, Berechnungen und Modellversuche tiber Druckwellen grosser Amplituden in Auspufflei- 
stungen, Dissertation ETH (Ziirich 1949), S. 74-83. 
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} 
3 
a 
ie 


J Langs den Stérlinien kénnen sich die drei abhangigen Variablen , @, w und 
p die zwei unabhangigen x, ¢ nur nach der Vertraglichkeitsbedingung : 


dw = 4 Faw (In fedt— pepo wt [LF &- Yl a 
| (6) 
d sd cee eed 


andern. Mit den passend gewahlten Bezugsgrossen fiir Druck f,, Lange L, 
Schallgeschwindigkeit a) und a), = 4 (po/p)\*-)2* werden folgende dimensions- 
lose Gréssen eingefiihrt : 


; pa eee eZ foe ne; hae 
| Le eee A ep Se Ae as ae Z= aoe 
ys . &L ts by L 
F,= (in f),L; ——— 2d’ 0, = Po 


Ay; ist die Schallgeschwindigkeit, die durch adiabatische Expansion eines Gases 

von #, a auf f, erreicht wird. Darum geben a,, und damit A, mittelbar die 
- Entropie des Gases an. Die Gleichungen (5) und (6) kénnen sich dann in dimen- 
 sionsloser Form umschreiben: 


adZ i 
i ae AS (Sa) 
iw== Z A, pepsi") a Ay pe-virx W E, aves 5 ci we 


“u—1 


[ W “— 1 —(* 2% 
x [LE Ge) Gopecmi|| 8 + —1 A, P-O*NB* Q,dZ } (a) 


2 x—1)/2% 
=F — 4A.4[ Pt yen) 4+ Af+ Ar+ Aq. 


Die letzten drei Glieder auf der rechten Seite der Gleichung (6a) sind gegentiber 
dem ersten Glied nur Korrekturen, bedingt durch Querschnittsdnderung, 


Aussere Reibung und Warmezufuhr. 
Fiir die Lebenslinie gilt die Energiegleichung (3) als Vertraglichkeitsbedin- 


gung. Durch einige Umformungen wird sie 


Ped ee Up Py one feel 7 (7) 


2 0 (pl-D2% 4 )? EW ‘ 


womit die Entropiednderung entlang einer Lebenslinie bestimmbar ist. 
Werden die drei Korrekturglieder in der Gleichung (6a) zuerst ausser acht 
gelassen, so lasst sich die Beziehung zwischen W und P in einem Koordinaten- 
system (W; P&—?*) durch eine Schar unter + 1/A, geneigter, geradliniger 
Charakteristiken darstellen, wenn der Mafstab fur die P™-1)2*-Achse [2/(% — 1)]- 
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mal grésser als fiir die W-Achse gewahlt wird!). Fiir A, = 1 sind dann die Cha- 
rakteristiken unter + 45° geneigt. Die so gewonnenen Charakteristiken werden 
nachtraglich mit Rticksicht auf die drei Korrekturglieder nach JENNY folgen- 
dermassen korrigiert”): Wenn zum Beispiel eine von links kommende Welle 7 
eine von rechts kommende 3 kreuzt, dann wird die resultierende Welle im 
W-P-Diagramm in erster Annaherung durch den Schnittpunkt 4’ der beiden, 
unter + 1/A, geneigten, durch die Punkte 7 und 3 gehenden Charakteristiken 
dargestellt (vgl. Figur 2). Die entsprechenden Stérlinien 74 und 34 im X-Z- 
Diagramm werden mit Hilfe der fiir das jeweilige A, giiltigen Geraden A = + W 
konstruiert. 45 ist die Lebenslinie. Somit konnen die zur Berechnung der Kor- 
rekturglieder A/, Mv, Aq benétigten Gréssen aus den beiden Diagrammen ent- 


Figur 2 
Konstruktion der Charakteristiken unter Berticksichtigung der Querschnittsanderung des Rohres, 
der dusseren Reibung und der ausseren Warmezufuhr. 


a)Aussere Reibung 5) Querschnittanderung c) Warmezutuhr 
ie ee 
. 


Warmeabfuhr 


SAR 


Richtung der waagrechten Verschiebung der Charakteristiken infolge: a) der dusseren Reibung 
b) der Querschnittsanderung des Rohres und c) der dusseren Warmezufuhr. 


Figur 3 


Abe P.DEHALLER, Uber eine graphische Methode in der Gasdynamik, Sulzer Tech. Rev. 1, 6-24 
*) E. Jenny, a. a. O., S. 83-85. 


x 
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nommen werden. Hierbei sind fiir P, W, A, und F, die Mittelpunkte der be- 
treffenden Linienelemente massgebend. Um die Summe von 


D4 = Af + Ar+ Aq 


verschiebt man nach Figur 3 die Geraden 74’ und 34’ im W-P-Diagramm in 
horizontaler Richtung. Der neue Schnittpunkt 4 stellt dann den korrigierten 
Zustand der resultierenden Welle dar. Weicht 4 stark von 4' ab, so muss die 
Berechnung nochmals verbessert werden. 


NP ee Ce ee AA NE te ee 


2. Stationire Randbedingung 
fiir das Strémen eines Gases durch eine Diise 


In diesem und den folgenden Abschnitten wird nur der Fall reibungsfreier, 
gegen aussen wdrmeisolierter Stromungen behandelt. Die Diffusion der Gase 
in den verschiedenen Schichten wird nicht beriicksichtigt. 

Beim Strémen eines Gases durch eine in einem Rohr befindliche Diise 
- 4ndern sich die Strémungsquerschnitte der Dtisenzone in der x-Richtung 

_ sprunghaft. Die nach x genommenen Differentialquotienten der drei Gleichun- 
gen (1), (2), (3) tiberwiegen deshalb die nach ¢ genommenen. Wenn diese ganz 
vernachlassigt werden, kann die Stromung in der Diisenzone als stationar be- 
handelt werden. Diese stationare Randbedingung ist sehr einfach und wird 
deshalb in den Berechnungen instationarer Stromungen vielfach verwendet. 
Man denke sich die Diise in Figur 4 zunachst durch einen Schieber ver- 
schlossen. Auf seiner linken Seite befinde sich ein Gas mit sehr hohem Druck 5, 
auf seiner rechten Seite eines mit sehr geringem Druck p,, so dass pa<B bo 
(B = kritisches Druckverhaltnis) sei. Wird die Diise plotzlich gedffnet, so 
strémt das Gas mit einer sehr hohen Geschwindigkeit w, durch die Diise von 
links nach rechts. Fiir eine stationare Randbedingung an der Diise ist w, gleich 
der Schallgeschwindigkeit a,. Vor der Diise entsteht dann eine nach links wan- 
dernde Verdiinnungswelle 0-7 und nach der Diise ein nach rechts wandernder 
Verdichtungsstoss s’ a. Im W-P-Diagramm sind die entsprechenden Punkte 
gezeichnet. Von den in Figur 4 angegebenen Zustanden des Gases sind die- 
jenigen fiir das von der Verdiinnungswelle noch nicht gestorte Gebiet des linken 
Rohrteils mit dem Index 0, diejenigen direkt vor dem Diiseneintritt mit dem 
Index 1, diejenigen am Diisenaustritt mit dem Index /, diejenigen fiir den im 
rechten Rohrteil durch die eindringenden Massen ausgelésten Verdichtungs- 
stoss mit dem Index s’ und diejenigen fir das vom Verdichtungsstoss noch 
nicht erfasste Gebiet mit dem Index a gekennzeichnet. Fir die Verdiinnungs- 
welle 0-7 und die Strémung in der Diise 1-# gelten nach EICHELBERG") fol- 


gende Gleichungen: 


1) G. EICHELBERG, Instationdre Stromungsvorgange in Motoren, Forsch. Geb. Ingenieurwesens 


14, H. 2, 41-47 (1948). 
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05 


Patt 
Figur 4 


W-P-Diagramm fiir diestationare Randbedingung in der Diisenzoneeiner Rohrs tr6mung mit Wy = Ay. 


a) Adiabatische Expansion 


oe (2-Die TT, Oh. 
VA ink a ean vic 
b) Zustandsgleichung der Verdiinnungswelle 


Gen eat oder w,—w,—=—~ 
dw 2 1 OS 


Be NS — 4); 
oa (ae a) bzw. rire : 


c) Fortpflanzungsgeschwindigkeit der Welle 


C1 
Ao 


im 


wt aes 
d 
a ay 


d) Energiegleichung fiir die quasistationiire Randbedingung 


w; — wi 


rg = %(T- T) 


Cgc reertlyg, sacaeal «| 


oder 


NE ee 


7a 
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__ e) Kontinuitatsgleichung 
4 F w, 0, = Fy Wy; @; « 


aa Aus diesen fiinf Gleichungen kénnen $,/9, p;/po, Dy/py und ¢,/a) als Funk- 
tionen von F,|F gelést werden. Die Tabelle 1 zeigt den Rechenvorgang fiir den 


0 Ormur ees Ce O05 08 07 SOs 09 10 is 


Figur 5 


Die Verhiltnisse p,/Po, Pz/Po, PzlPy, &1/%q) Ps’ minlPo UDA Ws’ max/% in Abhangigkeit von F,/F fir die 
! stationdare Strémung in der Diisenzone mit wy = ay und * = 1,4. 


_vorliegenden Sonderfall w,; = a, mit x = 1,4, und die Figur 5 stellt die Ergebnisse 


als Kurven dar). 
Fs bleibt noch der Verdichtungsstoss s’-a, also der Punkt P, in Figur 4 


zu bestimmen. Nach der stationaren Randbedingung gelten fiir das Umset- 
zungsgebiet der kinetischen in Druckenergie des sich nach der Diise bildenden 


- Strahls folgende Gleichungen: 
a) Kontinuitatsgleichung 
W; Of Fy = Wy es F; 


b) Impulsgleichung 
$F — &) +p, — by f= we F Og i w# Fy es 


mit p, = Druck auf der Diisenriickwand ; 


1) Y.N.CHEN, Druckwellenspiilung bet Zweitaktmotoren, Mitt. Inst. Thermodynamik u. Ver- 
brennungsmotorenbau, ETH Zirich, Nr. 12, 57-61 und 32-49 (1953). 
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Cc) Energiegleichung 


w? — wp = —*— (a? — a2) . 


-Zunachst wird der extreme Fall behandelt, bei dem eine méglichst hohe Ge- 
-schwindigkeit am Ende der Umsetzung zu erzielen ist. Dann muss f, = 0 sein. 
Somit wird die Impulsgleichung 

3 bs Fy — by F =; Fy Os (Wy — Wy). 

‘Mit 2, = a; 0,/x und py = a% @,/% ergibt sich weiter 


a} of Fy— 43, Oy F 


y, = WF; oy (Wy — Wy) , 


oder 
a2 ={ (1+ x) a,—x wy} wy. 


Nach Einsetzen dieses Ausdruckes in die Energiegleichung folgt 


2 2 (1+ x) x 
as vg ee 2 
ES eee ast eg bet ad oe ree 
oder 
w2 + a?—2a,v,=—0 
und 


Ws = Wy max = 4¢ = Wz 


Der zu dieser Geschwindigkeit gehérende Druck geht aus der Impulsgleichung 


hervor: pF p p a 
b., ESS simin _ Pr , *'t 
Py — Purmcin F oder Po Do F - 


Die héchste erzielbare Geschwindigkeit wy ,,a, kann die kritische Geschwindig- 
keit a, bestenfalls erreichen, kann aberniemals Uberschallgeschwindigkeit sein’). 
Die Werte von y min/Py sind ebenfalls in Tabelle 1 und Figur 5 eingetragen. 
In Figur 4 ist das Wertpaar wy max/4 und (py min/po)~”* durch den Punkt 
(Ws maxi Ps’ min) Aargestellt. Es ist ersichtlich, dass der Druck von py, auf py min 
sinkt, ohne dass die Geschwindigkeit sich erhéht. Der Vorgang ist rein iso- 
thermisch, und die dazu benétigte Warme riihrt ausschliesslich von den Car- 
notschen Stédssen im Strahl her. 

Zu jedem anderen, von Null verschiedenen Druck #, gehdrt ein anderes 
Wertpaar w, und p,, aus welchem sich die Kurve R der Figur 4 bildet. Sie geht 


1) Gemass der Messung von Boris N.CoLE und Brian Mits in The Theory of Sudden Enlarge- 
ments Applied to the Poppet Exhaust-Valve, with Special Reference to Exhaust-Pulse Scavenging, 
Proc. (B) I. Mech. Eng. 1 B, 372 (1952/53), kann p,, = 0 erst bei kleinem F/F angenahert erreicht 
werden. Darum ist die in diesem Abschnitt entwickelte Theorie um so genauer, je kleiner P/F 


als 1 ist. 
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vom Punkt P,iq,, der den Stauzustand des Punktes (w,; #,) darstellt, bis zum 
Punkt Vs mars Ps min) 

Der Zustand des im rechten Rohrteil durch das einstrémende Gas ausge- 
lésten Verdichtungsstosses ergibt sich aus dem Schnittpunkt der Charakteristik 
P,P, mit der Kurve R. Diese Charakteristik weicht um so starker von der 
Geraden ab, je grésser das Druckverhiltnis p,//A, ist. 

Liegt der Punkt P, sehr tief, zum Beispiel bei P,,, so kann es vorkommen, 
dass die Charakteristik P,, P, des Verdichtungsstosses und die Kurve R sich 
nicht schneiden kénnen. Die Stetigkeitsbedingungen von Geschwindigkeit und 
Druck am Ende des Strahls lassen sich dann nicht mehr gleichzeitig erfiillen. 
Sollte nur die Bedingung wy = wy may allein beriicksichtigt werden, so wiirde 
dies zu pyr< po min fahren. Es bestiinde also ein Drucksprung y min — dsr 
zwischen den einstrémenden und den bereits im rechten Rohrteil befindlichen 
Gasmassen. Dieses unwahrscheinliche Ergebnis ist eine unvermeidliche Folge 
der stationaren Behandlung der fiir die Diisenzone geltenden Randbedingung. 


1) Ee7 


Anfangszustand: W4,= 0,98 


C Vakuumn vordemHolben 
Ps PF 


fee 
zy) rs O7131 


Po 
Kein Fekuum vordemtolben 


2-0 0826p, 


Figur 6 


Stationarer Druckverlauf fiir ein in einem Rohr befindliches Gas bei der plétzlichen Exp 


ansion 
ins Vakuum, x = 1,4. 


ae wah ay 
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__ In einem weiteren Beispiel der Figur 6 befinde sich ein Gas mit dem Zu- 
-stand #, und a, im linken Teil eines langen Rohres, dessen rechtes Ende durch 
einen Kolben abgeschlossen sei. Der Kolben bewege sich plétzlich mit der Ge- 
-schwindigkeit 0,9 a) nach rechts fort. Wenn die Randbedingung stationar be- 
-trachtet werden kénnte, so liessen sich die ergebenden Driicke und Geschwin- 
digkeiten aus der Tabelle 1 oder der Figur 5 entnehmen, wie dies in Figur 6 
fiir die drei Falle K/F = 1; 0,6; 0,131 geschehen ist. Das Ergebnis ware: Bei 
F,/F = 1 kénnten die einstromenden Gasmassen den Kolben nicht erreichen, 
“und es entstiinde ein Vakuum dazwischen. Ware die Diisenéffnung vom Ver- 
hiltnis &/F = 0,6, so wiirde das Vakuumgebiet kleiner. Dieses wiirde bei 
_ F;/F = 0,131 vollstandig verschwinden. 


3. Ungenauigkeit der stationdren Randbedingung 
fiir das Strémen eines Gases durch eine Dise 


Sowohl der Drucksprung ymin — Psy nach dem Strahl in der Figur 4 als 
auch das Vakuum vor dem Kolben in der Figur 6 sind aber unvereinbar mit den 
-instationdren Strémungserscheinungen. Die Druckénderungen in den Gasen 
_ werden immer vom Ursprung der Stérung aus durch Druck- oder Verdiin- 
_nungswellen fortgetragen. Der Druck und die Geschwindigkeit entlang dem 
_ Fortpflanzungsweg der Stérung konnen sich nur stetig 4ndern. Ein Druck- 

sprung am Ende des sich nach der Diise bildenden Strahls ist deshalb unwahr- 
_ scheinlich. Das Vakuum stellt sich vor einen sich bewegenden Kolben auch nicht 
bei seiner in der Figur 6 gezeigten Geschwindigkeit ein, sondern bei viel hoherer 
Geschwindigkeit. Beispielsweise betragt diese 5a, fiir ein Rohr mit konstantem 
Querschnitt ohne Diise (/7/F = 1). Die im Abschnitt 2 erlangten Ergebnisse 
weichen nur deshalb von der Tatsache ab, weil sie auf Grund vereinfachter, 
- stationdrer Randbedingungen an der Diise abgeleitet sind!). Solche Rand- 
bedingungen gelten nur fiir jene Unterschallstro6mungen, bei denen das Zeit- 
intervall des zur Behandlung der Stromungsgleichungen benttzten Differenzen- 
verfahrens so gross gewahlt ist, dass ein geniigender Ausgleich der in der 
Diisenzone hin- und herwandernden Wellen erzielt wird. Da bei Uberschall- 
stromungen keine Welle gegen die Stromrichtung zuriickzuwandern imstande 
ist, muss hier die Randbedingung immer instationar behandelt werden. Die im 
rechten Rohrteil der Figur 4 befindlichen Gase werden durch die aus der Diise 
kommenden Druckwellen tiber den Punkt Py hinaus verdichtet, aber héch- 
stens bis zum Punkt P,, dem Schnittpunkt der Verlangerungen der beiden 
Linien P, P, und Py Py (siehe Figur 4). Dieser Punkt gilt fur Rohre mit 


1) Aus stationaren Randbedingungen lasst sich ferner ein gtinstigstes Offnungsverhaltnis 
F ylF ableiten, um in dem an der Duse anschliessenden Rohr einen méglichst starken Verdichtungs- 
stoss zu erzeugen (vgl. Y.-N. CHEN, a. a. O., S.63). Dies gilt natiirlich nach der instationaren Be- 


trachtung der Randbedingungen nicht mehr. 
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konstantem Querschnitt ohne Diise. Er stellt den maximal méglichen Verdich- 
tungsstoss fiir die beiden Anfangszustande 0 und a” dar. 


4. Instationare Randbedingung 
fiir das Str6men eines Gases durch eine Diise 


Eine exakte Behandlung dieser instationdren Randbedingung ist wegen der 
komplizierten, unbekannten Zustande des nach der Diise sich bildenden Strahls 
unméglich. Deshalb muss die Str6mung etwas vereinfacht dargestellt werden. 
Man denke sich den Strahl von einem Diffusor umgeben, der die gleiche Form 
habe wie die freie Grenze des Strahls. Der mit dem Strahl verbundene Carnot- 
Stoss wird dann als Reibungsverlust des Diffusors aufgefasst. Die instationare 
Strémung in der Diise diirfte damit durch eine solche in der gedachten Laval- 
diise reproduziert werden. 

Der Carnot-Stossverlust in dem an der Diise anschliessenden Strahl lasst sich 
angenahert nach der stationaren Energiegleichung: 


id - dw? ap 
lca 


berechnen. Die Lange des Strahls, in der die Umsetzung der kinetischen in 
Druckenergie als vollendet betrachtet werden kann, ist sehr schwer abzu- 
schatzen. Aus der Messung von VIKToRIN}) kann der Offnungswinkel des 
Strahls ungefahr auf 12° geschatzt werden, woraus sich die Strahllange / bestim- 
men lasst. Aus / kann der Ausdruck F, fiir den gesamten Strahl nach dem 
Differenzenverfahren berechnet werden: 


er oF at He 
Te (dt Sly) [Zo bya 


Da der Stossverlust ag als Reibungsverlust aufgefasst wird, kann der Aus- 
druck € w?/(2 d) in den Gleichungen (6) und (6a) durch g az/l ersetzt werden. 
Somit sind dann die Korrekturglieder Ay und 4} fiir den gesamten Strahl nach 
einer der beiden Gleichungen bestimmbar, wobei fiir a, w bzw. P, W, A, die ent- 
sprechenden, im Strahl bestehenden Mittelwerte zu setzen sind. Auf diese Weise 
lasst sich die instationdre Strémung im Strahl berechnen. 

In der obigen Betrachtung wird die gesamte Strahllange / als ein einziges 
Langenintervall des Differenzenverfahrens genommen. Es ist unméglich, das 
Verfahren durch Aufteilung des Strahls in mehrere Langenintervalle zu ver- 


feinern, da sich der Zustand des. Gases innerhalb des Strahls nicht berechnen 
lasst. 


1) K. Vixtorin, Untersuchungen turbulenter:M ischvorgdnge, Forsch. Geb. Ingenieurwesens 12, 
H. 1, 16-30 (1941). 
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E Wie genau dieses Differenzenverfahren ist, wird im folgenden Beispiel (Figur 
7) gezeigt. Die Anordnung der Strémung sei genau gleich der in Figur 4 ge- 
e zeigten. Die Rohrleitung sei in beiden Richtungen unendlich lang. Die urspriing- 
: lichen Zustinde o und a der auf beiden Seiten der Diise befindlichen Gase sind 
_ so gewahlt, dass die Charakteristik des Verdichtungsstosses ASimin Zerade durch 
_ den Punkt Sinin(Ws’ max; Ps’ min) Zeht (siehe Figur 7). In gentigender Zeit nach 
dem Offnen der anfangs geschlossenen Diise stellt sich ein Beharrungszustand 
; in der Diisenzone ein. Der Druck #, auf die Diisenriickwand ist in diesem Fall 
_ gleich Null. Die Zustande der Gase vor dem Diiseneintritt 7, am Diisenaustritt/ 


Figur 7 
W-P-Diagramm fiir die Strémung in der Dusenzone bei p, = 0 und py = Po min: 


und am Ende des Strahls s’,,;, kénnen dann nach Tabelle 1 oder Figur 5 abge- 
~ lesen werden. Diese Zustande sind aber nichts anderes als diejenigen der letzten, 
‘unendlich schwachen Wanderwelle, die sich nach dem vollstandigen Wellen- 
ausgleich in der Diisenzone bildet und durch diese fortpflanzt. Eine solche 
Welle lisst sich natiirlich auch nach der instationaren Auffassung durch die 
Gleichung (6) oder (6a) berechnen. Hierzu zieht man zundchst durch den 
Punkt f eine unter —1/A, geneigte Gerade } B mit A, = Mittelwert von 4, 
zwischen den Punkten f und s‘,;,. Dann wird zu ihr die Parallele CD gezogen, 


mit einem waagrechten Abstand von 
aS Afar: 


Ist das Verfahren ganz genau, so sollen die beiden Punkte C und D mit jenem 


Sinen Zusammenfallen, Wenn nicht, so stellt das Verhaltnis 
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die Ungenauigkeit des Verfahrens dar. Tabelle 2 zeigt ein Beispiel fiir &/F = 
0,910; 0,6 und 0,131. 

Die Umsetzung der Energien im Strahl geht im vorliegenden Fall wegen 
Pr = 0 isothermisch vor sich. Der Carnotsche Stossverlust wird somit nach der. 
Formel 


Ss 


dwt | fap 
lige ae 


f 


Why 


= RT,1n oa f__ (R= Gaskonstante) 


s’ min 


berechnet. Die Korrekturglieder Af und Ar sind dann 


a a af dt a et | 
Af= | — w (In fle dt = rahe meee ag OE EE Rg = (8) 
oder 
A a, a eee 
AT= aoe Lk EL (wegen a = W = a; = ws) 
und 
& tr nee? w 
Ane at a [2 («—1) >] ae 
€ 1 |w| w l 
= (aah) tare 1! OD gl 
Ue ea} w 1 
Ar = (ap 8) 3-75 [1—- #-) 2] (9) 
oder 
. aR fee) 72 a xe 1 sap ts ae 
Ar ae Pile "Bled See ACHCR (wegen a= w=a,= w,). 


Wie aus der Tabelle 2 zu ersehen ist, ist die neue Berechnungsmethode recht 
genau. Erst wenn J+/F viel kleiner als 1 ist, wird die Genauigkeit schlechter. 

Somit haben die Gleichungen (8) und (9) ihre Brauchbarkeit zur Berechnung 
instationarer Strémungen in Diisenzonen gezeigt. 

Hier wird das Beispiel der Figur 4 mit sehr tiefem Anfangsdruck #, (dem 
Punkt P,, entsprechend) nochmals behandelt, allerdings vom Gesichtspunkt 
der instationdren Randbedingung aus. Die Zustande der Gase vor und nach 
der Diise seien mit den gleichen Buchstaben und Indizes gekennzeichnet wie in 
der Figur 4. Nach dem plotzlichen Freigeben der Diise wird der an ihrer Aus- 
trittsmiindung sich augenblicklich einstellende Zustand durch den Schnitt- 
punkt m der beiden durch 0 und a” gehenden Charakteristiken dargestellt 
(siehe Figur 8). Da es sich hier um den engsten Querschnitt der gedachten 
Lavaldiise handelt, sind df = 0 und ag = 0. Darum werden die beiden Charak- 
teristiken Om und a”m so konstruiert wie fiir instationdre Strémungen in 
einem Rohr mit konstantem Querschnitt. Die Druckwelle m-—a" wandert in 
den gedachten Diffusor hinein. Wenn sie an den geraden Rohrteil gelangt, ist 


A 
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Tabelle 2 (% = 1,4) 


Bemerkungen 


(ParmentPo)” 7" | 0,829 | 0,799 | 0,673 Py miniPy aus Tabelle 1 


(Pelpy) "YI? * 0,84 0,86 0,90 b;|Pp aus Tabelle 1 
ripe 0,7045 | 0,7402 | 0,8099 | aus Tabelle 1 
az|ay 0,84 0,86 0,90 = (Ppl Po) "~ D?* 
In (byl Ps min) 0,096 0,513 2,032 
ap|R Ty 0,0676 | 0,38 1,645 = (Ty[To) In (PsP 5 min) 
Ar 0,0173 | 0,0945 | 0,391 = (ap/R Ty) (2 — x)/(2 % as/a) 
Af 0,0396 | 0,215 0,692 = (a/a,) (1 — F/F)/(1 + B/F) 
DA 0,0569 | 0,31 1,083 =BD 
A, 1,012 1,075 1,34 = (ap|a) (Pol Per min) 2” 
As, mittol 1,006 1,04 14% = (1+ 4,)/2 
SC gee: 0,057 0,317 1,33 vgl. Figur 7 
y ~0 2% 19%, ee oe DS eB) 


sie auf die Intensitat n-a” vermindert. Der Punkt x ergibt sich aus der Ver- 
schiebung der Geraden 0m um 2'A yn = Af + Ar. 

Die Wanderwelle im Strahl ist, strenggenommen, nicht linearer, sondern 
raumlicher Natur. Darum stimmt der Punkt m nur angendhert. Seine richtige 
Lage liegt vielmehr etwas héher auf der durch 0 gehenden Charakteristik. Die 

_ Abweichung ist um so grésser, je kleiner F;/F gegeniiber 1 ist. 

Die Verdiinnungswelle 0-m kann nur mit ihrem Unterschallteil 0-f von 
der Diisenaustrittsmiindung aus stromaufwarts in den linken Rohrteil wandern. 
Sie wird am Diiseneintritt zum Teil reflektiert. Die reflektierte Welle wandert 
nach dem Diisenaustritt zuriick und wird hier wieder zum Teil reflektiert usw. 

Nach dem Ausgleich solcher hin- und herwandernden Wellen stellen sich die 


SS ESE eee E 
‘gfe 5 
o {Oe 
f 
(. ach P Smin 
f s S yee. 


Zon a’ 
0(w,=0) 
Figur 8 
W-P-Diagramm fiir die instationare Randbedingung in der Diisenzone einer Rohrstromung mit 
We = Ay. 
if f 


ZAMP VI/19 


Mag, 7 
290 Y1an-NIAN CHEN ZAMP_ 


stationaren Zustinde 7-/’ zwischen Diisenein- und -austritt ein, welche aus : 
der Tabelle 1 bzw. Figur 5 abgelesen werden kénnen. In der gleichen Zeit wan-_ 
dert der Uberschallteil fm als Verdiinnungswelle durch den Strahl strom-— 
abwarts ab, wahrend das Wellenelement mit Schallgeschwindigkeit am Diisen-_ 
austritt stehenbleibt. Die Verdiinnungswellenelemente mit staérkerer Uber- 
schallgeschwindigkeit verlassen sehr bald den Strahl und treten in den Rohrteil — 
mit konstantem Querschnitt ein. Ihnen folgen dann die Verdiinnungswellen-_ 
elemente mit schwacherer Uberschallgeschwindigkeit. Nach gentigender Zeit — 
bleiben im Strahl praktisch nur Verdiinnungswellenelemente mit Schallge- 
schwindigkeit zuriick, da sie Fortpflanzungsgeschwindigkeit gleich Null besit- 
zen. Inzwischen hat sich die Schallgeschwindigkeit am Diisenaustritt von f in 
f’ verwandelt. Man zeichne nach Tabelle 1 bzw. Figur 5 den der Schallgeschwin- 
digkeit /’ zugeordneten Punkt si,in (Wy maz) Ps’ min): Dann bedeutet die Gerade 
f'Simin Cie Zustande der Gase im Strahl. Man suche den mittleren Wert fiir A, 
zwischen den beiden Zustanden /’ und s‘,,,, und ziehe durch den Punkt /’ die 
Charakteristik /’B mit der Neigung —1/A,. Der waagrechte Abstand s7,,;,-B 
bedeutet dann die Summe der Korrekturglieder von Af + Ar fiir den Strahl. 
Die endgiiltige, unendlich schwache Druckwelle, die am Diisenaustritt den 
Zustand /’ erlangt hat, wandert mit dem Zustand sj,,,,, aus dem Strahl in den 
geraden Rohrteil hinein. Sie trifft zuerst die aus der Diise eingedrungene Gas- 
masse und schreitet dann tiber die urspriinglich im Rohr befindliche, durch die 
vorangehenden Verdichtungsstésse erfasste Gasmasse hinweg. Diese beiden 
Gasmassen besitzen verschiedene Werte von Entropien und damit verschiedene 
A,. Man muss daraus ein mittleres A, ermitteln, um den Verlauf der Charak- 
teristik Sjnz,—S Zu bestimmen. Der Schnittpunkt s dieser Charakteristik mit 
derjenigen des Verdichtungsstosses as stellt dann den endgiiltigen Zustand 
der Front des nach rechts wandernden Verdichtungsstosses dar. Die beiden 
Kurven unter der Skizze der Rohrleitung in Figur 8 zeigen den f- und w-Ver- 
lauf der Str6mung nach dem Erreichen des Beharrungszustandes in der Diisen- 
zone. c, ist die Fortpflanzungsgeschwindigkeit des Verdichtungsstosses, |a,—w,| 
diejenige der ihm folgenden Verdiinnungswelle und w, die Eigengeschwindigkeit 
des Gasteilchens. Der Punkt s mit dem Pfeil w, stellt die Trennschicht der er- 
wahnten beiden Gasmassen dar. cy ist die Fortpflanzungsgeschwindigkeit der in 
den linken Rohrteil schreitenden Verdiinnungswelle 0-7. 

Nach der gleichen Methode kann die Berechnung des in Figur 6 dargestellten 
Beispiels revidiert werden (vgl. Figur 9). Die parallel zur P-Achse, durch 
w/a) = 0,9 gehende Gerade K bedeutet die Randbedingung am sich bewe- 
genden Kolben. Ihr Schnittpunkt m mit der Charakteristik 0-/-m stellt den 
Zustand des Gases am Diisenaustritt fiir den Augenblick dar, in dem der 
Kolben sich plétzlich von der Diise weg mit der Geschwindigkeit 0,9 ay bewegt. 
O-m ist die dabei am Diisenaustritt entstehende Verdiinnungswelle. Im Bild a 
geht der Unterschallteil 0-/ stromaufwarts in die Diise zuriick und wird in 


wap VD 
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dieser hin- und herreflektiert, bis die Beharrungszustande 1 und /’ am Diisen- 
ein- bzw. -austritt bestehen. 

Man suche den zum Punkt /’ gehérenden Punkt s;,,,, berechne sein A, 
und ziehe durch sj,,;, die unter —1/A, geneigte Charakteristik s},;,-s. Die 
Gerade ms stellt die Anderung des Zustandes des Gases an dem sich bewegenden 

_ Kolben dar. s bedeutet seinen endgiiltigen Zustand. 
Fiir das Offnungsverhaltnis F;/F = 0,6 sind nach den Tabellen 1 und 2: 


4 _ 0,96; —Pemin _ 0,208: A,(s’nin) = 1,075, 
ao Po 


woraus sich ergeben: 


(gz) = 0,82, Pm = 0,25 Bo 


oft =0137 


06 


ies 
Smin 


Figur 9 


W-P-Diagramm fiir das Stromen eines Gases ins Vakuum, wobei die im Randgebiet stehende 
Strémung instationar behandelt wird. 
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und : 
pg \(*%-D/2* _ ba i ee = 
i = 0,799 — 0,04 75-5 = 0,791, p,=0,147 py 
0-1 ist die endgiiltige, in den linken Rohrteil schreitende Verdiinnungswelle. — 
Fir den Fall /;/F = 1 (Bild 6) schrumpft die Gerade ms zu einem einzigen _ 
Punkt m zusammen und die Punkte 1, /’ und / zu einem einzigen f. 
Fir F,/F = 0,131 sind die beiden Punkte s und sj,,,,, gleich. Somit wird 


Shean Ls (Pemen\ = 0,673 


und 
Ps = Pirmin == 0,0626 Po- 


Fir i/F < 0,131 kommt die Kurve R mit der Geraden K zum Schnitt. 
Der Schnittpunkt s stellt in diesem Fall den endgiiltigen Zustand des Gases am 
Kolben dar. 


5. Erzeugung eines méglichst starken Verdichtungsstosses 
durch Verkleinerung des Rohrquerschnittes 


In Figur 10 grenzt eine Membran M das im linken Rohrteil auf hohen Druck 
komprimierte Gas von dem im rechten Rohrteil befindlichen unter niedrigem 
Druck ab. Die Buchstaben 0 und a kennzeichnen die Zustande der beiden Gase. 
Wird die Membran plétzlich zum Bruch gebracht, so wird das Gas an der 
Grenzstelle M augenblicklich auf den Zustand 3 expandiert bzw. komprimiert. 
Der Unterschallteil 02 von der Verdiinnungswelle 03 (durch gestrichelte Linien 


x 


Figur 10 
Pl6étzliche Expansion eines in einem Rohr befindlichen Gases auf sehr niedrigen Druck. 


Tee 


> 
yy 
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dargestellt) wandert nach links und ihr Uberschallteil 23 nach rechts, wahrend 
der Verdichtungsstoss 3a (durch ausgezogene Linie k dargestellt) sich natiirlich 


- nur nach rechts fortpflanzt. Das Wellenelement 2 mit Schallgeschwindigkeit 


_ bleibt aber immer an der Stelle stehen, wo die Membran war. Die Str6mung 
- wird durch diese feststehende Schallgrenze in Uber- und Unterschallgebiete 


_ aufgeteilt. 


In Figur 11 wird der Querschnitt des rechten Teils des Rohrs von F aut Fy 
reduziert. Die Membran befindet sich etwas rechts vom konischen Rohriiber- 


Figur 11 


_ Ein Verdichtungsstoss wandert ab Stelle M nach rechts und eine Verdiinnungswelle nach links. 


Diese schreitet mit Unterschallgeschwindigkeit durch den konischen Rohrtibergang. 


gang yz. Es soll untersucht werden, wie diese Querschnittsverengerung die 
Strémung und damit die Starke des nach rechts wandernden Verdichtungs- 
stosses 3a beeinflusst. 

Nach dem Bruch der Membran geht der Unterschallteil 02 der Verdiinnungs- 
welle 03 von der Membran aus durch den konischen Ubergang zy in das linke 
Rohr hinein, wahrend der Uberschallteil 23 nach rechts wandert. Durch Hin- 
und Herreflexion der Welle 02 zwischen den beiden Enden z und y des Uber- 
ganges werden schliesslich in diesem die Ausgleichszustande 1 — f (w, = ay) 
hergestellt, welche aus der Tabelle 1 bzw. der Figur 5 fiir das betreffende Off- 
nungsverhiltnis F;/F abgelesen werden konnen. Der Schnittpunkt s der beiden 
Charakteristiken fs und as stellt dann den Endzustand des nach rechts schrei- 
tenden Verdichtungsstosses dar. Dieser erhalt also gegentiber der Figur 10 einen 
Zuwachs entsprechend dem Charakteristikenstiick s3. Die Schallgrenze der 
Strémung kann aber ihre Lage nicht mehr wie in Figur 10 bei der Membran M 
behaupten, sondern sie muss stromaufwarts nach dem engsten Querschnitt z des 
konischen Uberganges yz wandern. 

Wird andererseits die Membran in den Rohrteil urspriinglichen Quer- 
schnittes F verlegt (vgl. Figur 12), so wandert der Verdichtungsstoss 3a nach 
dem Bruch der Membran in den Ubergang yz hinein. Nachdem der Stoss diesen 
durchwandert hat, wird sein Zustand von 3 in 4 verwandelt, indem die Refle- 
xionswellen « und f an der Stelle z bzw. y verursacht werden. Wenn die von y 
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ausgehende Reflexionswelle f die Stelle z erreicht hat, wird ihr Zustand in 4’ 
verwandelt. Die’ beiden Reflexionswellen « und 8 wandern weiter in den geraden 
Rohrteil hinein und werden von der Welle 1, welche die Reflexion der an die . 
Stelle y gelangten Verdiinnungswellenfront 3 darstellt, iiberholt. Die eigentliche © 


p-Verlauf : 


Figur 12 
Eine hinter einem Verdichtungsstoss mit Uberschallgeschwindigkeit wandernde Verdiinnungs- 
welle schreitet durch den konischen Rohriibergang. Fiir den Druckverlauf Z, ist die Wegachse nicht 
maBstablich gezeichnet. Die Linie t —-— -— stellt die Trennschicht zwischen den beiden, ur- 
sprunglich durch die Membran M getrennten Gasmassen dar. 


Front 3 geht als Welle y in den geraden Rohrteil mit dem Zustand 4” hinein, 
welcher durch die an der Stelle y reflektierte Welle m des ndchsten Verdiin- 
nungswellenteils 3’ bestimmt ist. 

Die Punkte 4, 4’ und 4” werden durch horizontale, parallele Verschiebung 
der entsprechenden Charakteristiken um die Korrektur A f ermittelt. Der Index 


ny Sa 
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ais 


hee 


_von Af in Figur 12 kennzeichnet die Strecke, auf welche sich Af bezieht. Af 
“ muss zundchst angenommen, alsdann nachgepriift werden, ob es dem Zustand 
_ des Mittelpunktes der bezogenen Strecke entspricht. 
P Ein weiteres Verdiinnungswellenelement 3” (als Welle ¢) kreuzt die aus dem 
_ nichsten Element 3” stammende Welle m an der Stelle z zu einem resultierenden 
_ Zustand entsprechend dem Punkt 4”, welcher just auf der Geraden A = + Ww 
liegt. Das heisst, die Geschwindigkeit w,' ist gleich der Schallgeschwindigkeit. 
Die der Welle » folgenden Wellen — eine davon ist zum Beispiel die Welle f — 
_ bringen dann die Geschwindigkeit der Welle ¢ unter die Schallgeschwindigkeit, 
- so dass die Welle ¢ einen Wendepunkt an der Stelle z hat und nach dem Punkt 
_ 4” stromaufwarts biegt. Die Figur 12 zeigt, wie die aus dem Verdiinnungswellen- 
~ element 73 stammende Welle p jene ¢ beim Punkt 74 kreuzt, der oberhalb der 
Geraden A = + W liegt. 
Das Unterschallgebiet breitet sich mit der riickschlagenden Front der Ver- 
' diinnungswelle ¢ ab Stelle z stromaufwarts aus, wodurch ein Wellenausgleich 
in der Ubergangszone yz erméglicht wird. Nach geniigender Zeit werden die 
- gleichen Beharrungszustande «7» an der Stelle y und «f» an der Stelle z er- 
reicht wie die in der Figur 11. Die Schallgrenze der Stromung muss wieder ihre 
_ Lage bei der Membran verlassen und nach dem engsten Querschnitt des koni- 
schen Ubergangs abwandern. 
Man zieht durch die gewonnenen Punkte 4’, 4", 4" und f die Parallelen zur 
Charakteristik 03 und erhalt die Schnittpunkte 41 4 4Ul und s auf der Ge- 
 raden 4a. (Die Punkte 7, f und s sind der Ubersichtlichkeit halber nicht in das 
W-P-Diagramm der Figur 12 eingezeichnet. Sie haben die gleichen Lagen wie 
in Figur 11.) Diese Punkte stellen die Zustande der dem Verdichtungsstoss k 
folgenden Verdiinnungswellen l,m,n, p dar. Gleichzeitig wandert eine end- 
giiltige Verdiinnungswelle o? ab Membran stromaufwarts. 

Es sind drei Zeitquerschnitte durch den X-Z-Plan der Figur 12 gemacht, 
naimlich Z = Z,, Z = Z, und Z= Z;, um den mit der Zeit veranderlichen 
Druckverlauf zu veranschaulichen. Fiir Z; ist praktisch der Beharrungszustand 
7 f im konischen Ubergang yz erreicht. 

Unter sonst gleichen Bedingungen ist hier der infolge der Verkleinerung des 
Rohrquerschnittes erzielte Verdichtungszuwachs der Stossfront um die Grésse 
ds starker als im vorangehenden Beispiel. Die stromabwarts der Membran ge- 
machte Verkleinerung des Querschnittes kann also einen viel heftigeren Ver- 
dichtungsstoss erzeugen. Ist das rechte, aussere Ende des Rohres geschlossen, 
so wird der Verdichtungsstoss dort gestaut, wodurch Druck und Temperatur 
noch weiter ansteigen. Auf diese Weise ist ein Mittel geschaffen, um Gase auf 
sehr hohe Temperaturen zu bringen und ihre dabei auftretenden Phanomene Zu 


studieren?). 


1) In E.L. RESLER, S.C. Lin, A. KANTROWITZ, The Production of High Temperature Gases in 
Shock Tubes, J. appl. Phys. 23, H. 12, 1390-1399 (1952), wird dieses Thema behandelt. 
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Zusammenfassung 


Die durch Querschnittsanderung gebildete Ubergangszone einer eindimen-— 
sionalen, instationaren Gasstr6mung wird im allgemeinen als stationar behandelt, 
weil eine exakte Lésung der betreffenden Strémungsgleichungen nicht existiert. 
Dies ist nur fiir jene Strémungen mit Unterschallgeschwindigkeit zulassig, bei 
denen der Zeitabstand zweier aufeinanderfolgenden Beobachtungen so gross ge- 
wahlt ist, dass ein gentigender Ausgleich der in der Ubergangszone hin- und her- 
pendelnden Wanderwellen sichergestellt wird. Anders ist es aber bei Str6mungen 
mit Uberschallgeschwindigkeit, wo ein solcher Wellenausgleich nicht md6glich ist. 

In der vorliegenden Arbeit wird versucht, die Randbedingung in der Uber- 
gangszone instationar zu behandeln. ie 

Es zeigt sich, dass die Schallschranke, die am Austritt der einfachen Diise in 
der stationdaren Strémung uniiberwindbar ist, in der instationdaren Str6émung 
nicht mehr besteht. Sowohl am Diisenaustritt als auch in dem nach der Diise sich 
bildenden Strahl kénnen Uberschallgeschwindigkeiten auftreten. 


Summary 


The transition zone of a one-dimensional non-steady gas flow formed by 
crossectional variation is generally treated as steady, because no exact solution 
of the gas equations concerned exists. The results are valid only for subsonic gas 
flow, in which the time interval between two successive observations must be so 
chosen as to ensure sufficient equalization of the oscillating waves in the transition 
zone. Under supersonic conditions, however, it is impossible to attain such a wave 
equalization. 

In this paper it is endeavoured to treat the boundary condition as non-steady 
in the transition zone. Furthermore, the sound barrier which cannot be overcome 
at the exit of a simple nozzle with steady flow, no longer exists with non-steady 
flow. 


(Eingegangen: 30. Juni 1954.) 


Rechnungen zur Beurteilung der Flattersicherheit 
von Flugzeugen’) 


Von EUGEN GRUSCHWITZz, Altenrhein?) 


Die Schwingungen, deren ein Flugzeug als elastisches Gebilde fahig ist, 
k6énnen je nach seinen mechanischen und aerodynamischen Eigenschaften ober- 
halb einer gewissen Fluggeschwindigkeit oder auch in einem begrenzten Ge- 
schwindigkeitsbereich durch Aufnahme von Energie aus dem Luftstrom ange- 
facht werden. Wenn diese Energieaufnahme durch Strémungsvorgange erfolgt, 


*) Dieser Aufsatz ist auf Wunsch der Schriftleitung als Einftthrung zu der nachfolgenden 
Arbeit von U. Hocusrrasser verfasst worden. 
2) Flug- und Fahrzeugwerke AG. 
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e die naherungsweise Potentialstrémungen darstellen, spricht man von Flattern. 
_ Eine Energieaufnahme durch Strémungsvorgange, die nicht durch Potential- 
_ strémungen angendhert werden kénnen, findet zum Beispiel beim Schiitteln 
- (englisch: «buffeting») statt, das bei hohen Fluggeschwindigkeiten auftreten 
kann. Es wird durch lokale Uberschallgebiete mit instabilen Verdichtungsstés- 
_ sen verursacht: 
Da beim Flattern die bei weitem gefahrlichste Anfachung erfolgt und da 
die Energieaufnahme durch Potentialstromungsvorgange in gewissen Grenzen 
der Berechnung zuganglich ist, sucht man die Flattersicherheit eines Flugzeug- 
_musters auf rechnerischem Wege zu beurteilen. Man geht dabei nach der 
Methode der kleinen Schwingungen vor und sucht insbesondere diejenige Flug- 
geschwindigkeit festzustellen, nach deren Uberschreitung eine Anfachung még- 

- lich ist, die sogenannte kritische Geschwindigkeit. Die Ungenauigkeiten in den 
Luftkraften, in der Feststellung des elastischen Verhaltens und der Massen- 

 verteilung des Flugzeuges sowie die Idealisierung des Flugzeugs oder einzelner 
seiner Teile durch irgendwelche der Rechnung zuganglichen Systeme lassen die 
Resultate von Flatterrechnungen in mancher Hinsicht fragwiirdig erscheinen. 
Die Erfahrung hat jedoch gezeigt, dass solche Rechnungen trotzdem sehr gute 
Dienste fiir die Verhiitung des Flatterns leisten kénnen. Ihre quantitativen 
Resultate miissen natiirlich mit einem gewissen Sicherheitsvorbehalt bewertet 
werden. Aber sie liefern auch in qualitativer Weise wertvolle Einblicke in das 
Flatterverhalten und insbesondere lasst sich durch Veranderung einzelner Kenn- 
gréssen die Auswirkung konstruktiver Massnahmen auf das Flatterverhalten 
verfolgen. 

Zu den unsichersten der in die Flatterrechnungen eingehenden Kenngréssen 
gehéren die instationdéren Luftkrafte, und zwar in um so héherem Masse, je 
mehr (im Unterschallbereich) die Kompressibilitat der Luft in Betracht zu 

_ ziehen ist. Schon in stationarer Strémung ist die Berechnung der Auftriebsver- 
teilung iiber einem Fliigel endlicher Spannweite umstandlich. Die raumliche 
instationare Strémung um einen unverformten Fliigel kann mit grossem Re- 
chenaufwand behandelt werden, jedoch diejenige um einen sich elastisch ver- 
formenden Fliigel macht vorlaufig noch zu grosse rechnerische Schwierigkeiten. 
Man hilft sich deshalb mit der sogenannten Streifenmethode. Dabei nimmt man 
an, dass jeweils an einem beliebig schmalen, in Flugrichtung liegenden Streifen 
des Fliigels dieselben Luftkrafte angreifen wie an einem entsprechenden Streifen 
eines Fliigels in ebener Strémung. Die instationaren Luftkrafte an einem Fliigel 
mit angelenktem Ruder in ebener Strémung liegen tabuliert vor. Die Energie- 
iibertragung von dem Luftstrom auf den Fliigel oder umgekehrt kommt dadurch 
zustande, dass die durch die Verformungen des schwingenden Fliigels hervor- 
gerufenen Luftkrafte phasenverschoben zu den Auslenkungen sind. Dement- 
sprechend werden die instationdren Luftkrafte durch komplexe Gréssen dar- 
gestellt. Die Machsche Zahl der Anstrémung geht in die instationéren Luft- 
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krafte als Parameter ein. Man kann die Rechnung jeweils nur mit Luftkraften _ 


fiir eine bestimmte, von vornherein zu wahlende Machsche Zahl durchfihren. 


® 
Z| 


Genau genommen miisste man sie mit Luftkraften fiir andere Machsche Zahlen _ 
wiederholen, die den Bereich tiberdecken, in dem man Flatterfreiheit nachzu- | 


weisen hat. Beispielsrechnungen haben jedoch gezeigt, dass die kritischen Ge- 
schwindigkeiten, die man mit Luftkraften fiir verschiedene Machsche Zahlen 
erhalt, im allgemeinen nicht weit voneinander entfernt liegen. Man rechnet 
deshalb mit Luftkraften fiir eine Machsche Zahl, die in der Nahe der oberen 
Grenze des als flatterfrei nachzuweisenden Geschwindigkeitsbereiches liegt. Da 
nach der Prandtlschen Regel die Strémung um einen mit der Machschen Zahl M 
angestrémten Tragfliigel der inkompressiblen Str6mung um einen Tragfliigel 
entspricht, dessen Abmessungen in Spannweitenrichtung mit dem Faktor 
V1 — M? verkiirzt sind, miissen um so gréssere Abweichungen von den nach der 
Streifenmethode gerechneten Luftkraften erwartet werden, je naher die 
Machsche Zahl der Anstré6mung an 1 liegt. 

In der sehr umfangreichen Literatur tiber das Flattern werden viele Még- 
lichkeiten fiir den Ersatz des Flugzeuges oder eines Fliigels durch ein fiir die 
Rechnung geeignetes System behandelt. Da ein Flugzeug ein Kontinuum mit 
unendlich vielen Freiheitsgraden ist, kann man die Bedingung fiir die kritische 
Geschwindigkeit durch eine Integralgleichung darstellen, die das Gleich- 
gewicht der elastischen Krafte mit den Luft-,; Massen- und Tragheitskraften in 
jedem Punkt ausdriickt. Fiir die Rechnung ist eine Algebraisierung des Pro- 
blems, das heisst die Reduktion auf eine endliche, nicht zu grosse Anzahl von 
Freiheitsgraden, erforderlich. Die beiden Verfahren, die in der nachfolgenden 
Arbeit von HocHstTrAssER behandelt werden, gehen von der Integralgleichung 
aus, und bei beiden wird das Flugzeug durch ein «System vorgegebener Form- 
anderungsfahigkeit» ersetzt. Es wird dem Flugzeug dabei nur die Verformung 
in einer endlichen Anzahl von vorgegebenen Elementarformen und linearen 
Kombinationen daraus erlaubt. Die Anzahl » der vorgegebenen Elementar- 
formen ist die Anzahl der Freiheitsgrade des Systems, die Faktoren, mit denen 
die Elementarformen multipliziert werden, sind seine dynamischen Koordina- 
ten. An Stelle der Integralgleichung hat man dann ein System von m linearen 
homogenen Gleichungen fiir die dynamischen Koordinaten. 

Die beiden erwahnten Verfahren unterscheiden sich durch die Wahl der 
vorgegebenen Elementarformen. Bei dem einen wird eine Auswahl aus den 
Eigenfunktionen des Flugzeuges in ruhender Luft (bezichungsweise im Vakuum 
~— der Unterschied ist im allgemeinen unerheblich) als Elementarformen genom- 
men. Dieses Verfahren diirfte das jetzt am haufigsten angewandte sein. Man 
wahlt also aus den unendlich vielen Eigenschwingungen des Flugzeuges jeweils 
diejenigen aus, von denen man aus anschaulichen Betrachtungen annehmen 
kann, dass durch ihre Kombination die Schwingungsform eines Flatterfalles 
angendhert werden kann. HocustRassEr hat fiir die Ableitung der hierfiir giil- 
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tigen Flattergleichungen einen Weg gewahlt, der bei KUssNER?) vorgezeichnet 
ist. Die Gleichungen werden dabei direkt aus dem Gleichgewicht der Krafte her- 
d geleitet und nicht, wie das sonst in der Regel geschehen ist, aus dem Variations- 
-prinzip der Mechanik oder den Lagrangeschen Gleichungen. Auf diese Weise 
werden die Schwierigkeiten vermieden, die das Vorhandensein komplexer Krafte 
und komplexer Amplituden bei Anwendung des Variationsprinzipes oder der La- 
grangeschen Gleichungen verursacht. Fiir das Flattern kommen nur die unteren 
_Eigenschwingungen in Betracht, da hohe Frequenzen keine nennenswerten Luft- 
_krafte liefern. Dazu kommt, dass bei hohen Frequenzen mehrere Knotenlinien 
und somit Vorzeichenwechsel der Luftkrafte vorhanden sind, so dass ein teil- 
_weiser Ausgleich stattfindet. Fiir den Fall, dass die Luftkrafte klein sind im 
Vergleich zu den Tragheitskraften, ist ohne weiteres klar, dass die Eigenfunk- 
tionen des Flugzeugs in ruhender Luft die besten Elementarformen liefern. 
~Doch auch in Fallen, in denen diese Voraussetzung nicht erfiillt ist, kann damit 
eine brauchbare Annaherung erwartet werden. Eine Plausibilitatsbetrachtung 
_hieritber findet sich bei BoRKMANN?). 
‘Das zweite Verfahren ist die von RAUSCHER®) in die Flatterrechnung ein- 
_ gefiihrte Methode der Stationsfunktionen. Dabei wird die Integralgleichung in 
-n Punkten erfiillt, und zwischen diesen Punkten werden die Verformungen 
durch eine lineare Kombination von » Funktionen interpoliert, die so gewahlt 
sind, dass sie im wesentlichen die Randbedingungen erfiillen. 
Das elastische Verhalten des Flugzeuges wird bei der ersten Methode durch 
die Eigenfunktionen, die Eigenfrequenzen und die Massenverteilung wieder- 
gegeben: In der Integralgleichung wird der Deformationstensor (beziehungs- 
weise die Einflussfunktionen) nach den Eigenfunktionen entwickelt. In der 
Rauscherschen Methode bleiben die Einflussfunktionen in der Integralgleichung 
stehen. Bei der notwendigen Beschrinkung auf wenige Elementarformen 
-k6nnte, falls die Luftkrafte von der gleichen Gréssenordnung sind wie die 
Tragheitskrafte, die Darstellung der Einflussfunktionen durch wenige Eigen- 
funktionen ungeniigend erscheinen. Die Einflussfunktionen werden jedoch nur 
selten genauer bekannt sein, als sie sich durch wenige Eigenfunktionen darstel- 
len lassen. Wenn am fertigen Flugzeug gemessene Eigenfrequenzen und Eigen- 
funktionen («Standschwingungsformen ») vorliegen, ist naturgemass die erste 
Methode die gegebene. Wenn genaue Einflussfunktionen vorliegen und die 
Luftkrafte von gleicher Gréssenordnung wie die Tragheitskrafte sind, kann 
moglicherweise die Methode der Stationsfunktionen bei gleichem Rechenauf- 
wand etwas wahrscheinlichere Resultate liefern. 


1) H.G. Ktssner, Fliigelflattern und verwandte Erscheinungen. Deutsche Ausgabe der Fiat 


Review of German Science, Bd. 11, S. 151. 
2) K.BorKMANN, Einfihrung der Baugrossen eines Flugzeuges m dessen Flatterberechnung, 
Jahrbuch 1941 der Deutschen Luftfahrtforschung, S. I, 540 (insbes. S. 544 und 545). 
3) M. RAUSCHER, Stationfunctions and Air Density Variations in Flutter Analysis, J. Aeron. Sci. 


16,345 (1949). 
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Der rechnerische Aufwand fiir Flatterrechnungen ist betrachtlich. Rech- 
nungen mit mehr als drei Freiheitsgraden sind praktisch nur mit programm- 
gesteuerten Rechenmaschinen zu bewAltigen. In der Arbeit von HOCHSTRASSER 
wird deshalb das Flatterproblem vom Standpunkt des im Umgang mit solchen 
Maschinen Erfahrenen aus behandelt. 


Summary 


It is given a short introduction to some problems of flutter-calculations, 
especially concerning the reprensentation of the aeroplane by systems of given 
deformations. : 


(Eingegangen: 14. Mai 1955.) 


Flatterrechnung mit Hilfe von programmgesteuerten 
Rechenmaschinen 


Von Urs HocustrasseEr, Ziirich?) 


§ 1. Einleitung 


Die Flatterrechnung befasst sich mit den selbsterregten Schwingungen von 
Flugzeugen, die infolge einer Stérung des stationaren Flugzustandes auftreten 
kénnen. Je nach der Geschwindigkeit der Maschine sind diese Schwingungen 
durch Energieabgabe an den Luftstrom gedampft oder infolge Energieauf- 
nahme angefacht. Die Geschwindigkeit, bei der weder Energie abgegeben noch 
aufgenommen wird, bezeichnet man als kritische Geschwindigkeit. Diese kann 
ohne irgendwelche Warnungszeichen iiberschritten werden, so dass dann eine 
Stérung des Flugzustandes unter Umstinden zu einer so stark angefachten 
Schwingung fiihren kann, dass das Flugzeug explosionsartig zerstért wird, 
Deshalb ist man bestrebt, diese Geschwindigkeit schon vor der Flugerprobung 
auf theoretischem Wege mit der Flatterrechnung zu ermitteln. 

Die ausgedehnte Literatur iiber die Flatterrechnung zeugt von der Wichtig- 
keit und Kompliziertheit dieses aerodynamischen Problems. Die Schwierig- 
keiten sind sowohl physikalischer als auch mathematischer Natur. Auf phy- 
sikalischem Gebiet ist wegen der komplexen Struktur der Flugzeuge die Be- 
schaffung der fiir die Rechnung notwendigen Daten, wie Massenverteilungen, 
Einflussfunktionen, Tragheitsmomente und Luftkraftkoeffizienten, meist sehr 
schwierig. Auch ist die fiir die mathematische Formulierung notwendige Ideali- 
sierung des physikalischen Problems oft in ihrer Tragweite nicht leicht iiber- 


1) Flug- und Fahrzeugwerke AG., Altenrhein, St. Gallen. 
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sehbar. Auf mathematischem Gebiet sodann ist das Problem fast immer, 
selbst bei Verwendung sehr stark vereinfachender Annahmen, nicht exakt lés- 
bar, so dass Approximationsmethoden verwendet werden miissen. Diese Ver- 
fahren sind meistens ziemlich grob, da der Arbeitsaufwand fiir eine genauere 
_ Rechnung sehr rasch zunimmt. Es hat auch gar keinen Sinn, die Giite der 
Approximation allzuweit zu treiben, da der oft recht grosse Unsicherheitsfaktor 
in den Ausgangsdaten grosse Fehlerschranken fiir die Resultate zur Folge hat. 
Bis zur Einfiihrung von programmgesteuerten Rechenmaschinen ist jedoch 
diese Begrenzung der Genauigkeit von der numerischen Seite aus kaum er- 
reicht worden, da der Aufwand fiir die Rechnung mit Handrechenmaschinen 
sehr gross ist. Ohne einen gréssern Stab von geiibten Rechnern liess sich bis vor 
kurzem die Flatterrechnung kaum ausfiihren. Heute ist an die Stelle des Stabes 
von Rechnern die programmgesteuerte Rechenmaschine getreten, die viel 
rascher und meist auch zuverlassiger arbeitet als der Mensch. Dies erlaubt 
einerseits die Approximation des mathematischen Problems weiterzutreiben, 
andererseits die Formulierung des Vorganges besser der Wirklichkeit anzu- 
_passen. Dementsprechend werden wir in der vorliegenden Arbeit zuerst eine 
Ubersicht geben iiber Ansdtze und Methoden, die gestatten, den Flattervorgang 
ohne allzu grosse Idealisierungen beziiglich der Form und des elastischen Ver- 
haltens der schwingenden Korper zu erfassen. Weiter werden wir numerische 
Lésungsverfahren angeben, die sich speziell fiir die Verwendung mit programm- 
gesteuerten Rechenmaschinen eignen. Zum Schlusse soll dann noch kurz be- 
schrieben werden, wie die uns zur Verfiigung stehenden Rechenmaschinen zur 
Lésung des Flatterproblems eingesetzt wurden. 

Die numerischen Methoden, die bei der Rechnung mit der Maschine ver- 
wendet werden, sind nicht unbedingt die gleichen wie die fiir die Handrechnung 
verwendeten. Da die Maschine, verglichen mit dem Menschen, ein ausserordent- 
lich dummer Rechner ist, fiir den alle Rechenvorginge bis ins kleinste Detail 
festgelegt werden miissen, wird man fiir sie méglichst einfache Verfahren be- 
nutzen. Je nach der Rechengeschwindigkeit der Maschine spielt es keine so 
grosse Rolle, ob das gewahlte Verfahren am raschesten zur Lésung fihrt. Da 
neben der Rechengeschwindigkeit auch die Programmierungsméglichkeiten 
variieren kénnen, hangt die Wahl des numerischen Verfahrens auch noch vom 
Typus der zur Verfiigung stehenden Maschine ab. Die in dieser Arbeit beschrie- 
benen Verfahren sind deshalb nicht unbedingt in jedem Falle die besten oder 
raschesten, sie haben sich jedoch fiir die von uns beniitzten Rechenmaschinen 


bewahrt. 


§ 2. Mathematische Formulierung des Flattervorganges 


Das Flugzeug stellt ein mechanisches System mit unendlich vielen Frei- 
heitsgraden dar, das im stationdren Zustand eine gleichférmige translatorische 
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Bewegung ausfiihrt. Eine Stérung des stationaren Zustandes bewirkt das Auf- 
treten von zusatzlichen Luftkraften. Fiihrt man die Tragheitskrafte als Schein- 
krafte ein, so kann man die dynamische Aufgabe, die durch die Stérung ver- 
ursachte Bewegung zu erfassen, mit Hilfe eines elastostatischen Ansatzes als 
statisches Problem formulierent). Als Koordinatensystem verwendet man ein 
mit dem Flugzeug fest verbundenes rechtwinkliges System mit den Achsen 
%4, X, Xz. Die am Flugzeug im Zeitpunkte ¢ zusdtzlich angreifenden Luft- und 
Tragheitskrafte kann man formal durch eine Kraftdichte K,(P, t) (j = 1, 2, 3) 
darstellen, die im Punkte P wirkt. Bezeichnet man mit w,(P, t) die relative 
Verschiebung des Punktes P des Flugzeuges in der x,;-Richtung infolge der 
Stérung des stationdren Flugzustandes, mit 7;,(P, P’) den symmetrischen 
Deformationstensor, so bestehen nach dem Hookeschen Gesetz die folgenden 
Beziehungen: 


3 ; 
w,(P, t) = i b ‘hi Dy TilP, P) KulP’,t) dey dag dee j7=1,2,3. (1) 
(F) i‘: 


Die Integration ist itber das ganze Flugzeug zu erstrecken. 

Da eine Stérung des stationdren Flugzustandes nur dann ungefahrlich ist, 
falls die relativen Verschiebungen in Abhangigkeit von der Zeit gegen Null 
konvergieren oder héchstens eine harmonische Schwingung mit kleiner Ampli-— 
tude ausfiihren, wird fiir die zeitliche Abhangigkeit von w;(P, t) der folgende 


Ansatz gemacht: oe 
w,(P, t) = w,(P) e'”". (2) 


Die Zeitabhangigkeit der Kraftdichte K,(P, t) kann in derselben Weise ab- 
gespaltet werden, weil die Kraftdichte eine homogene Funktion ersten Grades 
beztiglich der Verschiebungen und ihrer zeitlichen Ableitungen ist. Da K,(P, 2) 
die Resultierende von Luft- und Tragheitskraften ist, kann man sie nach Ab- 
spaltung der Zeit in der folgenden Weise darstellen: 


K,(P, t) =»? {M(P) w,(P) + 2, L[w,(P), w}} eh 4s 1s ee 


wobei M(P) die Massenverteilung im Punkte P, w = w/v die «reduzierte Fre- 
quenz» und v die relative Geschwindigkeit des Luftstromes sowie J eine cha- 
rakteristische Lange des Korpers ist. Der erste Term auf der rechten Seite von 
Formel (3) riihrt von den Tragheitskraften her, wahrend der zweite die Luft- 
krafte darstellt. Da die Luftkrafte eine Phasenverschiebung gegeniiber den 
Schwingungen des angestrémten Korpers besitzen, ist L;{w,(P), w] eine kom- 
plexwertige Vektorfunktion. Die Abhangigkeit von L,{w,(P), w] von der redu- 
_ zierten Frequenz @ ist so kompliziert, dass sie nur in tabellarischer Form 


1) Vergleiche zum Beispiel: KUssner, Fliigelflattern, S. 151 der deutschen Ausgabe der Fiat 
Review of German Science, Bd. II (Chemie, WEINSTEIN). 
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E gegeben werden kann. Setzt man (2) und (3) in (1) ein, so erhalt man ein System 
_ von 3 Integralgleichungen, das die Zeit nicht mehr enthalt: 


; 2 w)(P) = Sif f [ra. Ps) (MiP) w,(P’) + = L,[w.(P’); «| (4) 
bi Fie 


: SOU ESR CORT La a ee A 


_ wobei A = 1/v? gesetzt wurde. 4 und w sind Parameter in diesem System, wobei 
_ nur fiir bestimmte Wertepaare das System (4) eine Lésung besitzt. Das System 
- (4) stellt also ein zweiwertiges Eigenwertproblem dar. In der Flatterrechnung 
- interessiert man sich speziell fiir reelle Wertepaare A und «, fiir die (4) eine 
Lésung besitzt, da diese reellen Frequenzen v entsprechen und deshalb die 

_ Grenze zwischen gedampften und angefachten Schwingungen bestimmen. Das 
_ System (4) kann nur unter Verwendung von vereinfachenden Annahmen geldst 
werden. So werden praktisch nur Teilkérper des Flugzeuges (Fligel, Rumpf, 

_ Seiten-und Hohenleitwerke usf.) untersucht, da gefahrliche Kopplungenzwischen 
solchen Teilkérpern nur bei Nachbarschaft der Frequenzen auftreten kénnen. 

- Zudem macht man noch besondere Annahmen iiber das elastische Verhalten 
_ der Teilkérper, die eine Vereinfachung der Integralgleichung erlauben. Wir 
miissen jedoch verzichten, im Rahmen dieser Arbeit darauf naher einzugehen. 
Auch nach der Einfiihrung all dieser Annahmen in (4) kann dieses System 
gewohnlich nur approximativ gelést werden. Wir wollen hier zwei Approxima- 
tionsmethoden darstellen, die beide die Integralgleichungen in ein homogenes 
lineares Gleichungssystem iiberfithren. Beide beruhen darauf, dass man den 
Deformationstensor naherungsweise durch einen Ausdruck der folgenden Form 


ersetzt: 
Cee es De X4(P) Y,,(P’) (nm ganze Zahl), jy lL 2 one me) 
r=1 
wobei vorausgesetzt wird, dass die X*(P) ein System von » linear unabhan- 
gigen Vektorfunktionen bilden und diese sowie die Y,,(P) quadratisch inte- 


grable Vektorfunktionen iiber dem betrachteten Integrationsbereich sind. Setzt 
man (5) in (4) ein, so erhalt man eine Entwicklung des Verschiebungsvektors 


in den Vektorfunktionen X3(P): 
n 3 
aw, (P) = 3 X5P) ff [Sf Yel?) 
r=1 (F) k=1 
x {mr w,(P") + 7 Lilw.(P’), co}lax dx, dx (6) 


= Sa xi | 
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wobei fiir die g, wegen der vorausgesetzten linearen Unabhangigkeit der X5(P) 
und der Eigenschaften von L;, die Beziehungen bestehen: 


1 : ; 
aaa ff Yl) 
(ey 2 H1 
ee , 1 if Vip, , ! , 
x M(P YS an Apes Ste late (Py, ol ax, Gx,'d%," (S192 a) 
ie 


rt 


oder, wenn man die g, mit gleichem Index zusammenfasst : 


hdr = DY) Ure (®) Ip (7) 


mit 


(F) (8) 
, 1 App , , , 
x {M( PY AUP) LAP), w}} dx! dx, dx. 


Man erhalt also ein homogenes lineares Gleichungssystem fiir die Entwick- 
lungskoeffizienten g,, das nur nichttriviale Lésungen besitzt, falls die Koeffi- 
zientendeterminante, die sogenannte « Flatterdeterminante», verschwindet: 


Lire, 


D,,(A, oo) = || re (co) =A O-r || =0, wobei Ory = (9) 


Dh assed | 


Die u4,,(@) sind wegen der Glieder mit L, in (8) komplexwertige Funktionen der 
reduzierten Frequenz w, deren Abhangigkeit von w nicht in einfacher Weise 
darstellbar ist. Zur Bestimmung der kritischen Geschwindigkeit miissen die 
reellen Wertepaare A, w gefunden werden, fiir die die Determinante D,(A, @) 
verschwindet. Man kann dieses Problem auf ein lineares Eigenwertproblem zu- 
rickfihren, indem man die reduzierte Frequenz w als Parameter auffasst und 
die n Eigenwerte 4;, die im allgemeinen komplex sein werden, in Abhangigkeit 
von w bestimmt. Die reellen Nullstellen 4* kénnen dann aus dem Verlauf von 
4,(m) aufgefunden werden. Mit Hilfe der zugehérigen reduzierten Frequenz 


ih eee ae 1 


kann die kritische Geschwindigkeit v* berechnet werden, fiir die die durch die 

Stérung verursachten Schwingungen weder gedampft noch angefacht sind. 
Die Genauigkeit der so bestimmten kritischen Geschwindigkeiten hangt 

nattirlich wesentlich von der Giite der Approximation (5) des Deformations- 
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tensors ab. Wie wir sehen werden, erhalt man die beste Approximation des 
4 Deformationstensors, falls man die Eigenfunktionen /;(P) und die Eigenwerte 
A, des einfacheren Schwingungsproblems ohne Luftkrafte beniitzt. Da dieses 
Problem so formuliert werden kann, dass die Kerne der Integralgleichungen 
_symmetrisch und positiv definit sind: 


U 


7, VM (P) f(P) 


=) | [Vm@ nue, Py VRP VME) KP) ax, axa, | 


k=1 (F) 


_kann man einen Satz von ScumipT?) verwenden. Auf das Integralgleichungs- 
system (10) iibertragen besagt dieser, dass die Kerne /M(P) T;,(P, P’) /M(P’) 


durch Summen der Form 3’ «3(P) f,;,(P’) fiir gegebenes 7 am besten durch 
r=1 
DV M(P) 4, fi(P) fil P') VP’) 
r=1 


angenahert werden. Das heisst also, dass die beste Approximation des Defor- 
mationstensors in diesem Falle durch 


TylP, P= VEAP MP) j,2=1,2,3, (1) 


gegeben ist. Die Verwendung dieser Approximation bedingt die Kenntnis der 

Eigenvektorfunktionen und Eigenwerte des Problems ohne Luftkrafte. In man- 

chen Fallen kénnen diese experimentell gemessen werden. Eine andere Méglich- 

keit ist ihre Berechnung auf Grund der Integralgleichungen (10). Dies ist ein 

viel leichteres Problem als das entsprechende Problem mit Luftkraften, da in 
diesem Falle die Kerne symmetrisch und positiv definit sind. 

Ein anderes Verfahren zur Approximation der Lésungen des Integralglei- 
chungssystems (4) liefert die von RAUSCHER?) vorgeschlagene Methode der 
Stationsfunktionen. Bei dieser Methode verlangt man, dass der Ansatz (5) fiir 
n Punkte P = P,, P,, ..., PB, exakt giiltig ist fiir jede Wahl von k und P’, Als 
System unabhangiger Vektorfunktionen Xj(P) wahlt man in diesem Falle 
3 m Vektorfunktionen g%(P), deren Komponenten Stationsfunktionen genannt 
werden. Diese erfiillen die folgenden Bedingungen: 

1 falls r=Rk4+(f7-1)n, vr=1,2,3, SSE MEY ie 

8i(Pi) = ; (12) 
0 falls r+k+(7—1)n, k=1,2,...,%, p= 1y2;3; 

1) E, Scumipr, Uber die Auflésung der nichtlinearen I ntegralgleichungen, Math. Ann. 65 ( 1908). 

2) M. RauscueER, Stationsfunctions and Air Density Variations in Flutter Analysis, J. Aeron. Sci. 


16, no. 6, 345-354 (1949). 
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aus denen die lineare Unabhingigkeit der g%(P) sofort ersichtlich ist. Die 
Y,,(P’) sind in diesem Falle Einflussfunktionen, da wegen (12) 


a 


das heisst Y,,,(;-1)n,x(P) die Einflussfunktionen fiir die Verschiebungskompo- 
nente w; im Punkte P, sind. Aus den Beziehungen (6) folgt dann, dass 9,4. (3-1) n 
in diesem Falle im wesentlichen gleich der Komponente w; des Verschiebungs- 
vektors im Punkte P. ist. Daraus folgt auch, dass bei dieser Methode die Inte- 
gralgleichungen des Systems (4) an den Stationen exakt erfiillt sind. Deshalb 
gehort die Methode der Stationsfunktionen zu den Kollokationsmethoden. 

Die Giite dieser Approximation hangt sehr von der geeigneten Wahl der 
Stationen P. und der Stationsfunktionen ab. Es ist offensichtlich, dass die Sta- 
tionen in solche Punkte gelegt werden sollten, die bei einer Deformation eine 
fiir das elastische Verhalten des Kérpers charakteristische Verschiebung erlei- 
den. Da die Stationsfunktionen zusammen mit den Einflussfunktionen den 
Deformationstensor méglichst gut approximieren sollten, ist es empfehlens- 
wert, die Stationsfunktionen aus den Einflussfunktionen durch Superposition 
unter Beriicksichtigung der Bedingungen (12) zu konstruieren. Die Methode 
der Stationsfunktionen eignet sich besonders zur Ermittlung der Eigenvektor- 
funktionen und Eigenwerte fiir das Problem ohne Luftkrafte, falls die Einfluss- 
funktionen bekannt sind. Fiir die eigentliche Flatterrechnung gibt jedoch die 
Methode mit den Eigenvektorfunktionen nach den vorausgegangenen Bemer- 
kungen im allgemeinen eine bessere Approximation. 

Die Materialdémpfung wurde bei der Formulierung des Flattervorgangs der 
Einfachheit halber weggelassen, da die durch sie verursachten Phasenverschie- 
bungen bei den Verschiebungskomponenten keine prinzipiellen Anderungen in 
der Flatterrechnung zur Folge haben. Des weiteren ist zu beriicksichtigen, dass 
das mit dem Flugzeug fest verbundene Koordinatensystem noch Bewegungen 
mit héchstens sechs Freiheitsgraden ausfiithren kann. Diese zusatzlichen Frei- 
heitsgrade kénnen jedoch mit Hilfe der Gleichgewichtsbedingungen (Summe 
aller Krafte und Summe aller Momente miissen verschwinden) sofort eliminiert 
werden, da die Gleichgewichtsbedingungen A nicht enthalten. An dieser Stelle 
kann gerade auch noch erwahnt werden, dass die Integralgleichungen eigentlich 
noch von einem dritten Parameter abhangen, da die Luftkrafte eine Funktion 
der Luftdichte sind, die mit der Flughohe variiert. Die Beriicksichtigung dieses 
Parameters geschieht auf analoge Weise wie bei der reduzierten Frequenz. 
Man bestimmt fiir eine Reihe diskreter Werte der Luftdichte die kritischen 
Geschwindigkeiten und ermittelt dann allfallige weitere interessante Werte 
durch graphische Interpolation. 


Ti(P,. P) SY age ee ee eee (13) 


Fah} 
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§ 3. Numerische Methoden 
fur die Ermittlung der kritischen Geschwindigkeiten 


Fiir die numerische Lésung des Problems spielt es keine Rolle, welches der 
_ beiden Approximationsverfahren verwendet wird; in beiden Fallen erhalt man 
_ nach § 2 das folgende lineare Gleichungssystem: 


é | 1G 3 95 | 
PB [4is(@) — 2 6;5] 9; = 0 Nala (14) 


Die Koeffizienten yu; ;(@) sind komplexwertig und haben die folgende Form: 
Uii(@) = M;;+ K;,(@) , 


wobei die M,; von den Tragheitskraften herriithren und deshalb reell sind; die 
K,,; hingegen sind durch die Luftkrafte gegeben und komplex. Beide Ausdriicke 
bestehen im allgemeinen aus dreifachen Integralen tiber das gesamte Volumen 
des betrachteten Systems. Durch Einfiihrung vereinfachender Annahmen 
beziiglich des elastischen Verhaltens und der Luftkrafte sind meist zwei Inte- 
grationen direkt ausfiihrbar, so dass M;; und K,; in der folgenden Form dar- 
gestellt werden kénnen: 


b 
ae 1 E,,(x) G(x) dx, E,,(x) = reelle Funktion, 

: wobei G(x) =reelle Funktion, (15) 
Ki; = [ Eulx) (x, w) dx, H (x, w) = komplexe Funktion. 


Zur Bestimmung der kritischen Geschwindigkeit hat man fiir eine Reihe von 
Ree Frequenzen die Nullstellen der Koeffizientendeterminante 
| “:;(@) — 2 6;;| von (14) zu untersuchen. Dabei miissen die Parameterwerte 
so a, werden, dass sie in dem Gebiet, fiir das eine kleine kritische Ge- 
schwindigkeit méglich ist, geniigend dicht liegen, so dass der Nullstellenverlauf 
in Funktion von w verfolgbar ist. Aus dem Verlauf jeder Nullstelle konnen dann 
fiir die Berechnung der kritischen Geschwindigkeit diejenigen Punkte ermittelt 
werden, wo / reell wird. Falls der Rang der Koeffizientendeterminanten nicht 
sehr hoch ist, wird die Untersuchung der Nullstellen am besten mit Hilfe der 
charakteristischen Polynome durchgefiihrt. Ein einfaches und rasches direktes 
Verfahren zur Berechnung des charakteristischen Polynoms aus der Deter- 


| 
‘ 4 
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ree 5 
minante ist von VOETTER?) angegeben worden. Bei diesem Algorithmus ersetzt 


man schrittweise eine Zeile der Determinante durch eine solche Kombination : 
simtlicher Zeilen, dass zum Schluss die Elemente unterhalb der Hauptdiago-— 
nalen alle verschwinden und das letzte Element der Hauptdiagonalen gerade das ~ 
charakteristische Polynom darstellt. Bei einem allgemeinen Schritt geht man 
in der folgenden Weise vor: 

Es sei die Koeffizientenmatrix in den vorhergehenden (rv — 1) Schritten auf_ 


die folgende Form gebracht worden: 


br,» Merja ee% bry — Aone . nn 
Mrs. 
Aceed Wee Hone Han —A 
ae 2b ene (Pde Sarg 2h in ere 
Oo 
O20 O» brate Ue msectsOpine te Oriaais be 


Dann erhalt man im nachsten Schritt die neue Matrix D, durch Multiplikation 
von D,_, mit der Matrix B,: 


Dp, oer 1a 1 Psa > 
wobei 
Os 0 
OFS Orel aeO 0 
Bite 0 
0 0) 1 
I Ay,» As, n—-L A 7 ap, n 


Die Elemente dieser Matrix sind durch die folgenden Formeln bestimmt: 


fir? Sm—r A, ,=—by,,, wobel by ¢= lng 


n-?r m 


me 255 = 
frase nim ty Apg sages aks yp Ay aa prep as, enh On etal ete 
’ j=1 


‘eet 
(m1, 2,3... 7) 5 


! Ag) ae ee Uber die numerische Berechnung der Eigenwerte von Sakulargleichungen, ZAMP 
: : 


Een 


Mace 
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we 


Wie man leicht verifiziert, wird durch fortgesetzte Multiplikation mit solchen 
_Matrizen die Koeffizientendeterminante von (14) in (m — 1) Schritten auf die 
_vorher beschriebene Form gebracht. Im letzten Schritt vereinfachen sich die 
' Formeln, da b,_1,1; gemeinsamer Faktor von allen A,,, ist, so dass man ihn 
_ wegkiirzen kann, womit 

5 n—2 

: Anan = Aga Moa, — pa A nas, 141 95, 0-5 


7 


eae 


wird, wahrend die itibrigen A,_,,;, nach den Formeln (16) berechnet werden. Die 
Fortpflanzung der Rundungsfehler geschieht ahnlich wie beim Gaufschen 
Eliminationsverfahren zur Auflésung von linearen Gleichungssystemen. Des- 
_halb kann bei umfangreichen Determinanten die Genauigkeit durch Run- 
dungsfehler stark beeintrachtigt werden. 
Falls man den Verlauf der Nullstellen in Abhangigkeit von w explizit 
-kennen will, muss man die Nullstellen der charakteristischen Polynome be- 
rechnen. Sobald der Grad der Polynome grésser als zwei ist, sind speziell auch 
im Hinblick auf die Verwendung von programmgesteuerten Rechenmaschinen 
die iterativen Verfahren zur Auffindung der Nullstellen die gegebenen Metho- 
_ den. Fiir langsame Rechenmaschinen, die das Wurzelziehen als eingebaute Ope- 
ration besitzen, eignet sich das Laguerresche Iterationsverfahren!) am besten, 
da es kubische Konvergenz besitzt. Bei dieser Methode lautet die Korrektur- 
formel zur Verbesserung einer Versuchslésung A, fiir eine Nullstelle eines Poly- 
noms P,(A) n-ten Grades: 


F,,(A1) Pia) = dog 5 (A) 
ae (Tm n n ’ 
AA, = P,,(A1) apae ae = wobei S 
1+ / — 4) [» —1-—n a q | Pi (a) = ae Pays 


Gegeniiber dem Newtonschen Verfahren ist der Aufwand zur Berechnung der 
Korrekturen grosser, dafiir hat man eine viel bessere Konvergenz. Oft kennt 
man die ungefahre Lage der Nullstellen des Polynoms P(A, w), da man schon 
diejenigen fiir eine benachbarte reduzierte Frequenz w+ Aw bestimmt hat, 
und diese in Funktion von nicht stark variieren. Liegen die Nullstellen nicht 
sehr nahe beisammen, so kann man die einfachere Formel fiir die Korrektur 
verwenden: 


— Prltr) 
Pi (hy) 
Ah, = ——q—Piti,) BA) iS) 
7 BAL 


Auch dieses Verfahren konvergiert kubisch, sobald man so nahe an der gesuch- 


1) H.J. MAuty, Zur iterativen Auflésung algebraischer Gleichungen, ZAMP 4, 261 (1954). 


(17) 


310 Urs HocHSTRASSER ZAMP 


ten Nullstelle ist, dass die benachbarten nicht mehr stéren. Bei beiden Verfah- 
ren kann die Rechnung stark abgekiirzt werden, wenn man jeweils das Polynom 
durch die gefundene Nullstelle dividiert, so dass ein neues Polynom niedrigeren 
Grades entsteht, das nur noch die unbekannten Nullstellen enthalt. Dabei ist 
wegen der Rundungsfehler darauf zu achten, dass die dem Betrage nach kleinste 
Nullstelle wenn méglich zuerst bestimmt wird. Hat man fiir eine gréssere Zahl 
von Parameterwerten w die Nullstellen der Flatterdeterminanten berechnet, 
so kann man durch Auftragung des Imaginar- und Realteiles jeder Nullstelle in 
Abhangigkeit von w die rellen Werte auf graphischem Wege ermitteln. 

Dieses Vorgehen zur Bestimmung der reellen Nullstellen erfordert schon bei 
niedriggradigen charakteristischen Polynomen einen erheblichen Rechenauf- 
wand. Dabei ist ein Teil der berechneten Nullstellen physikalisch nicht inter- 
essant, da sie fiir keinen Wert von w reell werden, ausser fiir w~ = oo. Falls man 
den Verlauf der Nullstellen in der komplexen Zahlenebene nicht ben6tigt, kann 
man mit einer anderen Methode den Rechenaufwand fiir die Ermittlung der 
reellen Nullstellen erheblich vermindern. Bei dieser stellt man das Polynom 
P(A, @) mit komplexen Koeffizienten durch Aufspalten der Koeffizienten in 
Real- und Imaginarteil als Summe zweier Polynome von n-tem sowie (x — 1)- 
tem Grade mit reellen Koeffizienten dar: 


n n n-1 
P, (A, w) = 3” aj(w) 4 =D) aj(w) ai +i >| 4; (@) #7 =G, (A, o) +7 Gye 
j=l j=1 j=1 


wobei a;(w) = aj(w) + 7 aj(w). Mit Hilfe dieser neuen Polynome kann man als ~ 
erstes feststellen, in der Nahe welcher der ausgewdhlten reduzierten Frequenzen 
uberhaupt eine reelle Nullstelle von P,(A, w) zu finden ist. Dazu benutzt man 
die mit diesen Polynomen gebildete Sturmsche Kette!) von Polynomen, die 
rekursiv durch die folgende Formel gegeben sind: 


[oe ae ae 


Das letzte Glied in der Kette ist eine Konstante, die auch den Wert Null anneh- 
men kann. In diesem Spezialfall besitzt P,(A, w) eine reelle Nullstelle. Wenn 
die Zahl der Vorzeichenwechsel der Werte von aufeinanderfolgenden Polyno- 
men in. der Kette fir {= — ‘co gleich 7, fiir 1 =co gleich R& ist, so legen 
n/2 + (j — k)/2 Nullstellen in der obern komplexen Halbebene. Mit Hilfe dieser 
Abzahlung der Nullstellen kann man diejenigen w-Werte ermitteln, fiir die 
P,(A, ) mindestens eine Nullstelle in der Nahe eines Uberganges von der 
untern in die obere Halbebene besitzt. Da man gewohnlich fiir die gewahlten 
diskreten Parameterwerte nicht direkt eine reelle Nullstelle bekommt, muss 


Pmi(A) ~ Pm—a1(A) [A i Dinat | i Pm-2(A) , Wwobei 


1) F.Wiiers, Methoden der praktischen Analysis, 2. Aufl. (W. de Gruyter, Berlin 1950), S. 288. 
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man diese durch Interpolation aus den Daten fiir die mittels der Sturmschen 
-Kette ermittelten benachbarten reduzierten Frequenzen bestimmen. Dazu 
_berechnet man mit einem der beschriebenen Iterationsverfahren die Nullstellen 
A, von G,_,(A, w) fiir diese Frequenzen. Da G,_,(A, w) reelle Koeffizienten 
“pesitzt, sind diese reell oder konjugiert komplex. Setzt man die A, in die ent- 
_sprechenden Polynome G,,(A, «) ein, so wird G,,(A;, w) im allgemeinen von Null 
_verschieden sein, ausser wenn A; mit der reellen N ullstelle von P,,(A, w) zusam- 
-menfallt. Durch Auftragung von G,(A,, w) und A, in Abhangigkeit von m kann 
“man graphisch die Nullstellen von G,,(A,, w) und dann auch die reellen Null- 
stellen 4, und die zugehérigen reduzierten Frequenzen bestimmen. Da bei dieser 
~Methode nur die Nullstellen von reellen Polynomen berechnet werden missen, 
ist der Aufwand fiir die Bestimmung der reellen Nullstellen wesentlich kleiner. 

Hat man jedoch in den Ansatzen die Materialdampfung nicht beriicksichtigt, 
so kann es geschehen, dass die Informationen, die dieses Verfahren liefert, nicht 
-ausreichen, so dass man auf die direkte Bestimmung der Nullstellen von 

P,(A, w) zuriickgreifen muss. Dies ist dann der Fall, wenn eine Nullstelle schon 

bei relativ grossen reduzierten Frequenzen in die obere Halbebene wandert und 
- fiir kleinere w-Werte wieder in die untere Halbebene zuriickkehrt. Physikalisch 
heisst das, dass zwischen zwei Geschwindigkeiten ein Anfachungsbereich liegt, 
- der jedoch je nach der Starke der Anfachung bei Beriicksichtigung der Material- 
_ dampfung verschwinden kann. Nimmt man namlich die Materialdampfung wie 

iiblich in Form einer Phasenverschiebung der Deformation um den Winkel » 
_an (der Einfachheit halber soll jeder Freiheitsgrad dieselbe Phasenverschiebung 

besitzen), so besteht zwischen dem Eigenwert A fiir die Schwingung ohne Mate- 

rialdimpfung und demjenigen /* fiir die Schwingung mit Materialdimpfung 

die Beziehung: 


A= A" (L-+itgg). (19) 
_ das heisst, falls A’ = A’ + 74", A* = 4*’ + i 4*" sowie tgy und A” kleine Gréssen 
gegentiber A’ sind, wird: Ae cw As aH ed — I tee. (20) 


Daraus sieht man, dass auch fiir positive Imaginarteile von A der Imaginarteil 
von 4* negativ bleiben kann, so dass also bei Beriicksichtigung der tatsachlichen 
Verhiltnisse mit Materialdimpfung keine angefachte Schwingung in diesem 
Geschwindigkeitsbereich méglich ist. Der Entscheid itiber die Verwendung des 
einen oder des anderen Verfahrens zur Bestimmung der reellen Nullstellen muss 
demnach von Fall zu Fall getroffen werden. 


§ 4. Der Einsatz von programmgégesteuverten Rechenmaschinen 
fiir die numerischen Rechnungen 


Wie in der Einleitung schon erwahnt wurde, ist der Aufbau der numerischen 
Auswertung fiir die Flatterrechnung und die verwendeten Rechenmethoden 
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von der Art der zur Verfiigung stehenden Maschinen abhangig, da die konstruk- 
tiven Daten den Zeitaufwand fiir die gesamte Rechnung wesentlich beeinflus- 
sen. Wir hatten die Méglichkeit, die programmgesteuerte Relaisrechenmaschine 


Z 4 des Institutes fiir angewandte Mathematik der ETH. sowie die Lochkarten- — : 


maschinen der International Business Machines Co. fiir unsere Rechnungen zu 
bentitzen. Die Rechenmaschine der ETH. kann nicht ganz siebenstellige Zahlen 
mit gleitendem Komma verarbeiten und besitzt 64 Speicherzellen. Das Rechen- 
programm darf beliebig viele arithmetische und auch bestimmte logische Ope- 
rationen enthalten. Das Programm wird auf Filmstreifen gelocht, die dann die 
Maschine steuern. Da das Rechenwerk mit Relais arbeitet, ist die Rechen- 
geschwindigkeit nicht sehr gross. Im Mittel werden 1000 Rechenoperationen 
pro Stunde ausgefiihrt. 

Es bestehen verschiedene Typen von IBM-Lochkartenmaschinen, wobei der 
' einzelne Typus nur wenige Funktionen ausiiben kann. Die Lochkarten spielen 
dabei die Rolle des Speichers, in dem die Ausgangsdaten und Zwischenresultate 
sowie in beschranktem Ausmasse auch einfache Steuerbefehle fiir einzelne 
Maschinen auf dem Weg von Maschine zu Maschine in Form von eingestanzten 
Léchern aufbewahrt werden kénnen. Fiir unsere Rechnung standen uns die 
folgenden Maschinen zur Verfiigung: 

Key Punch. Mit dieser Maschine kann man die Ausgangsdaten und Kenn- 
ziffern in die Karten einstanzen. Die Kennziffern dienen zur raschen Identifi- 
kation der auf den Karten eingelochten Zahlen. 


i 
é 


4 
: 
: 

| 
: 
’ 


: 


Sorter. Erlaubt die Sortierung der Lochkarten nach einer bestimmten _ 


Reihenfolge der Kennziffern oder das Herausgreifen einer oder mehrerer 
Karten. 

Reproducer. Besorgt die Reproduktion von Zahlen und Kennziffern von 
einer Lochkarte auf eine andere. 

Collator. Erméglicht die maschinelle Mischung von zwei Beigen von Loch- 
karten in eine bestimmte Folge. 

Multiplier. Fiihrt die in der Rechnung auftretenden arithmetischen Opera- 
tionen (Addition, Multiplikation, Subtraktion, Division) aus und locht die 
Resultate in Karten. Uns stand die Relaismaschine 602A und der elektronische 
Rechenlocher 604 zur Verfiigung. Beide Maschinen arbeiten mit festem Komma 
und besitzen interne Speicherméglichkeiten fiir 72 Dezimalziffern beim 602A 
bzw. 37 Dezimalziffern beim 604. Dazu kommen noch Zahler mit 30 Dezimal- 
ziffern beim 602A bzw. 13 Dezimalziffern beim 604, in denen Zahlen akkumu- 
liert werden kénnen. Bei diesen Maschinen kénnen im beschrankten Masse 
arithmetische Operationen zu einem kleinen Rechenprogramm kombiniert 
werden, das auf auswechselbaren Schaltbrettern mit Hilfe von Kabeln gesteckt 
wird. Die Rechengeschwindigkeit der 602A liegt im Mittel leicht iiber derjenigen 
der Z4 und ist vom Rechenprogramm abhangig. Die Rechengeschwindigkeit 
der elektronischen 604 hingegen ist so gross, dass die Zahl der verarbeiteten 


A 
7 
© 
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_Lochkarten nur von der Arbeitsgeschwindigkeit der Ablese- und Locheinheit 
-abhangt. Diese betragt 6000 Karten in der Stunde. _ 
4 Tabulator. Diese Maschine kann die gelochten Zwischen- und Endresultate 
in Tabellenform drucken und, wenn nétig, auch einzelne Additionen oder Sub- 
_ traktionen ausfiihren. . 

Mit Ausnahme der ersten zwei werden diese Lochkartenmaschinen durch 
_auswechselbare Schaltbretter gesteuert, auf denen mit Hilfe von Kabeln die 
_notigen Schaltungen gesteckt werden kénnen. Da die Erstellung jeder Schal-, 
_tung einigen Zeitaufwand -erfordert, muss die Rechnung so aufgebaut werden, 
_ dass sie aus méglichst gleichartigen Operationen besteht. Zudem soll jede Ope- 
ration auf eine grosse Zahl von Karten angewendet werden kénnen, damit sich 
der oft recht betrachtliche Zeitaufwand fiir den Aufbau der Schaltungen lohnt. 
_ Deshalb wird man kleine Nebenrechnungen gewohnlich von Hand ausfiihren. 
Aus den Charakteristiken der eigentlichen Recheneinheiten ersieht man, dass 
sich die IBM-Maschinen speziell fiir Rechnungen eignen, die sich aus Folgen 
von wenigen arithmetischen Operationen zusammensetzen, wobei die Grdéssen- 
-ordnungen der Zwischen- und Endresultate leicht tibersehbar sind. Bei der 

Flatterrechnung weist die Berechnung der Integrale, die in Formel (15) auf- 
treten, diese Eigenschaften auf. Auch das Erfordernis der grossen Anzahl von 
gleichartigen Ausdriicken ist schon bei wenigen Freiheitsgraden erfillt. Die 
Integranden bestehen aus Produkten von meist tabellarisch gegebenen Grdssen, 
so dass die Integrale sowieso nur numerisch integriert werden kénnen. Ein Teil 
dieser Gréssen ist gewohnlich nur auf einige Prozent genau bekannt. Deshalb 
geniigt fiir die numerische Integration die Verwendung der Simpsonschen Na- 
herungsformel, das heisst, die Integrale in (15) werden durch Summen ersetzt 
von der Form: 


b 
f Ej m(4) H(%, 0) dx = ——[Ejm(a) H(a, 0) + 4 Bj, (a + A) H(a + A, @) 


a 


+ 2Es_(a+ 2A) Hat 2.4, 0) +> 
+4E,,,(b — A) H(b — A, w) + E; (6) H(8, w)], 


wobei 2 die Zahl der gleich grossen Intervalle von der Lange A ist, die das 
Integrationsgebiet tiberdecken. Sieht man von den Integrationsfaktoren ab, so 
stellt die rechte Seite eine Art von Skalarprodukt dar. Da die IBM-Rechen- 
locher Zahler enthalten, in denen Produkte gleichzeitig gebildet und zueinander 
addiert werden kénnen, wahrend die Z4 diese Méglichkeit nicht besitzt, war die 
Verwendung dieser Lochkartenmaschinen aus zeitlichen und finanziellen Griin- 
den gegeben. Auf den detaillierten Aufbau dieser Rechnung kann jedoch im 
Rahmen dieser Arbeit nicht eingegangen werden. 
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Fiir die numerische Durchfithrung der Flatterrechnung nach den Methoden 
des § 3 sind die beschriebenen Lochkartenmaschinen jedoch nicht geeignet, da 


die Rechenprogramme zu kompliziert und zu lange sind. Infolge ihrer grossen — 
und einfachen Programmierungsméglichkeiten kann in diesem Falle hingegen : 


die Z4-Rechenmaschine mit Vorteil eingesetzt werden. Besonders das gleitende 
Komma erleichtert bei diesen komplizierten Rechnungen ausserordentlich die 
Arbeit mit der Maschine. Die beschraénkte Speicherkapazitat begrenzt aller- 
dings den Umfang der auf rationelle Weise noch zu bewdltigenden Flatter- 
determinanten. Praktisch stellen fiinfreihige Determinanten das Maximum dar, 
das die Maschine ohne spezielle Massnahmen verarbeiten kann. Der Zeitauf- 
wand fiir die Berechnung des charakteristischen Polynoms aus einer vierreihi- 
gen Determinante betragt beispielsweise etwa eine halbe Stunde. Die iterativen 
Methoden fiir die Nullstellenbestimmung kénnen sehr leicht auf dieser Maschine 
programmiert werden, da eine Wiederholung eines bestimmten Rechenpro- 
gramms durch das Zusammenkleben des Anfanges und des Endes des Film- 
streifens, auf dem die entsprechenden Befehle festgehalten sind, erreicht werden 
kann. Da keine Rechenmaschine fehlerlos rechnet, muss man der staéndigen 
Kontrolle der Rechnung spezielle Aufmerksamkeit widmen. Bei der Berechnung 
der charakteristischen Gleichung kann man die Zeilensummen der vom Faktor 
A freien Glieder als Rechenkontrolle mitfiihren, wahrend beim iterativen Auf- 
lésen dieser Gleichungen das Einsetzen der letzten Lésung in die urspriingliche 
Gleichung eine gute Kontrolle ergibt. 

Es ist vielleicht niitzlich, wenn hier auch einige Bemerkungen iiber die 
numerischen Eigenheiten der Methoden des § 3 aufgefiihrt werden. Beim Voet- 
terschen Verfahren ist zu beriicksichtigen, dass es versagt, falls eines der 
Glieder 6; , zufallig Null wird. In einem solchen Fall muss man von neuem be- 
ginnen und durch Umordnung der Zeilen und Kolonnen versuchen, das Auf- 
treten der Schwierigkeit zu vermeiden. Bei der iterativen Berechnung der Null- 
stellen sollte am Anfang die erste Nullstelle so genau wie méglich bestimmt 
werden, da wegen der Division des charakteristischen Polynoms durch die 
schon gerechneten Nullstellen allfallige Ungenauigkeiten die Berechnung der 
andern Nullstellen beeintrachtigen kénnen. Bei langsamen Rechenmaschinen 
kann man Zeit sparen, wenn man im Rechenprogramm die Moglichkeit einbaut, 
fiir jede Nullstelle eine geschatzte Lésung der Maschine mitzuteilen. Gewohnlich 
wird die gerade errechnete Nullstelle als Schatzung fiir die nachstfolgende ver- 
wendet, wobei dann unter Umstanden mehrere Iterationen benétigt werden, 
um nur in die Nahe der nichsten Nullstelle zu gelangen. Bei schnell rechnenden 
Maschinen lohnt sich allerdings ein solcher Eingriff nicht, da die Zeit, die 
gebraucht wird, um der Maschine die geschatzte Lésung einzugeben, grésser 
als die eingesparte Rechenzeit ware. Zu dem Verfahren fiir die direkte Bestim- 
mung der Nullstellen ist zu bemerken, dass die Methode versagt, falls die durch 
G,,(A, w) definierte Kurve die w-Achse nur beriihrt oder fiir zwei benachbarte 


ws 3 be 
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_w-Werte schneidet, zwischen denen kein berechneter Wert von G,(A, w) liegt, 
da man dann zu wenig Informationen fiir die graphische Interpolation der Null- 
stellen besitzt. Ein weiterer Nachteil dieser Methode liegt darin, dass keine ein- 
-fachen direkten Kontrollen der Rechnung existieren. Falls die Parameterwerte 
_q@ geniigend dicht liegen, gibt der Verlauf der Kurve G,,(A, w) eine gewisse Kon- 
trollméglichkeit, da G,,(A, w) in stetiger Weise von w abhangt. Auf diese Weise 
_kénnen jedoch nur grobe sporadische Fehler, hingegen keine systematischen 
_entdeckt werden. Die systematischen Fehler kommen nur dann zum Vorschein, 
-wenn die Lage der reellen Nullstellen, wie sie aus der graphischen Auftragung 
hervorgeht, im Widerspruch zu den Aussagen steht, die von der Sturmschen 
_Kette geliefert werden. Bei der Anwendung dieser Methode ist es also not- 
wendig, eine sorgfaltige Diskussion der Rechnung durchzufiihren, um Fehler 
in den Resultaten zu vermeiden. 
Mit Hilfe der beschriebenen Methoden wurde das Flattern der verschieden- 
-sten Teilsysteme untersucht. (Fliigel mit Querruder, Héhenleitwerk mit Ruder, 
Rumpfhinterteil mit Héhen- und Seitenleitwerk usw.) Durch den Einsatz der 
_programmegesteuerten Rechenmaschinen konnten diese Rechnungen in ver- 
haltnismassig kurzer Zeit und mit dem Einsatz von ein bis héchstens zwei 
Leuten durchgefiihrt werden. 
Der Verfasser méchte an dieser Stelle allen Herren, die ihm beim Zustande- 
kommen dieser Arbeit in irgendeiner Weise behilflich waren, seinen Dank aus- 
sprechen. 


Summary 


A set of three integral equations for the components of the relative displace- 
ment vector are used to describe the flutter oscillations of an aircraft. Two methods 
for the approximate solution are presented, one making use of the eigenvector- 
functions and eigenvalues of the problem without airforces, the other being 
RAUSCHER’S method of stationfunctions in a slightly generalised form. Then 
numerical methods for finding the zeros of the flutterdeterminants are given. 
The paper is concluded by a short description of actual numerical computations 
on automatic computing machines. 


(Eingegangen: 20. Dezember 1954. Veréffentlicht ausserhalb des Normalumfanges der Zeitschrift.) 
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Mathematical Analysis 
of a Simple Model for the Stripping Reaction 
By Res Josr4), Princeton, N. J., U.S. A-?) 


Introduction 


The mathematical problem analyzed in this papercorrespondstoanextremely 
simplified model of the stripping reaction. From the standpoint of the physicist, 
this investigation is, therefore, somewhat academic. Furthermore, we do not, 
in this paper, actually solve the problem, we only indicate the procedure which 
leads to the solution. This procedure, however, seems to us to be of sufficient 
interest in itself to merit publication. It shows, namely, that in this case, the 
wellknown method of N. WIENER and E. Horr can be successfully applied to a 
more complicated integral equation [compare (1.6)]. It is unknown to us how 
fruitful this generalization may be. Apart from somewhat trivial generalizations 
to systems of integral equations, no attempt has been made to enlarge the 
class of integral equations to which our method is applicable. 


1. Derivation of a Functional Equation 


We start with the following differential equation 


Get + et + ALE) (x) + 6(n) O8)] y+ Ey =O, | 
LES OG | (1.1) 
OD 46.510) 2 : 


(1.1) can be interpreted as the Schrédinger equation of two particles 
moving each in one dimension and acting on one another with a 6-function 
potential. The coordinates of the particles are (& + y) and (€ — 7). Furthermore, 
we assume that the particle with the coordinate (& — 7) is reflected by an 
infinite wall located at € = 7. (1.1), then arises by using the reflection principle. 


1) Most of the content of this paper grew out of discussions, which the author had with J. H. 
D. Jensen, J.M.Lurrincer, and T.H.BEeriin during their respective stays at the Institute for 
Advanced Study. The author’s only merit consists in having written down the results. 


*) Institute for Advanced Study. Now at Eidgenéssische Technische Hochschule, Ziirich, 
Switzerland. ; 


be 


; 

24 
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ys ‘ 

_ The symmetry condition 


ylE, n) + y(n, €) = 0 (1.2) 


‘is equivalent with the vanishing of y(&, 7) at the wall. 

We are interested in solutions of (1. 1) of which the incoming part corre- 
_ sponds to a bound particle. In order to get such a solution the energy E has to 
exceed the energy of two particles bound by a 6-potential. The condition 

A\2 


E=-(5) (1.3) 


03> 54) 
‘ 


* 


is easily verified. 
(1.1), (1.2), and the condition on the asymptotic behavior of y (with a 
suitable definition of outgoing waves) leads to the following integral equation: 


hit 
€, ») =—_| 4’ 
ven= J oe 


x [A® VEVE—£)2 + 9) — HP (VEVE+ (n—£))] OE), 
where 
@() = ylé,0); forE<0:ImVE>0; forE>0:VE>0. (1.5) 


If we make 7 zero in (1. 4) we get a homogeneous integral equation for ®(€). 
ge) = 4* / de’ [H® (VE |e — £'|) — H® (VEVe2 + &)] ®E). (1.6) 
0 


: Our main concern is to provide a method for solving this integral equation. 
The form of (1.6) suggests a Fourier transformation. For reasons, which will 
appear later, we are, however, forced to treat (1.6) (at least for the case E > 0) 
as a limiting case: let E have a small positive imaginary part: 
E=E,+ita («> 0) (1.7) 
and define VE such that ie 
Im [/E] > 0. (1.8) 
Finally, we will take the limit « > 0. We now impose on the solution ®(&) the 
condition @(é) = 0 (E+) (y>0) for >to. (1.9) 
As a consequence of (1.9) the integral in (1.6) becomes absolutely conver- 
gent, and furthermore, 


P(E) glive+el€ 50 foré > —oo and arbitrary e > 0. (1.10) 


| 
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Now the Fourier transformation can be performed. Following a procedure 


by N. Wiener and E. Horr?) we split the transform of ©() into half trans- 


forms: 
0° 0 


; 
’ 


’ 


tp) = [eM OE) a, yl) = — [ei OE as. (1.11) 


—0o 


According to (1.9) and (1.10) g(p) is regular in an open upper half-plane, 
which contains the point VE; y(P) in an open lower half-plane Im (pf) < Im [VE] : 
p(P) and »(p) have, therefore, a common strip of regularity. The Fourier trans- 
formation of (6) now leads to the functional equation 


o(6) — yb) = 7 [p) — 9(o)]. (1.12) 
where 
o(p) =VE— 2, Im[o]>0. (1. 13) 


(1.12) is derived in a #-domain such that # is in the common strip of regularity 


of y(f) and y(p) and w is in the half plane of regularity of y(p). By analytical 
continuation (1.12) holds everywhere. 
In parenthesis we mention the Fourier transform of (1.4) 


any 


ob) = 25 ape (Ph) — 9(60)] (1.14) 


where 


wtb ba) =f [ a dn e8l#*H i, n) (1.15) 


The asymptotic behavior of y(é, 7), which is the interesting quantity, is in a 
simple way related to the value of the square bracket in (1.14) for vanishing 
denominator p? + $3 — E. This quantity, however, is just 


I(P) = (2) — @(@) . (1.16) 


With the help of (1.12) and the function-theoretical requirements on 9() 
and y(f) we can infer more about the singularities of w(p) and »(p) 

(1) y() has no singularity at p=VE; 

(2) y(f) has at most a first order pole at p=+VE+ (A/2)?; 

(3) v(4) has no singularity at p= —VE,p = + VE + (4/2)? 


1x TircuMarsH, Fourier Integrals (Oxford University Press, London, 1937), Chapter XI, 
17, p. 339. 
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Branch cuts for @ = yz — p. 


PLL Common strip of regularity of g(p) and y(p). 


Figure 1 
The cut p-plane. 


With regard to (2), we can remark that () actually has the first-order 
poles mentioned, since asymptotically &(f) behaves as the superposition of two 
harmonic oscillations for § > + oe. 

Finally g(p) > 0 for |p| >co in Im[p] > Im VE + (A/2)? and y(p) >0 
for |p| >oo in Im[p]<Im/E. 


2. Derivation of a Difference-Equation from (1.12) 
(a) The Uniformization 


One of the difficulties which make (1.12) not transparent is the fact, that 
w(p) and, therefore, also p(p) and y() are multi-valued functions of #. It is not 
hard to convince oneself that g(p) and y(p) are, in general, infinitely many 
valued. We, therefore, try to uniformize alle these functions. It is clear that 
we can only hope for this with a transcendental uniformization of w(p). As such 
we choose 

p=VEsin= Z, wo =VE cos Z. (2.1) 


Now (1.12) reads in a somewhat sloppy notation 


WO) VO) = esa TIO Me +H (2.2) 
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or, making use of the abbreviation 


cos Co=+ ane 


2VE 


2z 2 


(0 < Re[{] <1, Im[g,] < 0), (2,3) 


cos. Zy(Z) = [cos Zs cal p(Z) + cos f,o(Z +1). (24) 


Next we have to translate the properties (1), (2), (3) of section 1 into proper- 
ties of the new functions g(Z) and (Z). Of special importance are the properties 
(1.1) and the first property (1.3). We normalize. the mapping (1.1) so that 


Z = 1 corresponds to VE. Z = — 1 to — VE of Figure 1. Then the fact that — 


y(p) has no singularity for p = /E reflects itself in a symmetry property for 
y(Z): 
g(1l+Z)=g(1-—Z). (2.5) 
Similarly 
y(—1 + Z) = y(—-1-—Z). (2.6) 


The rest of the discussion is quite simple and we indicate only the results: 

(1) p(*+7y) and y(% +7 y) tend to zero for | y| +00; 

(2) p(Z) can only have poles in the points + & + n(n = 0, —1, + 2...); 
it has no pole at + ¢,+ 1 and only a first order pole in 


+ to9—1,+ C4 & Oo +2, +0,4+3. 
(3) y(Z) has no pole in the points - f +&—-1,+6,-2. 


(b) The Difference Equation 


With (2.5) and (2.6) we are able to eliminate from (2.4) the function y(Z). 
Replacing Z by — 2 — Z in (2.4) and using (2.5) and (2.6) leads to 


—cos> Z y(Z) = [- cos Z — cos ca] w(Z + 4) + cos fy q(Z+3) (2.7) 
and addition of (2.4) and (2.7) gives 


cos, Z[p(Z + 4) — 9(2)] 


; | (2.8) 
+ cosy ol p(Z + 4) — 9(Z + 3) — 9(Z + 1) + o(Z)] =0. 


(2.8) is simplified by introducing the function /(Z) [compare (1.16)] which 
contains the relevant information: 


f(2Z) = (2) — 9(Z + 1) (2.9) 
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AVANT S Poe 


g According to (2.5) {(Z) has the symmetry 
1 
i+) 2)~0. an 


_ The equation for /(Z) can be written as 


cosy Cf (Z) + HZ +1) + (Z + 2)] | 
— cos Z[f(Z) + f(Z + 1)] + sin> Z[f(Z + 1) + NZ + 2)] 
= cos = Cd il2 ed) af + 2) 2 3) (250) 


— cost (Z+ 1)[f(Z+1) + (Z+2)] 


sin (2-4) (2 + 2) + fZ + 3)] 


2 


or else 


cos Co[f(Z) + HZ +1) + HZ + 2)]— cos Z[f(Z) + HZ + 1)] 
(2.12) 
+ sin> Z(f(Z+ 1) + f(Z + 2)]=2(Z). 
When x(Z) is a periodic function of period 1, it will turn out that 2(Z) = 0. 
To show this we first note that 


n(—Z) + x(Z) =0 (2.13) 


according to (2.12). 
From the property [a(1)] of the preceding section we have 


av +7 y) 
att a >0 for|y|>o. (2. 14) 


Using (2.4), (2.5), and (2.6), one finds furthermore, 
a(Z) = [cos y= cos Z| [ol [p(Z) — 9(Z + 2)]—siny Zy(Z—1) (2.15) 


and from the properties [a(2)] and [a(3)], one concludes that z(Z) has no pole 
at Z=-+ ¢,. This is, however, the only place where 2(Z) could have a pole 
according to its definition (2.12). Evidently, 2(Z) is a rational function of 
v = exp 21 Z which according to (2.13) has zeros at v= + 1 and no poles, 
because (2.14) excludes an infinity at v = 0 or v = oo. Therefore, (Z) vanishes 


ZAMP VI1/21 
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and we have 


F bol f(Z) )+ 42+ 1) + KZ + 2)]— cosy Z[KZ) + (Z + 1] 


+sin> Z[f(Z +1) + (Z + 2)]=0 
and still 


It is slightly more agreeable to transform (2. Ne to a normal form by the 
substitution 


f(x tiy)>0 for|y|>o. (2.17) 


h(Z) = [cos 5 by — cos Z| [4(2) 1) eee (2.18) 
Which has the inversion 
cost fy f(Z) = A(Z— 1) + A(Z). (2.19) 
The equation for h(Z) then reads 
h(Z +1) + A(Z—1) = ¢(Z) A(Z) , (2.20) 
cos (1/2) €5 — 2 cos(x/2) Z 
A es cos (7/2) €) — cos(x/2) Z (2.21) 
and (2.10) leads to 
h(Z) + h(—Z) = 0 (2.22) 
and (2.17) to 
LSA Ian i ead en (2.23) 


cos (7/2) (* + 7 y) 


In the next paragraph we need, however, a sharper condition for the beha- 
vior of A(x +7 y), which we are going to make plausible. Since 


C(x+ty) >—2 for|y] > co 


it is reasonable to assume (and can be proven) that A(x + 7 y) approaches in 
the limit y + + oo (or y > — oo), a solution of the equation 


Wea Lyte Bees lee ee (2. 24) 


The behavior of h(x + ty) for y oo (under reasonable restrictive assump- 
tions) is the following: 


hha) ide 2) gen (ja| + || >0), (2.29) 
where is an integer. Therefore, we expect 


h(Z) > (a+b Z)e™@"*V2 for y=Im (Z) ++ c. (2.26) 


SANA PL A ee 


Vol. VI, 1955 Mathematical Analysis of a Simple Model for the Stripping Reaction 323: 


_ According to (2.23), however, 7 = 0. Therefore, we have 


h(x+24)+>0 for|y| >. (2.27) 


3. The Uniqueness 


Let us now collect all the information about h(Z). 
(1) A(Z) satisfies the difference equation 


AZ + 1)4- kel) =O ZyR(Z), C(Z +4) = C(Z) . (3.1) 
(2) h(Z) is antisymmetric 
h(Z) + h(—Z) = 0. (3.2) 
(3) h(Z) is meromorphic and has (at most) first order poles in the points 
ety ek (r= 1,2, 3,4); 
(at most) second order poles in the points 
Gp cee (Ra= 5, 6,:7;.8) . 
These statements follow from [§ 2(a) (3)] and (2.9), (2.18). 
x h(x +iy)>0O for|y|>o. (3. 3) 


We will now show that these properties define 4(Z) up to a constant factor. 
If h,(Z) and h,(Z) are two nontrivial functions satisfying the above mentioned 
conditions then we first claim that their Wronskian vanishes: 


hy(Z) hy(Z + 1) 
W(Z) = =a). (3.4) 
h,(Z) ho(Z + 1) 
From (3.1) we have 
W(Z + 1) — Wiz) = 0 (3.5) 
from (3.2) 
W(-Z) + W(Z) = 0; (3.6) 


from (3.3) follows that W(Z) is meromorphic and has at most a first-order pole 
in the points + 6,+ »(n=0, +1, + 2,...). (3.4) finally tells us that 


W(x+iy)>0 for|y|>o. (3.7) 
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Therefore, W(Z) is a rational function of v = e**/4 with zeros at v = 0, v= +1, 
vy = oo and at most two poles. Such a function evidently vanishes. 


t 


| 
{ 
: 
5 | 
: 


Due to (3.4) A,(Z) and h,(Z) differ only by a factor, which is a rational 


function of U: 
hy(Z) = x(Z) ha(Z) . (3.8) 


If now z(Z) is not a constant, we can find a nontrivial linear combination of 


h,(Z) and h,(Z) with constant coefficients such that 


h(Z) = 4 h,(Z) + OX h(Z) (3.9) 
has at most 


a zero of first-order at Z=-+¢,, 
regular points at Lise tot ek, (= 1,2,3, 45 
a first-order pole at 42S ty he k= 5-6, 7.8). 
The function . 
h(Z) = h(Z + 4) — k(Z — 4), (3. 10) 
therefore, satisfies (1), (3), (4), instead of (2) we have 
A(Z) — h(—Z) =0. (3.11) 


Now we claim that the Wronskian 


| (3.12) 


vanishes. Obviously W(Z) is again a rational function of v which vanishes for 
v = 0 and v = oo is everywhere regular and hence zero. 
There is, therefore, a rational function of v, R(v) such that 


h(Z + 4) —h(Z — 4) = 2 Riv) W(Z). (3.13) 


On the other hand, one derives easily, that every solution of (3.1), hence also 
h(Z) satisfies the equation 


h(Z + 4) + h(Z — 4) =2 A(v) h(Z) , (3.14) 
where 
2 A(v) = ¢(Z) ¢(Z + 1) C(Z + 2) C(Z + 3) 
—[0(Z) (2 + 1) + (2 + 1) (2 + 2) + C(Z + 2) C(Z + 3) (3.15) 
+ €(2 + 3) €(Z)] +2. 


4 
a 
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; i otaing to (2.21) one finds 


v?@— 2av+4+1 


AS wo 2pvy 1’ 
a=—4 costS Co+ 8 cos* Carrel 3 (3. 16) 
p= 8 cost $5 — 8 cos? oti. | 
_ Addition and subtraction of (3.12) and (3.14) yield 
WZ +4) =(A+ R)K(Z), (Ss17) 
WZ +4) = =p WZ) (3.18) 
and these equations are only compatible if 
Revs 1. (3.19) 
In order that R be a rational function of v, « — f has to vanish, i.e., 
Afvy=1, (3.20) 
Ri(o)-= 0: (3.21) 
But in this case - 7 
h(Z + 4) =h(Z). (3-22) 


Since, however, (Z) has no singularities at Z = + €,+ k (k = 0, 1, 2, 3); it 
has no singularities at all and, due to (3.3), it vanishes. Thus we have reached 
a contradiction, and the z(Z) in (3.8) is a constant. 


4. Final Remarks 


After the preliminaries of the preceding sections, the reader could reasonably 
well expect us to solve the difference equation (2.20), (2. 21). This equation 
can actually be solved. The theory of a class of difference equations containing 
the one under discussion, has been given recently!). This theory contains a 
construction for the solution, which in our case may be manageable, though 
tedious. We will not attempt here to push it through. 

There is, however, one case for which the solution is trivial. This is the 
case where the auxiliary equation (3.14), (3.16) 


H(Z +4) +H(Z—4) =2 A(0) H(Z) (4.1) 


1) R. Jost, Comm. Math. Helv. 28, 173 (1954). 
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reduces to an equation with constant coefficients: 


4 


A(v)=1, cos*F y= 5. 


One verifies easily that every solution H(Z) of (4.1) in the general case gives rise — 
to the following solution of (2.20), (2.21): 


h(Z) = 
=H(Z+3) +0 (243) H(Z +2) +[6(2 + 3) o(2 + 2) 1 (Zz + 1) 
+[6(Z + 3) ¢(Z + 2) ¢(Z +1) — ¢(Z + 3) —¢(Z +1) ] AZ) 

+ H(Z — 3) + ¢(Z2 +1) H(Z — 2) +[¢(2 + 1) €(Z + 2) -1] H(Z — 1). 


(4.3) 


If one applies (4.3) to the special case (4. 2), where the auxiliary equation 
(4.1) has the fundamental system of solutions H,(Z) = 1, H,(Z) = Z one gets 
the following fundamental system of (2.20): 

1 + 3 cos (x/2) Z 


= 4.4 
me) [4 — 3 sin®(x/2) Z] [2 + 3 cos (x/2) Z] ° ee 


h(2) = [1 + 3 cos (x/2) Z] Z + V3 sin (x/2) Z (4.5) 
[4 — 3 sin? (x/2) Z] [2 + |/3 cos (x/2) Z] 
Evidently ,(Z) satisfies all the properties of Section 3 which are required 
for h(Z) and, therefore, solves the problem. 


Zusammenfassung 


Die Integralgleichung (1.6), die aus einem sehr vereinfachten Modell der 
«stripping reaction» hergeleitet ist, wird diskutiert. Es ist bemerkenswert, dass 
die bekannte Methode von N. WrENER und E. Hopr auf diese komplizierte Inte- 
gralgleichung angewendet werden kann. Die entstehende Funktionalgleichung 
(1.12) wird im §2 auf die lésbare Differenzengleichung (2.20) reduziert. § 3 
behandelt die Eindeutigkeit der Lésung und § 4 die explizite Lésung in einem 
Spezialfall. Die explizite Lésung des allgemeinen Falls wird nicht angegeben. 


(Received: June 24, 1954.) 
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= Kurze Mitteilungen — Brief Reports - Communications bréves 


Eine neue Berechnungsmethode der quadratischen Regelflache 


7 Von Kraus ANKE, Erlangen?) 


: Einleitung 


Bei der Untersuchung linearer stetiger Regelungen wird vielfach als Giite- 
_ mass fiir die dynamische Regelabweichung ~(t) die 


co 


« quadratische Regelflache » - w(t) dt 
0 


verwendet; dabei ist x(¢) die (durch gegebene Anfangsbedingungen festgelegte) 
Lésung einer gewdhnlichen linearen Differentialgleichung mit konstanten Koef- 
fizienten. Der Verlauf der Regelabweichung 4(¢) wird dann als besonders giinstig 

- angesehen, wenn durch geeignete Wahl der verfiigbaren Koeffizienten die quadra- 
tische Regelflache zu einem Minimum gemacht werden konnte [1], [2], [4]?). 

Zur Berechnung der quadratischen Regelflache scheidet der Weg iiber die 
explizite Lésung x(t) aus, da eben wegen der verfiigbaren Koeffizienten diese 
Lésung im allgemeinen gar nicht explizit angegeben werden kann. 

Bei Benutzung der Laplace-Transformation gelingt es, zu einem geschlossenen 
Ausdruck fiir die quadratische Regelflache zu kommen ; der Weg fiihrt entweder 

_ tiber den komplexen Faltungssatz [1] oder einen Algorithmus [3]. 

In der vorliegenden Arbeit soll gezeigt werden, dass die Berechnung der qua- 
dratischen Regelflache im Zeitbereich auch ohne explizite Kenntnis der Lésung 
«(t) gelingt. ® 

Wir formulieren das Problem genauer: Es sei #(¢) die durch die Differential- 


gleichung 
2 a" x arty dx 
a de a cs Sal al (1) 
mit den Anfangsbedingungen 
dx hes 
as eee ee 2 2 
#(0) = qo dt t=0 Gir++e5 dt} a0 In-1 ( ) 


bestimmte Funktion der Zeit. Wir setzen voraus, dass die Koeffizienten 4, 4, 
..., Gm der Hurwitz-Bedingung gentigen. Dann gilt 


dk x% 
im —— = ure Ol Ue 2 ices 3 
pe o 
und es existiert das Integral s 


fafa ae. 


0 


1) Siemens-Schuckertwerke-Aktiengesellschaft. 
2) Die Ziffern in eckigen Klammern [] verweisen auf das Liter 


Arbeit. 


aturverzeichnis am Schluss der 
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Berechnung von jJ, in einem Beispiel 


Wir demonstrieren das Verfahren zur Berechnung von J, zunachst am Beispiel 
n= 3. 
Die Differentialgleichung fiir die Funktion *%(¢) 


a®x a*4 dx 
ae Mr NRC 


mit den Anfangsbedingungen 


dx d*x 
BO) Toss ee 4 evan arenes 


und dem Verhalten im Unendlichen 
hmv = OF line == 6 
t>00 


multiplizieren wir nacheinander mit x, dx/dt und d?x/dt?; es entstehen so die 
Gleichungen 


aex a*x ax 
CD ae we rg a + Oy ae + dy & % = 0, 
7 a2x ax . AG. ax dx dx | ax ; rn 
3 dt gn We otto gpa nae” Ogee Cg Vos Mgr ee gee (4) 
d'x dx Cian Gin ax | d*% 


. { 5 ° = 
oUt ide ont die a fe i gee ih Oe 


Wir integrieren diese Gleichungen von Null bis Unendlich und beachten dabei 
die Anfangsbedingungen, das Verhalten im Unendlichen und die folgenden, durch 
partielle Integration entstehenden Relationen: 


[ee [oe 
a> x di dix "Ge Oe di : ie 
a” ee der ap Cee Se 
0 0 
pd? [oe oo 
a ax ax \2 ax \2 
iaere at” es, oe fo 41 ie i 
0 0 0 0 
[oe) 
ax & 
[Gru--4 q » 
0 
ee) co co 
lan = Gtx «ae “* a2x oo. ag 7 
ais ag > ay dt? red teee (Gr) e 
0 0 0 
[o.e} [oe) 
(HEA (alee il da ax Ne eae 
|e ie Bite IG faa ane | 
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or 


NER 


‘Mit den Abkiirzungen | 


= 


PS oo ; co co 7 
s f= 2f x ae, j= 2] (G) ae, h=2f (Ga) a 
0 0 0 
; entsteht schliesslich aus (4) das lineare Gleichungssystem 
Ay Jo + 2 Sy ia + a4 Qip + 4, Qo + 43 O30 (= Qo) » | 
ay Jy + a3 Jz = — 4% Qor — a, Qo1 — 43 O31 (= —Q)) » (5) 
Ay Jy + 42 Jp = +4 Qo2 + 41 Qie + ds Q32 (= Qs) - | 


; mit den Unbekannten J, J,, J, und den Abkiirzungen 
Q5=1= % > Qo = Vox = 2% M1 > Q.=Qn = Wi,» 
Q39 = 2409—-T> V1 =2n 4; Qs. = 93 - 


Fiir den allgemeinen Fall ist es vorteilhaft, die Matrixschreibweise zu ver- 
wenden; wir gehen daher an dieser Stelle dazu iiber und setzen 


a 

| Ay Ay O ie | ds 0 Qiw Qs Qo 

A=j||0 4 4;,||, @= i ; P= » Q=|/-Qo1 9 —-Qn —Qs1 

0 ay ay Ss | Jo Qo2 Qin 0 Qs2 | 
’ 3 
Das Gleichungssystem (5) schreiben wir damit 
AJ=Qa. (6) 
Wir lésen (6) nach J auf und erhalten 

Ja47Q a. (7) 


_ Fiir die uns hier durch die Aufgabenstellung allein interessierende Unbekannte Jy 


ergibt sich aus (7) 
ae Qo , a (az Qi + 43 Qe) 


ao Ay (Ay Fy — Ay A) 


: Das System (6) ist genau dann eindeutig auflésbar, wenn die Determinante 
| A | der Matrix A von Null verschieden ist; dies wird aber durch die Erfiillung 
der Hurwitz-Bedingung gew4hrleistet, denn | A | ist die Hurwitz-Determinante. 


Allgemeiner Fall 


Fiir den allgemeinen Fall der Differentialgleichung (1) mit den Anfangsbe- 
dingungen (2) und dem Verhalten im Unendlichen (3) verfahren wir analog. Wir 
multiplizieren also (1) der Reihe nach mit 

a’ x 
i 

Die entstehenden -Gleichungen integrieren wir von Null bis Unendlich unter 

Beachtung der Anfangsbedingungen, des Verhaltens im Unendlichen und der 


y=0,1,...,”—1. 


ad 
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folgenden, durch partielle Integration entstehenden Relationen: 
P dx Gx ee \ —Quv falls w+ v = 1 (2) : 
di dt ~ (Que + (—1)" Jape falls w+ » = 0 (2) 4 
0 . 
Btyv=0,1,...4,2%— 1; 
mit den Abkiirzungen 
(u—v—1)/2 
2 a Se cae Tu—o Wyto-1 a pts Geutv—1yr2 lad ae! (2) 
o=1 : 
One Ore ea 
(u—»)/2 . 5 
(u> ») 2 ah 1p ee keg er w+ vy =0 (2) 
o=1 
foe) ee 
a* x 
sen apa f (SE 
0 
Mit den Matrizen 
ONC re | Ay || Jo || 
0 a a as i a | St | 
0 2 a, % | el - |] 
a5 y; @=l]- Is F=f. |p 
| i I]. 1] 
I} \| || | 
| ae Sah et 
0. Gn—1 || || %n || Jn—a 
0 Qi Qao --- Qno | 
—Q 0 = Osi eareny =O. 
02 Ore SY pip ne 
Q => 
| 
: | 
(= Ajeet Qon-1 As (> Blea: oe, n—1 


schreiben wir das aus (1) durch die angegebenen Operationen entstehende Glei- 
chungssystem 
ea) = Ola 


J=A7Qa. (8) 


mit der Lésung 


Damit ist die gestellte Aufgabe auf die Multiplikation bekannter Matrizen zuriick- 
gefiihrt und wird als gelést angesehen. 


Minimalproblem 


Es ist nicht angebracht, fiir den allgemeinen Ausdruck von Jo, Wie er sich aus (8) 
ergibt, das Minimum aufzusuchen. Die entstehenden Resultate wirden wegen 
der Allgemeinheit der Anfangsbedingungen und der jeweils moglichen verschie- 
denen Anzahl der freien Koeffizienten a, zu uniibersichtlich werden. 


“7 


eR INES 


< 
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~ 


Ne, 


__Fiir einen zur Orientierung bei den Anwendungen recht gut brauchbaren Fall 
“geben wir jedoch die Minimalbedingungen an. 

__ Wir setzen zunachst go = 1, 9) = dg = *** = Yn_1 = 0 und, ohne Beschrankung 
der Allgemeinheit, a) = a, = 1; sodann sehen wir alle iibrigen Koeffizienten 
1 4g,..-, Ay_, als verfiigbar an. Die Durchfiihrung der notwendigen Rechnungen 
ergibt, dass unter den genannten Voraussetzungen J) ein Minimum wird, wenn 


: Pa | 


Ao ay as a3 a, as a ay ag 
a 1 1 
2 1 all al 
ee 1 Zz, i! 1 
n= 4 ie 2 3 df 1 
=" 5 1 3 5 4 al al 
n= 6 1 3 6 a 3 iT it 
N= 7 al 4 6 10 3 6 of 1 
n= 8 1 4 10 10 i ys 6 a 1 1 


nd mae ais We) ai wl pi rar le Tele le oe, 2 tte. i @ © coe 6 0. 6.0%. .0 6) 8 80 'b. 0. 0 y780.6ie.\e) te) #18) 4! \8! (8) .911'6! (0): 


fiir Differentialgleichungen, die bei den genannten Anfangsbedingungen die 
‘Lésung mit minimaler quadratischer Regelflache besitzen. 
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Summary 


When analyzing linear and steady controls the so-called quadratic control 
area is frequently used as a measure of quality for the dynamical deviation of the 
control. All familiar methods used for calculating the quadratic control area are 
based on the Laplace transform. A method is presented by means of which the 
quadratic control area can be calculated in the time domain without explicit 
knowledge of the control’s dynamical deviation itself. 


(Eingegangen: 6. Januar 1955.) 
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A New Derivation and Interpretation of Yule’s ‘Characteristic’ K 


By Gustav HERDAN, Bristol (England) ?) 


YuLe’s ‘Characteristic’ K of the word-frequency distribution of a linguistic 
text is derived under the assumption that the occurence of a word in such texts 
was governed by a law of chance (the Poisson law)?). This assumption, and with 
it the use of K as a characteristic of the text, has been attacked by linguists), 
without foundation, as the writer believes. : 

However, the constant K can be derived without such an assumption, which 
has not only the advantage of obviating adverse criticism of the kind referred to 
above, but of showing K to be an easily interpretable, useful and interesting 
characteristic of a linguistic text. 

The number of words that occur X times in the sample from a given text is in 
statistical terms called the frequency of X (X being the ‘variable’) and the 
manner in which the frequencies are distributed over the scale of X is called the 
frequency distribution. Such distributions are peculiar in that 


the arithmetic mean M, = tee 
a 2 lage 3 
the standard deviation o,, = N M2. 
and 
the coefficient of variation v, = pe (1) 
x 


all vary with the size N = 7, + f, + ---+ fy, of the sample, increasing as N in- 
creases, but each doing so at a different rate. 

However, the quantity v,/ N is found to remain sensibly constant for samples 
of any size. The reason is that it is no longer dependent upon N, as can be seen 
when it is written explicitly as 


ve _ oa/N _ Ef X'N*— MBN _ SfpXYN? 1 Sf,X® 1 


N M2 M2 Te Teo XING | ON ef ee (2) 
since N cancels out in the first term on the right and 1/N becomes for great N 
negligibly small. 

Expressing (2) by the 1st and 2nd moment of the variable about zero, S,= Ds fy 
and S,= >’ }, X#, we write the supposedly constant characteristic as 


ug Sg 1 


Vo 2p 8) 
and for great N as 
S 
Na pee 
z ae 136 = F (3a) 
1) University of Bristol. 
‘ . Upny Yutr, A Statistical Study of Vocabulary (Cambridge, 1948). 


.S.C. Ross, Philological Probability Problems, J. Roy. Statist. Soc. [B], 12, Nr. 1,39 (1950). 
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_ which, but for the factor 104 which YULE used so as to avoid very small decimals, 
is YULE’s ‘characteristic’!), K = 104 (S,/S}— 1/S,) for great samples when 1/S, 
_ becomes negligibly small. 

¥ - We have thus derived the ‘Characteristic’ without any assumption about the 
_ variate being governed by the Poisson law, and indeed without any assumption 
_ about a stochastic process, as implied in YULE’s derivation, for the generation of 
_ word occurrences. 

a Moreover, the ‘Characteristic’ when written in the form 


o6lVN_ 
My 


Da) oi | | 


_ is seen to represent the coefficient of variation of a mean, or the relative fluctua- 
tion of a mean, as distinct from that of a single value, o,/M,. We shall distinguish 
_ it from the latter by the symbol v,, and understand by it the coefficient of varia- 
tion of the sampling distribution of means. 

: To test the deduction that v,, is a constant for a linguistic text, a Russian 
word count was used. JossErson of Wayne University, Detroit, has made a 
- complete count of the approximately 30,000 words in PusHKIN’s ‘ The Captain’s 
_ Daughter’). He then made a count for a sample of approximately 10,000 words 
; composed of all the words in 16 two-page groups spread uniformly throughout 
_ the book; finally he withdrew at random from the entire book two independent 
_ samples of approximately 5,000 words each. The exact number of occurrences are 
_ given below. For each of these word counts, I calculated the quantities M,, o,, 
Uy, Um as given in the following table. 


4 ° 


Occur- | Vocabulary M 
rences 


0-00603 
0-00656 
0-00663 
0.00673 


0-00583 
0-00610 
0-00605 
0-00606 


We note that whereas neither M,, nor ¢, nor v, can be used as a characteristic 
of the text, since they vary with the number of occurences and with the size of 
vocabulary (the sample size), rising as the latter increases, v2, remains constant as 
predicted and is sensibly equal to Ew a... 

This brings out an important property of word counts, and indeed of linguistic 
texts in this respect: such texts are characterised by the ratio 


=. +. = 


Standard Error of the Mean 
Mean 


being constant throughout the text irrespective of the sample size, and also 
irrespective of the mean, the standard deviation and the coefficient of variation, 
each taken by itself. Graphically this means that the standard error of the mean 


o/VN plotted against the mean in a sample of size N must be a straight line, the 


1) G. Upny Yue, A Statistical Study of Vocabulary (Cambridge, 1943). 
2) H. J. Josserson, The Russian Word Count (Detroit, U.S.A., 1953). 


> 7 
a 
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slope of which is v,,. Applied to our data, the relative fluctuation of the average } 

frequency of use of a word by PUSHKIN was about 8% of the mean frequency in| 

word aggregates, no matter what the size of these aggregates. 3} 
In general, we can say that a style with small relative fluctuation is charac 
terised by words being rare whose frequency expressed as deviation from the 
mean frequency represents an appreciable fraction of the latter, and a style with) 
a large relative fluctuation by such words being numerous. This implies that a 
particular style may be characterised by a constant relation between uniformity’ 
and diversity in the number of iterations of the items of vocabulary. 


: 
Zusammenfassung | 
| 
Worthaufigkeitsverteilungen sind dadurch charakterisiert, dass Mittelwerte, , 
Streuungen und Variabilitatskoeffizienten nicht konstant sind, sondern mit der’ 
Stichprobengrésse wachsen, aber jedes in verschiedener Weise. Es wird nun ge-: 
zeigt, dass der Variabilitaétskoeffizient des arithmetischen Mittels unabhangig ist: 
von der Stichprobengrésse und deshalb geeignet erscheint, als statistischer Para-: 
meter solcher Verteilungen zu dienen. Diese Grésse ist im grossen und ganzen. 
identisch mit G.UpNy YuLes ‘Characteristic’ K, wird aber hier abgeleitet ohne: 
die Annahme einer vollstandig zufallsmassigen Wahl der Worte, wie sie Y ULE} 
statuierte. 
Vom Standpunkt der theoretischen Statistik aus gesehen, ist es interessant, | 
dass es Verteilungen gibt, in denen nur der Quotient 


mittlerer quadratischer Fehler des Mittels 
Mittelwert 


konstant ist und unabhangig vom Werte des mittleren quadratischen Fehlers und 
vom Werte des Mittels. 


(Received: February 28, 1955.) 


On the Effect of Lubricant Inertia in Hydrodynamic Lubrication 


By FLETCHER OsTERLE and EpWarp SAIBEL, Pittsburgh, Pennsylvania, U.S.A.) 


The Reynolds equation of hydrodynamic theory, modified to take lubricant © 
inertia into approximate account, is applied to the steady-state operation of the 
inclined slider-bearing to determine the effect of lubricant inertia on load capa- | 
city. A simple relationship results, relating this ‘inertial’ load capacity to the 
physical properties of the lubricant and the operating characteristics of the slider- 
bearing. It is found that the inertia effect increases the load capacity by an | 
amount which can be significant in the laminar regime. 


Introduction 


The laminar hydrodynamic lubrication problem is governed by a differential 
equation which relates the pressure and viscous forces acting on a particle of the 


1) Departments of Mechanical Engineering and Mathematics, Carnegie Institute of Technology. 
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fan 


i 


_ lubricant to the lubricant inertia in accordance with Newron’s second law of 


- motion. In solving this equation for the steady-state load-supporting capacity of 
_ hydrodynamically lubricated bearings the inertia term is always neglected on 


ey 


‘the grounds that under normal operating conditions it is quite small in compa- 


_ rison with the maximum values of the pressure and viscous force terms. Therefore, 


me eae 


: 
of 
i 
; 


at the present time it is impossible to evaluate quantitatively the effect of 


- inertia on load capacity under a given set of operating conditions. Such an 
evaluation, if available, would perform the useful function of establishing the 


range of operating conditions under which neglect of the inertia effect is a permis- 


_ sable approximation. It is the purpose of this paper to determine quantitatively 
~ the effect of lubricant inertia on the load capacity of the slider-bearing by a 


method based on the approximation procedure suggested by SLEZKIN and TarRG 


_ (1946). From the results of this investigation the range of validity of the negligible 


- inertia assumption can be readily determined. 


Bs 
a 


Theory 
The differential equation governing the steady-state laminar hydrodynamic 
lubrication of a slider-bearing without side-flow is 
0?u dp Ou Ou 
a ASE pees 1 
Voy? ax efuse ros}, (1) 


where 
p(x) is the lubricant pressure, 
wis the lubricant viscosity, 
o is the lubricant density, 
u(«, y) is the lubricant velocity in the » direction, and 
v(x, y) is the lubricant velocity in the y direction, Figure. 


Jy 
B 


ie 


Slider-bearing configuration. 


In this paper the lubricant viscosity and density will be assumed constant. The 

first term in equation (1) originates from the viscous force, the second term from 

the pressure force, and the third term from the lubricant inertia. ot 
In addition to equation (1), the lubricant must also satisfy the continuity 


equation which for an incompressible fluid is 


Qa thay 7 
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and the following boundary conditions ; 
u(x,0)=—U, wu(x,h)=0, v(¥,0)=0, v(x, hk) =0, (3a) 
p(0)= 0, p(B)=0, (3b) 


where U, h, and B are as shown in the Figure. : 

An exact solution of equations (1) and (2) for the lubricant pressure subject 
to the boundary conditions (3) is extremely difficult. However, owing to the 
relative smallness of the inertia term, an approximate solution would seem to be 
justifiable. SLEzKIn and Tare (1946) suggest averaging the inertia over the film 
thickness. If this is done equation (1) becomes 


0?u | 
bra aet “t/l uae to ae) ay (4) 


The right hand side of equation (4) is thus a function of ¥ alone. Equations (4) 
and (2) together with the boundary conditions on velocity (3a) are then readily 
solvable for the pressure gradient yielding 


db 6uU(H—h) , 9U* Jk . | (HA) (HE — h)2 Wh’ 
ie he w adde le Sh 7 Sey 


where H is the film thickness at which the right hand side of equation (4) vanishes, 
and the primes indicate differentiation with respect to x. The boundary conditions 
on pressure will be used to determine the constant H and the constant (C) which 
will arise when equation (5) is integrated for the pressure. Equation (5) is a 
reduced form of equation (7) in the Slezkin and Targ paper which is derived in 
somewhat more detail in that reference. 

To proceed from equation (5) to the pressure distribution, the film-thickness 
variation must be specified. In this paper we will consider the film-thickness 
variation of exponential type employed by CHARNES and SaIBEL (1952) and 
shown by them to approximate very closely the plane slider-bearing. This film- 
thickness variation can be written 


hah (6) 


where hmi_ is the minimum film-thickness (henceforth referred to as h,,) and 
is the ratio of minimum to maximum film-thickness. From equation (6) we see that 


,— */B 
min * , 


1 
is > h. (7) 
Therefore equation (5) becomes 
dP (H—h) hi, In 7 Hy2 
CP ee hs 4 30 {2 (5) et 4, (8) 
where 
Bie 
ait (Syn Ueney (9) 
Bs 
ig Cling thi 
A Gu Bl oo 6B Re (11) 


with Ie the Reynolds number based on the minimum film-thickness. 


Nia a 


4, 
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With the Reynolds number set equal to zero, equation (8) reduces to the familiar 
equation (which neglects inertia) as we would expect. If we let P, be the solution 
of the Reynolds equation and H, be the value of H corresponding to this case, 
"we can write for the solution of equation (8) 


Pee Po Ae, (12) 


UE NERA te 


‘ with 


ig = Ar (13) 


_ where the second terms on the right of each equation are small corrections to 
account for inertia. Substituting equations (12) and (13) into equation (8) and 
- neglecting small quantities we obtain 


dPy _ (Ho —h) hm 


and 
oP Ht hes Inv ee 2 
| dX We paves eae - 
_ Equation (14) integrates to 
5 3X 2X 
: P, Hy + A ae. (16) 


“The  sinr 2inr 


en 


- The constants C, and H, can be determined from the boundary conditions on 

- pressure (3b), namely that P, vanish when X equals either zero or one. Thus 
Se 

1 Ay 1 3 (1 — r?) (17) 


Pie oe One Hy = Ibm 7 = 9) 


C= 


Now equation (16) becomes 


(18) 


Va dy 9 maaan 


Po= "Fin y ,1-%7 1— 7 


Equation (18) expresses the pressure distribution obtained by neglecting the 


- inertia. ; 
Treating equation (15) the same as equation (14) we obtain 


27\Inv (/1— 77 \2 2inyr (beep hee 
Oaks a) 15 i‘ tt. oie 


Equation (19) expresses the pressure distribution which must be added to Py 
in accordance with equation (12) to correct for the effect of inertia. 


The load-supporting capacity per unit area is the average pressure, defined by 


B 
| p=[p dx (20 
, 0 
or in dimensionless form 
1 
P=/[Pax, (21) 
0 
where from equation (9) is 
) pelle (22) 
re ee 


ZAMP VI1/22 
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Substituting equations (18) and (19) into equation (21) we obtain 


— » C= (lait= Dities, ie ea (23) 
y 2Iny 1— 73 1— 7 
and / 
a 27 Inv (1— 7?\2 — Inv een 24 
ni Ee (Lanne Beat agemn al 


respectively. Now from equation (12) we have that 


P=P,+A1P, (25) 
or i 

p P, 

Po Po 


which expresses the ratio of the average pressure to the average pressure neglect- 
ing inertia. 


Discussion and Conclusions 


By equations (11), (23), and (24), equation (26) can be written 


5 o U hi, 
a at euene (27) 


where 7 has been written for P,/P, and is recognized as being a function only of ¢. 
It is given by 
1 — (v3 — 1)/3 Inv 


27 Inv (1 — v?\2 2 (iy) iia 1— 78 
ierreg, (=) 15 (v?— 1) ) 1— (#8 —1)/3Inr 1— (7®— 1)/2Inr 
1— 8 Al 1— 7? 
(28) 
Representative values of this function are tabulated in Table I. 
Table I 

Y é 

1/3 0-6740 

1/2 0-2099 

2/3 0-0727 

1 0 


It is apparent from equation (27) that the inertia effect increases the load 
capacity. 
In terms of the Reynolds number, equation (27) can be written 


—=1+/—" Re. (29) 


If we consider a slider-bearing with a minimum to maximum film-thickness ratio 
of 1/2 and a minimum film-thickness to slider length ratio of 0-001, the inertia 


a 


4 
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_ effect will increase the load capacity one percent at a Reynolds number of about 
_ 300. The critical Reynolds number for parallel flow (i.e., y= 1) is about 2000. 
_ Thus is would appear that the inertia effect can be significant in the laminar 
_ hydrodynamic lubrication problem. 

2 The effect of inertia on load capacity can be expressed in terms of fewer 
_ variables if equation (26) is written 


F 


om S539 2 
eats, 2. (30) 


Po Po 


where & has been written for P, and is a function only of 7. Representative values 
; of k are tabulated in Table IT. 


L\ 


Table II 
Y k 
1/3 0:01777 
1/2 0:00568 
2/3 0-:00160 
1 0 


It is apparent from equation (30) that the inertia effect is most pronounced in 

' the case of high speed, lightly loaded bearings. Further examination of equation 
(30) reveals that under normal operating conditions the inertia effect is indeed 
negligible. For example, with an ¢ of 1/3, a lubricant of 0-90 specific gravity, an 

_ average pressure of 100 lbs. per sq. in., and a bearing velocity of 500 in. per sec., 
the load capacity is increased only about 1/3 percent by the inertia effect. 
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Zusammenfassung 


Die Reynoldsche Gleichung der hydrodynamischen Schmierungstheorie, unter 
' ungefahrer Beriicksichtigung der Trigheit des Schmiermittels, wird fiir ein 
schiefes Gleitschuhlager im stationaren Zustand angewendet, um den Einfluss 
der Tragheit des Schmiermittels auf die Belastungskapazitat zu ermitteln. Man 
findet einen einfachen Zusammenhang sowohl zwischen der Tragheitswirkung 
und den physikalischen Eigenschaften des Schmiermittels wie den Bedingungen, 
welche das Arbeiten des Gleitschuhlagers spezifizieren. Man stellt fest, dass der 
Einfluss der Tragheit die Lastkapazitat um einen Betrag erhoht, der wichtig sein 
kann in der laminaren Region. 


(Received: July 27, 1954.) 
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Tasgung iiber elektronische Rechenmaschinen und Informations - 
verarbeitung in Darmstadt vom 25. bis 27. Oktober 1955 


Die Gesellschaft fiir Angewandte Mathematik und Mechanik (GAMM) und 
die Nachrichtentechnische Gesellschaft (NTG) im Verband Deutscher Elektro- 
techniker (VDE) planen eine Fachtagung «Elektronische Rechenmaschinen und 
Informationsverarbeitung». Sie ist in Aussicht genommen fiir Dienstag, 25., 
Mittwoch, 26., Donnerstag, 27. Oktober 1955, in Darmstadt, Technische Hoch- 
schule. Behandelt werden sollen die logisch-mathematische und die technisch- 
konstruktive Seite der Maschinen, die Ausfiihrungsformen und die Anwendungen. 

Interessenten sind gebeten, ihre Teilnahme spatestens bis Donnerstag, 15. 
September 1955, an den Ortlichen Tagungsleiter, Prof. Dr. A.Walther, Institut 
fiir Praktische Mathematik der Technischen Hochschule Darmstadt, anzuzeigen. 
Anmeldungen fiir Vortrage mit héchstens 15 Minuten Sprechzeit und in einer der 
Tagungssprachen, Deutsch, Englisch, Franzésisch, sind bis Sonnabend, 6. August 
1955, an den Tagungsleiter zu richten; sie miissen, um in Erwagung gezogen zu 
werden, von einer kurzen Inhaltsangabe des Vortrags von héchstens 20 Schreib- 
maschinenzeilen Lange (etwa 1500 Buchstaben) begleitet sein. A. WALTHER 


Errata 


Correction to the paper Note on the Boundary Layer on a Rotating 
Sphere. By Swami Dayar Nicam, Wageningen, Netherlands?). 


On page 153 NiGAm?) the functions F,, F;, F;, G,, G3, G; should read as follows: 


2 


2 


Fy=as(1—s)? (1+ 2s) — s? (1 — s)*, -G, = (1/2) (2+ s) (1 — s)? 
Fy = bs (i —«s)2 (Ll + 2\s), Gg =cs(1—s)?(1+ 2s) 
Fi=ds(1—s)? (1+ 2s), Gs; =es(1—s)?(1+ 2s) 


Numerical integration of the equations gives 
@ = 1:518, b= 0-373, c = 0-426, d = 0-253, e = 0-195, and 6 = 2-794 


The values of H,, H;, H; and other functions can be calculated as indicated in the 
paper. The general conclusions of the paper remain unaltered. 


2nd correction to the paper Boundary Layer on Rotating Spheroids. 


By BuaskKar SapAsuiv Fapnis, Kharagpur, India‘), 
On page 158 Fapnis*) equations (3.15) — (3.18) should read 


F=as(t—9*(1+25)— (>) st— 9) (3.15) 


1) Ship Model Basin. 

2) S.D. Nicam, Note on the Boundary Layer on a Rotating Sphere, ZAMP 5, 151 (1954). 
3) Department of Applied Mathematics. Indian Institute of Technology. 

4) B.S. Fapnis, Boundary Layer on Rotating Spheroids, ZAMP 5, 156 (1954). 


4. 

“ 
= 
2: 
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(aS (3) (l—s)*(24+ 5) (3. 16) 
; Rate@—s(re2s) = st (= (3.17) 
Bath aes aie ais) (3.18) 


| The corrected values of the constants are 
@ = 1:5183, b = 0°4158, c = 0:3304, 6 = 2:794 
: Page 159 the functions H,, H, are 
f —H, = 4:2822 s?— 7-271 53+ 4-5432 st — 0-9789 s°, 
—H, = 2:5358 s? — 0:1091 s* — 3-6402 s* + 1-9632 s°. 
Oblate Spheroid. Page 161. The expressions for F,, G,, G, are same as given 
above, while F, is 


2 
i O8Sy (i — ys)2 (let 2 s)l ot : 


fal See 
4 cosh? a» st ) 


The values of the functions can be calculated as shown in the paper. The general 
trend of curves in Figures 2, 3, and 4 remains same. The conclusions of the paper 
remain unaltered. 
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Microwave Spectroscopy. By W. Gorpy, W. V. Situ, and R. F. TRam- 
BARULO (John Wiley & Sons, Inc., New York 1953). 446 pp. with more than 90 
Figures and Tables; $8.-. 

Die Mikrowellenspektroskopie ist in der kurzen Zeitspanne ihres Bestehens zu 
einem der reizvollsten Gebiete der Kernphysik, Hochfrequenztechnik und physi- 
kalischen Chemie geworden. Die stiirmische Entwicklung besonders der letzten 
fiinf Jahre liess oft den Mangel eines Standardwerkes fiihlbar werden. Das vor- 
liegende Werk bekannter amerikanischer Mikrowellenspektroskopiker schliesst 
diese Liicke und bietet gleichzeitig auch in bezug auf theoretische Zusammen- 
hinge eine detaillierte Darstellung des gegenwartigen Standes der Mikrowellen- 
spektroskopie und einiger verwandter Gebiete der Radiospektroskopie. 

Das erste Kapitel bringt eine gedrangte Zusammenstellung der technischen 
Hilfsmittel der Mikrowellenspektroskopie, wobei naturgemass die experimentelle 
Methodik der Autoren im Vordergrund steht. Von besonderem Wert ist dabei ein 
Abschnitt iiber die theoretische Empfindlichkeit verschiedener Spektrographen- 
typen. Die folgenden Kapitel tiber die Theorie der Mikrowellenspektren von 
Gasen, des Stark- und Zeeman-Effekts und der Intensitat der Rotationslinien 
reprasentieren eine kondensierte Zusammenfassung von in der Originalliteratur 
weitverstreuten, theoretischen und experimentellen Arbeiten. Der praktisch ar- 
peitende Mikrowellenspektroskopiker vermisst hier vielleicht die explizite Wie- 
dergabe der King-, Hainer-, CroBschen-Tabellen der Energieniveaus des asym- 
metrischen starren Kreisels. 


a a , Oo", hee) aac 2 
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Den Interessen des Festk6rperphysikers und des Kernphysikers kommen die 
Kapitel iiber para- und ferromagnetische Resonanz in Festkérpern bzw. tiber 
Kernquadrupolspektren entgegen. Fiir den physikalischen Chemiker enthalten 
die letzten Kapitel eine Fiille von Angaben tiber den Zusammenhang der Kopp- 


lung von Kernquadrupolmomenten und molekularen Dipol- und Quadrupolmo- — 
menten mit Zustanden der Valenzelektronen und modernen Problemen der — 


chemischen Bindung. 

Besonders wertvoll sind auch die als Anhang wiedergegebenen zahlreichen 
Tabellen experimenteller Resultate und die umfangreiche Zusammenstellung von 
Literaturangaben. 

Zweifellos wird das Buch von GorDy, SMITH und TRAMBARULO dem Mikro- 
wellenspektroskopiker zur Einfiihrung und als Nachschlagewerk von grésstem 

Nutzen sein. Hs. H. Giinthard 


Hohere Mathematik, Teil VI: Integration und Reihenentwicklung im Kom- 
plexen. Gewohnliche und partielle Differentialgleichungen. Von R. RoTHE und 
I. Szané (Verlagsgesellschaft B. G. Teubner, Leipzig und Stuttgart 1953). 251 S., 
34 Abb.; DM 17.60. 

Zu den bewahrten Leitfaden Héhere Mathematik von R. RoTHE ist nun eine 
im Grunde schon langst notwendige Erganzung erschienen. Dieser von I. SzaBé 
verfasste Teil beginnt mit ausgewahlten Fragen aus der Funktionentheorie (Inte- 
gration, Reihenentwicklung, asymptotische Reihen). Dann wird die allgemeine 
Theorie der gewohnlichen linearen Differentialgleichungen weiter ausgebaut und 
das Verhalten der Lésungen in der Umgebung isolierter Singularitaten unter- 
sucht (Fuchssche Theorie). Ein grosser Teil des Bandes ist den wichtigsten spe- 
ziellen linearen Differentialgleichungen und den Funktionen der mathematischen 
Physik gewidmet. Den Abschluss bilden einige grundsdtzliche Bemerkungen zu 


| 


; 


den partiellen Differentialgleichungen erster und zweiter Ordnung (diese sollen — 


im noch folgenden Teil VII ausfiihrlich behandelt werden). 

Die Auswahl des Stoffes ist sehr gliicklich getroffen; die Darstellung klar und 
mathematisch einwandfrei. Jeder Abschnitt enthalt eine Anzahl Ubungsaufgaben, 
die zum Teil auch die Theorie weiterfiihren. Das handliche Buch wird den Stu- 
dierenden der angewandten Mathematik, der Physik und der Ingenieurwissen- 
schaften sowie den mathematisch interessierten Physikern und Ingenieuren sehr 
gute Dienste leisten und kann daher vorbehaltlos empfohlen werden. 

E, Roth-Desmeules 


Small Particle Statistics. By G. HERDAN; with a Guide to the Experi- 
mental Design of Particle Size-Determination, by M. L, SmrtH (Elsevier Publishing 
Co., Amsterdam, 1953). 520 pp.; 70s. 

Obwohl in jiingster Zeit umfassende Fach- und Lehrbiicher iiber moderne 
Statistik und ihre verschiedenen technischen Anwendungsgebiete, vor allem in 
der amerikanischen Literatur, erschienen sind, fehlt es an einer einschlagigen zu- 
sammenfassenden Darstellung der statistischen Verfahren im Gebiete der kleinen 
Stoffteilchen von siebbarer und untersiebbarer Grésse (von 1cm bis 10-5 cm 
Durchmesser), das heisst fiir einen Teilchengréssenbereich, der sich an die kol- 
loidalen Dimensionen unmittelbar anschliesst und fiir welche die Gesetze der phy- 
sikalischen Statistik nur bedingt iibertragbar sind. Es musste dies als ein emp- 
findlicher Mangel angesehen werden, indem der Technologie kérniger Stoffe eine 
beachtliche Bedeutung zukommt und ihre Weiterentwicklung ein verfeinertes 
Verstandnis der Verfahren zur Bestimmung von Korngrésse und ihrer Wechsel- 
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beziehungen (Verteilungsgesetz usw.) verlangt. Deshalb ist es ein grosses Ver- 
dienst des Verfassers, das umfangreiche Material, das in einer ansehnlichen Zahl 
von Abhandlungen in verschiedenen Zeitschriften zerstreut liegt, iiberarbeitet 

_ und in einem Gesamtwerk zusammengefasst zu haben. 

Das Buch zerfallt in die folgenden Hauptabschnitte: I. The Determination of 
the Fineness of Solids; II. The Technological Importance of the Fineness of 
Solids; III. Attainment of a Specified Fineness and of Homogeneity in a Sub- 
‘stance; IV. Principles of Statistic of Coarse Disperse Systems; V. Experimental 

Design and Experimental Errors of Particle Size Determination (verfasst von 
M. L. SuiryH). 

: Vom Leser werden keine naheren Vorkenntnisse iiber physikalische Statistik 
verlangt, indem die Statistik der siebbaren und untersiebbaren Stoffteilchen 
unabhangig als neues Erkenntnisgut in unmittelbarem Kontakt mit den experi- 

-mentellen Beobachtungen und Erfordernissen aufgebaut wird. Eine wertvolle 
Bereicherung bildet der letzte Abschnitt, der von einem erfahrenen Mitarbeiter 
verfasst wurde und eine kritische Wiirdigung der heute gebrauchlichen Verfahren 
der Teilchengréssenbestimmung darstellt. Das Buch wird sich als willkommenes 

_Hilfsmittel erweisen fiir alle Industrie- und Hochschullaboratorien, die sich mit 
dem leidigen, aber sehr interessanten Gebiet der dispersen Systeme tiberkolloi- 
daler Teilchendimensionen auseinandersetzen miissen. R. Sdnger 


. 
, 


Introduction to Solid State Physics. By CuarLes KitTeLt (John Wi- 
ley & Sons, Inc., New York 1953). XII + 396 pp. 

Der Autor schreibt diese Ver6ffentlichung fiir Kreise, die sich einen ersten 
und doch nicht allzu summarischen Uberblick iiber das Gebiet « Grundfragen der 
Physik des festen Kérpers» erwerben wollen. Das Ergebnis leistet aber mehr: die 
Fiille des gebotenen Stoffes, alles in praziser und knapper Form, stempelt das 
Werk zu einem handlichen Vademekum fiir Leser, die verschiedene Fragen aus 
der angewandten Physik des festen Kérpers bearbeiten miissen und nicht ohne 
ein gewisses Mass an tieferer Einsicht in den grundsatzlichen Sachverhalt bleiben 
wollen. Eine gute Ausniitzung des Buches fordert allerdings Mittel, wie sie etwa 
im bekannten Joos, Theoretische Physik, zur Darstellung gelangen, und das ist 
erheblich mehr als die im Vorwort etwas optimistisch veranschlagten Minimal- 
kenntnisse. 

Speziell bemerkenswert ist die Aufnahme neuester Messmethoden und Mess- 
ergebnisse sowie die den einzelnen Kapiteln beigegebene Aufgabensammlung und 
Literaturliste von buch- und referatmdssigen Publikationen tiber die beziiglichen 
Sachgebiete. 

Die Kapitel 1 bis 5 handeln uber strukturelle und thermische, die Kapitel 6 
bis 10 iiber dielektrische und magnetische Eigenschaften von Kristallen. Elek- 
trische Leitfahigkeit wird in den Kapiteln 11 bis 14, der Begriff «Stérstelle» in 
den Kapiteln 15 und 16 diskutiert. Ein reichhaltiger Appendix erganzt in stoff- 
licher und didaktischer Hinsicht die einzelnen Abschnitte. W. Baumgartner 


Einfiihrung in die Mikrowellen und ihre wissenschaftlichen Anwen - 
dungen. Von H. H. KiincGeER (Verlag S. Hirzel, Stuttgart 1954). 117 S:, 92 Abb.; 
DM 9.60. 

Die Erschliessung eines neuen Frequenzgebietes elektromagnetischer Wellen 
bedeutet zugleich den Beginn einer Bereicherung der Wissenschaft durch neue 
Méglichkeiten und Erkenntnisse. Diese Tatsache hat sich in héchstem Masse bei 
den Mikrowellen (Zentimeterwellen, Millimeterwellen) bewahrheitet, wobe1 die 
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Leistungsfahigkeit und Organisation der modernen Forschung das Ihrige ail 
einer beinahe sprunghaft raschen Entwicklung beitrug. H. H. KLINGER hat sich — 
in verdienstvoller Weise in seinem Buch zum Ziel gesetzt, einen weiten Kreis von — 
Lesern — Studierende, Physiker und Naturwissenschafter, welche nicht Hoch- — 
frequenzfachleute sind — mit dem reichen Tatsachenmaterial in einer einfachen, ~ 
2 _ tibersichtlichen und klaren Form vertraut zu machen. Diese Zielsetzung, verbun-— 
Exe den mit der notwendigen Beschrankung des Gesamtumfanges, gestattet kein ver-_ 
an tieftes Verweilen bei Einzelheiten; der interessierte Leser findet aber wertvolle 
_  Hinweise in dem ausfiihrlichen Literaturverzeichnis. Der Stoff ist reichhaltig und — 
i ‘zweckmassig gegliedert: Hohlleiter, Hohlraumresonatoren, Erzeugung von Mikro- 
Pea’ wellen, Detektoren, Messtechnik, Mikrowellenspektroskopie, Mikrowellen-Reso- 
nanz-Absorption in Kristallen, anomale Dispersion .und Absorption von Mikro- 
. wellen in polaren Fliissigkeiten, elektrische Leitfahigkeit von Metallen bei hohen 
_... Frequenzen, biologische und medizinische Wirkungen, Mikrowellen in der Astro- 
-nomie und Astrophysik, Mikrowellen in der Experimentalphysik, Mikrowellen in 
- der Kernphysik. : 
Das kleine Werk liest sich gut und kann zur allgemeinen Orientierung oder 
zur Einfiihrung in das Gebiet der wissenschaftlichen Anwendungen der Mikro- 
_ wellen bestens empfohlen werden. F. Tank 


St ee Electro chemical Data. By B. E. Conway (Elsevier Publishing Company 
oe Amsterdam, 1952), 374°pp. 3) 55's. 
: ‘ Das Buch gibt in Tabellenform einen Uberblick tiber die wichtigsten elektro- 
chemischen und physikalisch-chemischen Daten unter Beriicksichtigung der 
a. neuesten Methoden und Ergebnisse. Der Stoff ist in folgende Kapitel eingeteilt: 
a I. Universelle Konstanten und Umrechnungsfaktoren. — II. Physikalische Da- 
| ten, wie Entropiewerte von Ionen in wassriger Lésung, Polarisierbarkeit von 
Ionen und Molekiilen, Ionenradien und kovalente Abstande. — III. Wechsel-_ 
ie , wirkungskrafte bei Ionen und Molekiilen in wassriger Lésung, wie Aktivitaten, 
pe Hydratationskonstanten, Aussalzkonstanten usw. — IV. Transport in Lésungen 
; starker Elektrolyte, wie Leitfahigkeiten und Ionenbeweglichkeiten verschiedener _ 
Salze in verschiedenen Lésungsmitteln. — V. Dissoziationskonstanten, Léslich- 
keiten und Pufferlésungen. — VI. Eigenschaften elektrischer Doppelschichten an 
Grenzflachen. — VII. Allgemeine Eigenschaften kolloidaler und makromoleku- 
larer Elektrolyte von biologischem Interesse. — VIII. Elektrochemie von Schmel- 
zen bei hohen Temperaturen. — IX. Reversible Elektrodenprozesse. — X. Kinetik 
elektrolytischer Vorgange. 
Ausfihrlicher Quellennachweis ermoglicht zusatzliche Orientierung in der 
einschlagigen Literatur. Das Tabellenwerk ist fiir Chemiker und Physikochemiker 
sehr wertvoll, und die gedrangte, klare Darstellung erméglicht ein rasches Auf- 
a finden der gewiinschten Daten. Allen Laboratorien, die auf elektrochemischem 
Gebiet tatig sind, kann das Werk empfohlen werden. F. Held 


Tables des fonctions de Legendre associées, calculées pour le centre 
national d’études des télécommunications, (Edition de la Revue d’optique, 
Paris 1952). 291 pp. 

Diese Tafeln geben die zugeordneten Legendreschen Funktionen P™(cos @) 
fir n= —0,5, (0,1), 10, fiir m= 0, (1), 5 und fiir 0 = 0, (1°), 90° mit durch- 
schnittlich fiinf bedeutsamen Stellen und ohne Differenzen: sie bilden ein wert- 
volles Hilfsmittel fiir die numerische Lésung von Potentialaufgaben in sphari- 
schen Koordinaten. E. Stiefel 
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Neuerscheinung aus dem Birkhiuser Verlag 


WILLIAM PRAGER ri he 
PROBLEME 
DER PLASTIZITATSTHEORIE 


2 Von Professor Dr. WILLIAM PRAGER, i 
Brown University, Providence, R. I., USA. i 


100 Seiten mit 52 Figuren, Ganzleinen Fr. 12.50 (1955) 
(Sammlung «Lehr- und Handbiicher der Ingenieurwissenschaften», Band 17) 


INHALT 


I. Mechanisches Verhalten plastischer Stoffe. II. Mechanisches Verhalten 
plastischer Tragwerke. III. Traglastverfahren. IV. Endliche plastische 
Formanderungen. 


Das Werk stellt im wesentlichen eine Niederschrift von Gastvorlesungen 
dar, die der Verfasser im November und Dezember 1954 an der Eidgendés- Nora! 
sischen Technischen Hochschule in Ziirich hielt. Die Beschrankung auf 
ausgewahlte Probleme der Theorie ideal-plastischer Stoffe ermdglicht es 
dem Verfasser, den Leser bisan die Front gegenwartiger Forschung zu fiihren, 
ohne vorherige Bekanntschaft mit der Plastizitatstheorie vorauszusetzen. 
Die Mehrzahl dieser Probleme entstammen der Fachliteratur der letzten 
Jahre und sind im vorliegenden Band wohl zum ersten Male in Buchform 
behandelt. Die Darstellung zielt hauptsachlich darauf ab, im Leser ein 
intuitives Verstandnis fiir die Besonderheiten des plastischen Verhaltens 
zu erwecken. Zahlreiche Literaturverweise geben dem Leser die Méglich- 
keit, dieses Verstandnis durch weiteres Studium zu vertiefen. 
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